def test_autobin(self):
"""The 'autobin' command."""
bam_fname = "formats/na12878-chrM-Y-trunc.bam"
target_bed = "formats/my-targets.bed"
targets = tabio.read(target_bed, 'bed')
access_bed = "../data/access-5k-mappable.hg19.bed"
accessible = tabio.read(access_bed, 'bed').filter(chromosome='chrY')
for method in ('amplicon', 'wgs', 'hybrid'):
(cov, bs), _ = autobin.do_autobin(bam_fname, method,
targets=targets,
access=accessible)
self.assertGreater(cov, 0)
self.assertGreater(bs, 0)
def test_autobin(self):
"""The 'autobin' command."""
bam_fname = "formats/na12878-chrM-Y-trunc.bam"
target_bed = "formats/my-targets.bed"
targets = tabio.read(target_bed, 'bed')
access_bed = "../data/access-5k-mappable.hg19.bed"
accessible = tabio.read(access_bed, 'bed').filter(chromosome='chrY')
for method in ('amplicon', 'wgs', 'hybrid'):
(cov, bs), _ = autobin.do_autobin(bam_fname, method,
targets=targets,
access=accessible)
self.assertGreater(cov, 0)
self.assertGreater(bs, 0)