New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4035360	NA (NA)	8 (0.260)	8/234	2.4	100%	intergenic (+207/‑171)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon
	?	NC_000913	= 4035384	0 (0.000)					intergenic (+231/‑147)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon

ACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAA  >  NC_000913/4035352‑4035525                                 |                                                                                                                                              acacCCCGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtt                                         >  1:104784‑M2/1‑33 (MQ=255) acacCCCGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtt                                         <  2:104784‑M2/135‑103 (MQ=255)     cccGCGCCGCTGAGAAAAAGCGAAGCGGCacggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgc                                                                                          >  3:10024‑M2/1‑29 (MQ=255)     cccGCGCCGCTGAGAAAAAGCGAAGCGGCacggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgc                                                                                          <  4:10024‑M2/82‑54 (MQ=255)        gcgcCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttg                                 >  7:27633‑M2/1‑26 (MQ=255)        gcgcCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttg                                 <  8:27633‑M2/136‑111 (MQ=255)          gccgcTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttattttt                                                                                            >  1:125420‑M2/1‑24 (MQ=255)          gccgcTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttattttt                                                                                            <  2:125420‑M2/75‑52 (MQ=255)                            agcgGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggag  >  4:7019‑M2/1‑6 (MQ=255)                             gcggcgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgt                                   >  7:141233‑M2/1‑5 (MQ=255)                             gcggcgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgt                                   <  8:141233‑M2/113‑109 (MQ=255)                                 |                                                                                                                                              ACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAA  >  NC_000913/4035352‑4035525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ 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.