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

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

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

the documentation in
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 )

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
Senior Member
Location: Germany

Join Date: May 2010
Posts: 101

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

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