Hey Huys,
How can I generate an alignment statistics report from bwa or samtools? I used bwa for alignment and samtools for the rest. Any other tool to be used once I generate the bam file? I basically want fields like below generated from a commercial software.
MeanInsertSize
Total read#
Uniquely paired read#
Non-uniquely paired read#
Uniquely mapped read1#
Non-uniquely mapped read1#
Uniquely mapped read2#
Non-uniquely mapped read2#
Unmapped read1#
Unmapped read2#
Uniquely paired read%
Non-uniquely paired read%
Uniquely mapped read1%
Non-uniquely mapped read1%
Uniquely mapped read2%
Non-uniquely mapped read2%
Unmapped read1%
Unmapped read2%
2 answers
I think what you want can be achieved using Picard, which is an implementation of Samtools, but in Java. There are also a few tools in Picard which are absent in Samtools. Namely, you want to use CollectAlignmentSummaryMetrics. http://picard.sourceforge.net/javadoc/net/sf/picard/analysis/CollectAlignmentSummaryMetrics.html
You can run this on your SAM/BAM, and you will have most of what you describe above.
EDIT: If you have memory issues, you can increase / decrease the amount of memory available to the program using these flags: "-vmargs -Xms###[mg] -Xmx###[mg]"
java jvm-args -jar PicardCommand.jar
becomes
java -vmargs -Xms64m -Xmx4000m -jar PicardCommand.jar
The first Xms is the minimum amount of memory available to the program when it starts, and the second Xmx is the max memory limit. You can specify the amount of memory using either megabytes (m) or gigabytes (g).
Please see new post as I could not fit everything as a comment to this post:
Edited my original answer to include information about increasing initial memory available to a java app.
Log in to answer this question.
-
Why Is My Uniquely Mapped Read Count Higher in Samtools Compared to Bowtie2?
written by daffodil •
-
Total number of trimmed reads less than number of mapped reads
written by kspata
-
Getting uniquely aligned properly paired reads
written by nadiabeg.comsats •
-
samtools stats result different from samtools flagstat
written by snishtala03 •
-
Realign Aligned.sortedByCoord.out.bam to Aligned.toTranscriptome.out.bam for RSEM input
written by mkvalluru •
-
difference between mapq=0 and samtools flag 0x900
written by bitpir
-
Why The Properly Paired Results From Bwa Mem Is An Odd Number
written by liupfskygre
-
Insert Size For Illumina Gaiix Paired-End Library From Sam/Bam File
written by bioinfo
-
Output From Picard And Samtools Flagstat Not Very Useful From My Sambam Files:
written by Angel
-
Should Samtools Pileup Be Performed On Uniquely Mapped Reads Or All The Reads?
written by Doctoroots
http://biostar.stackexchange.com/questions/9552/basic-bam-file-annotation