Entering edit mode
2.2 years ago
Berk
•
0
Sometimes there is more than one gene in the same region, for example paralog genes. I want to find the reads that are multi-mapped in the same region. I have a genome aligned BAM file as input and score the reads to the genes with pysam.