SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools flagstat output nguyendofx Bioinformatics 23 01-22-2014 04:45 AM
Samtools flagstat - no duplicates? Orr Shomroni Bioinformatics 3 11-25-2011 12:46 AM
Samtools flagstat Anelda Bioinformatics 0 09-26-2011 03:55 AM
mapped vs. aligned in TopHat and flagstat PFS Bioinformatics 2 04-07-2011 09:34 AM
samtools flagstat bair Bioinformatics 3 05-28-2010 06:15 AM

Reply
 
Thread Tools
Old 08-19-2011, 05:01 AM   #1
fabrice
Member
 
Location: paris

Join Date: Oct 2009
Posts: 86
Default samtools flagstat and ViewSam ALIGNMENT_STATUS=Aligned

I am using and samtools flagstat to statistc how many reads mapped in my s_4.bam (produced by eland, convert from export files).

samtools flagstat s_4.bam

Then I used picard.

java -jar ViewSam.jar INPUT=s_4.bam ALIGNMENT_STATUS=Aligned PF_STATUS=All VALIDATION_STRINGENCY=SILENT > s_4_Aligned_reads.sam

I expected that the line number of s_4_Aligned_reads.sam will be equal to mapped number report by samtools. But it is biger than the number report by samtools.

Does any one know why?

samtools flagstat s_4.bam

43055574 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 duplicates
8834640 + 0 mapped (20.52%:nan%)
43055574 + 0 paired in sequencing
21527787 + 0 read1
21527787 + 0 read2
0 + 0 properly paired (0.00%:nan%)
5208944 + 0 with itself and mate mapped
3625696 + 0 singletons (8.42%:nan%)
117012 + 0 with mate mapped to a different chr
112688 + 0 with mate mapped to a different chr (mapQ>=5)

wc -l s_4_Aligned_reads.sam
11789604 s_4_Aligned_reads.sam



Thanks.

Last edited by fabrice; 08-19-2011 at 05:11 AM.
fabrice 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 12:30 PM.


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