New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3426897 =	NA (NA)	17 (0.120)	7/168	NT	NA	intergenic (‑114/+61)	rrlD/alaU	23S ribosomal RNA of rrnD operon/tRNA‑Ala
	?	NC_000913	3426933 =	NA (NA)					intergenic (‑150/+25)	rrlD/alaU	23S ribosomal RNA of rrnD operon/tRNA‑Ala

TGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGAAGTATTTTTTATTTAATCACTACAGAGATGGTGGAGCTATGCGGGATCGAACCGCAGACCTCCTGCGTGCAAAGCAGGCGCTCTCCC  >  NC_000913/3426873‑3427016                                                             |                                                                                    tgctCCAGATTTTTAAAGAAAAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGAAGTATTTTTTATTTAATCACTACAGAGATGGTTGAGCTTTGCg                                              >  2:1866999/4‑100 (MQ=1)                                             agttatacTCTGAAGTATTTTTTATTTAATCACTACAGAGATGGTGGAGCTATGCGGGATCGAACCGCAGACCTCCTGCGTGCAAAGCAGGCGCTCTccc  <  2:736338/94‑1 (MQ=0)                                                             |                                                                                    TGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAACCTTTGCAGGTACACTCTGAAGTATTTTTTATTTAATCACTACAGAGATGGTGGAGCTATGCGGGATCGAACCGCAGACCTCCTGCGTGCAAAGCAGGCGCTCTCCC  >  NC_000913/3426873‑3427016

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.