SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why MAQ consensus seq better than SAMtools consensus ?? av_d Genomic Resequencing 5 10-18-2015 04:44 AM
Problem calling consensus sequence with Samtools mpileup? Lspoor Bioinformatics 3 10-24-2011 08:43 AM
meaning of every column of accepted_hit.bam from tophat chenyao Bioinformatics 7 08-12-2011 09:26 AM
MACS SAMfiles khb Bioinformatics 2 12-19-2010 10:37 AM
question on SAMtools consensus calling orionzhou Bioinformatics 9 11-16-2010 03:42 PM

Reply
 
Thread Tools
Old 07-20-2010, 03:39 AM   #1
KevinLam
Senior Member
 
Location: SEA

Join Date: Nov 2009
Posts: 203
Default Column order in samfiles in consensus calling

the documentation in
http://samtools.sourceforge.net/cns0.shtml
is slightly confusing at this portion

When consensus calling is switched on, pileup command will insert the consensus base, consensus quality, SNP quality and maximum mapping quality of the reads covering the sites between the `reference base' and the `read bases' columns.

as 'the number of reads covering the site' is actually between the `reference base' and the `read bases' columns. (see http://samtools.sourceforge.net/pileup.shtml )

I believe the order of the columns should be as per this csv line?

"chromosome"," 1-based coordinate"," reference base","consensus base"," consensus quality"," SNP quality ","maximum mapping quality "," the number of reads covering the site"," read bases"," base qualities"

Last edited by KevinLam; 07-20-2010 at 03:50 AM.
KevinLam is offline   Reply With Quote
Old 07-28-2010, 10:05 AM   #2
epigen
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 101
Default

You're right. It seems that the docu is outdated. Several people have also noted the presence of two additional columns for indel lines in pileup that are not explained in the manual or FAQs.
epigen 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 02:46 PM.


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