New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4035360	NA (NA)	10 (0.010)	7/266	NT	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

TTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAA  >  NC_000913/4035346‑4035509               |                                                                                                                                                      ttatGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCacggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgt                <  2:1625367‑M2/150‑112 (MQ=255)       acacCCCGCGCCGCTGAGAAAAAGCGAAGCGGgtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtcagt           <  1:11141270‑M2/150‑118 (MQ=255)       acacCCCGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg          <  1:2622044‑M2/150‑118 (MQ=255)       acacCCCGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg          <  2:1669569‑M2/150‑118 (MQ=255)         acCCCGCGCCGCTGAGAAAAAGCGAAGCGGCacggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtc        <  1:8489483‑M2/150‑120 (MQ=255)            ccGCGCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtcagt     <  2:4105668‑M2/150‑123 (MQ=255)             cgcgCCGCTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtcagtg    <  1:1211599‑M2/150‑124 (MQ=255)              gcgcCGCTGAGAAAAAGCGAAGCGGCacggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtcagtgg   <  2:2721077‑M2/150‑125 (MQ=255)              gcgcAGCTGAGAAAAAGCGAAGCGGCacggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtcagtgg   <  2:8867956‑M2/150‑125 (MQ=255)               cgccgcTGAGAAAAAGCGAAGCGGCgcggggtgtgcataatacgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgtcagtggt  <  1:12238317‑M2/150‑126 (MQ=255)               |                                                                                                                                                      TTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAA  >  NC_000913/4035346‑4035509

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