How can I find reads that have multi-mapped with multiple genes overlapping in the same region
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.
• 774 views
•
link
0 answers
No answers yet.
Log in to answer this question.