def test_alignment():
s = BAM(datatest)
# no need to call reset but does not harm and reminds us that it shoudl be
# used in general to make sure we start at the beginning of the iterator.
s.reset()
a = Alignment(next(s))
a.as_dict()
def test_alignment():
datatest = sequana_data("test.bam", "testing")
s = BAM(datatest)
# no need to call reset but does not harm and reminds us that it shoudl be
# used in general to make sure we start at the beginning of the iterator.
s.reset()
a = Alignment(next(s))
a.as_dict()
def test_bam(tmpdir):
s = BAM(datatest)
assert len(s) == 1000
assert s.is_sorted is True
assert len(list(s.iter_unmapped_reads())) == 2
s.reset()
assert len(list(s.iter_mapped_reads())) == 998
s.reset()
# call this here before other computations on purpose
with TempFile(suffix=".json") as fh:
s.bam_analysis_to_json(fh.name)
assert s.get_read_names()
s.get_mapped_read_length()
s.get_stats()
s.get_full_stats_as_df()
with TempFile(suffix='.png') as fh:
s.plot_bar_flags(filename=fh.name, logy=True)
s.plot_bar_flags(filename=fh.name)
with TempFile(suffix='.png') as fh:
s.plot_bar_mapq(filename=fh.name)
with TempFile() as fh:
s.to_fastq(fh.name)
from sequana import FastQ
ff = FastQ(fh.name)
len(ff) == len(s)
s.get_gc_content()
s.get_length_count()
s.plot_gc_content()
try:
s.plot_gc_content(bins=[1, 2, 10])
assert False
except:
assert True
