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.
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy