Closed:Difference between samtools unmapped reads and picard metrics file unmapped reads
1
0
Entering edit mode
6.4 years ago
DVA ▴ 630

I have noticed that the count of unmapped reads from:

samtools view -f 4 $sample.sam > $sample.unmapped.sam

is LESS than (about ~90% of) what I got from the "UNMAPPED_READS" column in the metrics file generated by MarkDuplicates via picard:

java -Xmx45g -jar /dir/picard.jar MarkDuplicates INPUT=$sample.sam OUTPUT=$sample.sort.mdup.bam METRICS_FILE=$sample.metrics.mdup.txt

I checked picard definition of UNMAPPED_READS, and it is: "The total number of unmapped reads examined. (Primary, non-supplemental)". I don't understand why the two numbers do not match. Any suggestion please? Thank you.

samtools picard • 119 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2861 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6