Extracting percent identity of mapped reads and the coordinates on the genome the read mapped from a BAM file?

1

Hello,

I am working with bam files that have environmental reads mapped back to an isolate genome. I am looking for a way to extract both the percent identity of all mapped reads in the BAM file and the coordinates on the genome where they mapped. I was looking at some combination of maybe bedtools and pysam or parsing the CIGAR string but nothing I have tried has worked so far.

I was trying to find a pythonic answer, but if people have other solutions feel free to share.

Thanks!


Bedtools


PySam


python


BAM


Samtools

• 60 views



Source link