SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
No output from Samtools pileup pravee1216 Bioinformatics 1 12-05-2012 03:59 AM
output info from samtools pileup? shuang Bioinformatics 1 09-26-2011 04:35 PM
SAMtools pileup output zippered_ohio Bioinformatics 1 09-08-2011 01:27 PM
samtools pileup output file charlie_sequencing Bioinformatics 1 02-10-2011 12:51 PM
samtools pileup output for SNP detection Mansequencer Bioinformatics 1 09-08-2010 06:36 PM

Reply
 
Thread Tools
Old 03-30-2011, 12:03 PM   #1
mr_boourns
Junior Member
 
Location: New Haven

Join Date: Apr 2010
Posts: 2
Default pileup output different using Maq and Samtools commands

Dear all,

We are trying to detect allelic specific expression with RNA-seq data by mapping against a masked genome (with the reference allele substituted for N's at known SNP positions). We have followed the approach of the Degner/Pritchard paper (2009) where they perform the mapping step using Maq. Our approach has been to use the pileup command in Maq to derive the allele counts at each position. We thought we were getting good results. However, we recently realized that the counts from the Maq produced pileup file do not agree with the counts from a visual inspection of the BAM file produced from the mapping. (This BAM file was produced by converting the .map file from Maq to a .bam file using Samtools). However, when we use the Samtools pileup command on this BAM file, the counts in the pileup file do agree with a visual inspection of the counts when viewing the BAM file in a genome browser.

For instance, for the Maq produced pileup file we get 24 T's and 0 C's at a particular base, whereas in the Samtools produced pileup file we get
14C's and 35T's.

Does anyone have any idea of what is happening here? Should we trust the Maq or the SAM tools results?

Thanks so much for the help,

John
mr_boourns is offline   Reply With Quote
Old 03-31-2011, 03:05 AM   #2
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 415
Default

I would use Samtools as it is under active development and Maq, to my knowledge, is not any more.

However for SNP calling try the mpileup pipeline described on the Samtools webpage instead of using pileup. It's important to generate the I16 info fields for variant calling in our experience.
colindaven is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:52 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO