Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • mpileup and pileup results

    Hi everyone.

    1. I tried to call the snps using both mpileup and pileup.
    There is more differences in snps calling numbers.
    I found that mpileup give me less snps compared too pileup which produce lots of snps calling.
    eg pileup = about 10 K snps while mpileup = about 800 snps

    i run using this command
    samtools mpileup -uf ref.fa aln1.bam aln2.bam | bcftools view -bvcg - > var.raw.bcf

    is using mpileup is more thorough compared to pileup?

    btw, i call the snps from rice data.

    2. can anyone here explain how can i make use the genotype data that present in vcf 'FORMAT' ? I read already in the documentation but not really understand.

    Thanks all!

  • #2
    According to the documentation from Samtools pileup is now depreciated. Use mpileup.

    Apparently the algorithm is more complex, takes into account local realignment of reads around putative SNP positions and also can modify base qualities around these positions. See the part on BAQ in the mpileup docs.
    Probably mpileup is reducing the number of false positives here, but have a look at some of the data. I'm sure you'll see some patterns emerging.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin




      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
      04-22-2024, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 11:49 AM
    0 responses
    15 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-24-2024, 08:47 AM
    0 responses
    16 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    61 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    60 views
    0 likes
    Last Post seqadmin  
    Working...
    X