samtools view -X <bam file>
Convert bitwise FLAG to human readable FLAG.
samtools flagstat <bam file>
Print number of mapped reads and so on.
samtools view -H <bam file>
check if the bam files is sorted.
http://biostar.stackexchange.com/questions/5273/is-my-bam-file-sorted
Manual Reference Pages - samtools
how to extract unique hits from bam files
samtools view -q 1 <bam file>
You can use samtools view -q 1 to only output uniquely mapped reads
(and omit both unmapped and ambiguous ones). Ambiguous reads that can
be mapped to different locations with the same score have mapping
quality 0 in BWA.
samtools view bwa.bam | grep "XT:A:U"
BWA outputs have XT tag set to "U".
Hide Comments