Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Impact on quality of SNP calls using samtools mpileup

    I have a question on the impact on quality of SNP calls using samtools mpileup: calling each sample individually vs calling all samples altogether
    I basically follow the samtools URL:
    http://samtools.sourceforge.net/mpileup.shtml and
    http://samtools.sourceforge.net/samtools.shtml (except using -D 2000 or -d 10 for filtering)

    the commands I used basically like below:
    Call all samples together (each bam file is from exome-seq of one sample )
    samtools mpileup -ugf ref.fa aln1.bam aln2.bam aln3.bam ... *.bam| bcftools view -bvcg - > var.raw.bcf
    bcftools view var.raw.bcf | vcfutils.pl varFilter -d 10 > var.flt.vcf

    Call one samples a time (each bam file is from one sample)
    for loop (for each bam file x.bam)
    samtools mpileup -ugf ref.fa x.bam| bcftools view -bvcg - > var.raw.bcf
    bcftools view var.raw.bcf | vcfutils.pl varFilter -d 10 > var.flt.vcf
    end of loop

    We observed that if we use samtools to call each sample (bam file) individually one at a time or call all samples (multiple bam files) together, then we seem get a better SNP calls if we call individually than if we call them altogether (of course we get one vcf file for each sample/bam file if call individually, whereas just one single vcf file for all SNPs of all samples for calling them altogether. We do lose homo reference SNPs in vcf file if calling individually, since for single sample, we only see hetero or homo variant, no homo reference unless we call multiple sample altogether, then some samples would have homo reference calls for sites that other samples were called as hetero or homo variants)

    Does any one have similar observation or any comments or insight for such observation (get better SNP calls if we call individually than if call together)? Since I know most people would call all bam files altogether and would not think that calling each sample individually would get a better quality SNPs (using samtools mpileup SNP call). Any interpretation from algorithm side of samtools? We all thought it shall not matter, but we do have bench mark data to show that is true at least in our hand. We use -d 10 or -D 2000 in the filtering step for both cases.

  • #2
    Have you asked on the samtools mailing list?



    I'd be interested to hear the reply.
    Homepage: Dan Bolser
    MetaBase the database of biological databases.

    Comment


    • #3
      Thanks, Dan! yes, I did, but not much luck and input from there...

      Comment


      • #4
        mpileup vs single calling

        Did you obtaned any answer for this??
        I am wondering why I obtain a higher total amount of variants during single callings than if i do mpileup with samtools even filtering in the same way for both approaches

        for several files
        Originally posted by myi View Post
        I have a question on the impact on quality of SNP calls using samtools mpileup: calling each sample individually vs calling all samples altogether
        I basically follow the samtools URL:
        http://samtools.sourceforge.net/mpileup.shtml and
        http://samtools.sourceforge.net/samtools.shtml (except using -D 2000 or -d 10 for filtering)

        the commands I used basically like below:
        Call all samples together (each bam file is from exome-seq of one sample )
        samtools mpileup -ugf ref.fa aln1.bam aln2.bam aln3.bam ... *.bam| bcftools view -bvcg - > var.raw.bcf
        bcftools view var.raw.bcf | vcfutils.pl varFilter -d 10 > var.flt.vcf

        Call one samples a time (each bam file is from one sample)
        for loop (for each bam file x.bam)
        samtools mpileup -ugf ref.fa x.bam| bcftools view -bvcg - > var.raw.bcf
        bcftools view var.raw.bcf | vcfutils.pl varFilter -d 10 > var.flt.vcf
        end of loop

        We observed that if we use samtools to call each sample (bam file) individually one at a time or call all samples (multiple bam files) together, then we seem get a better SNP calls if we call individually than if we call them altogether (of course we get one vcf file for each sample/bam file if call individually, whereas just one single vcf file for all SNPs of all samples for calling them altogether. We do lose homo reference SNPs in vcf file if calling individually, since for single sample, we only see hetero or homo variant, no homo reference unless we call multiple sample altogether, then some samples would have homo reference calls for sites that other samples were called as hetero or homo variants)

        Does any one have similar observation or any comments or insight for such observation (get better SNP calls if we call individually than if call together)? Since I know most people would call all bam files altogether and would not think that calling each sample individually would get a better quality SNPs (using samtools mpileup SNP call). Any interpretation from algorithm side of samtools? We all thought it shall not matter, but we do have bench mark data to show that is true at least in our hand. We use -d 10 or -D 2000 in the filtering step for both cases.

        Comment

        Latest Articles

        Collapse

        • 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
        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 04-11-2024, 12:08 PM
        0 responses
        13 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        19 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 09:21 AM
        0 responses
        15 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-04-2024, 09:00 AM
        0 responses
        43 views
        0 likes
        Last Post seqadmin  
        Working...
        X