getMT: grab the mitochondrial reads from a BAM and estimate their...

Description Usage Arguments

Description

grab the mitochondrial reads from a BAM and estimate their fraction nb. this could probably be done faster for a list of BAMs but it's not nb. nb. returns NuMt-depleted mitochondrial GenomicAlignments

Usage

1
getMT(bam, chrM = "chrM", mtGenome = "hg19", plotMAPQ = FALSE)

Arguments

bam

a BAM filename, which should have an index (else index it)

chrM

what the mitochondrial contig is called. Default is "chrM"

mtGenome

what mitochondrial assembly was used (default is hg19)

plotMAPQ

plot distribution of mitochondrial mapping quality? (FALSE)


RamsinghLab/ATACseeker documentation built on May 8, 2019, 8:05 a.m.