New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4035245 =	NA (NA)	6 (0.060)	3/168	NT	100%	intergenic (+92/‑286)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon
	?	NC_000913	4035272 =	0 (0.000)					intergenic (+119/‑259)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon

CGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAA  >  NC_000913/4035173‑4035371                                                                                                    |                                                                                                    cgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg                                                                                                     <  2:967834‑M2/1‑1 (MQ=255) cgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg                                                                                                     <  1:2573930‑M2/1‑1 (MQ=255) cgccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtg                                                                                                     >  1:845904‑M2/100‑100 (MQ=255)  gccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagatcaccgctgacccggcggcgtgtttgccgttgttccgtgt                                                                                                    <  1:1722729‑M2/2‑1 (MQ=255)  gccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtgt                                                                                                    <  2:52443‑M2/2‑1 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   >  1:1036837‑M2/98‑100 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   <  1:1513129‑M2/3‑1 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   >  1:194046‑M2/98‑100 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   >  1:2277008‑M2/98‑100 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   <  1:227720‑M2/3‑1 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   <  2:2762388‑M2/3‑1 (MQ=255)   ccttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTc                                                                                                   <  1:2372969‑M2/3‑1 (MQ=255)    cttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCa                                                                                                  >  1:1183476‑M2/97‑100 (MQ=255)     ttctcgctacggagttatgcatttatttttgctttctctgcgggagttctgaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCAg                                                                                                 <  2:1811992‑M2/5‑1 (MQ=255)     ttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCAg                                                                                                 <  1:331408‑M2/5‑1 (MQ=255)     ttcccgctacagagtcaagcatttatttttgctttctctgccggagttctcaggagaaccccgctgacccggcggcgtgtttgccgttgttccgtGTCAg                                                                                                 <  1:383297‑M2/5‑1 (MQ=255)                                                              accccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAg                                        >  1:378892‑M2/39‑100 (MQ=255)                                                               ccccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGc                                       >  1:1840128‑M2/38‑100 (MQ=255)                                                                cccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGc                                       >  1:3031743‑M2/37‑99 (MQ=255)                                                                cccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCg                                      >  1:78325‑M2/37‑100 (MQ=255)                                                                cccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCg                                      >  2:2373539‑M2/37‑100 (MQ=255)                                                                cccgctgacccggcggcgtgtttgccgttgttccgtGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCg                                      >  1:1361442‑M2/37‑100 (MQ=255)                                                                                                    GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGaa  <  1:677359‑M2/100‑1 (MQ=255)                                                                                                    GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGaa  >  2:2220818‑M2/1‑100 (MQ=255)                                                                                                    GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGaa  >  2:263950‑M2/1‑100 (MQ=255)                                                                                                    GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGaa  >  1:2301632‑M2/1‑100 (MQ=255)                                                                                                    GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGaa  >  2:39002‑M2/1‑100 (MQ=255)                                                                                                    |                                                                                                    CGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAA  >  NC_000913/4035173‑4035371

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