This is a test version of Biostars. For the public version, visit https://www.biostars.org.
mapped reads on target from samtools and Ion Torrent

so basically I'm trying to replicate a pipeline to get metrics similar to Coverage Analysis Report plugin from Ion Torrent S5 and I think I'm missing something here because when I run this command

samtools view -F 0x40 -L target.bed final.bam | cut -f1 | sort | uniq | wc -l

I only get 2042179 reads

but when I go to the report than I get from the plugin I get this

enter image description here

the total assigned amplicon reads is what I'm trying to get.

basically I have a bed file with every amplicon but when I use it to get the reads from that target I get less reads than I should

samtools

0 answers

No answers yet.

Log in to answer this question.