New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 438152	18 (0.700)	4 (0.190)	4/106	1.7	22.0%	intergenic (‑17/+163)	yajO/dxs	2‑carboxybenzaldehyde reductase/1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring
	?	NC_000913	= 438293	14 (0.680)					intergenic (‑158/+22)	yajO/dxs	2‑carboxybenzaldehyde reductase/1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring

ACAGCCGAGGCAAAGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTGCTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/438076‑438152 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccactcctgctatGCTTAAGAAATTATTCATAGACTCTAAATAATTCGAGTTGCAGGAAGGCGGCAAACGAGTGA  <  NC_000913/438293‑438232                                                                                                                                             ACAGCCGAGGCAAAGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCAC                                                                         >  1:116386/1‑68     CCGAGGCAAAGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCT                                                                     <  2:116386/68‑1     CCGAGGCAAAGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCT                                                                     <  2:178787/68‑1      CGAGGCAAAGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTG                                                                    >  1:533786/1‑68       GAGGCAAAGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTGC                                                                   <  1:3858/68‑1       GAGGCAAAGTCGGGAAACGCGAAGGTCGGGTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTGC                                                                   <  2:1728/68‑1              AGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTGCTATGCTT                                                            >  2:765093/1‑68                      AACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTGCTATGCTTAAGAAATT                                                    >  2:914957/1‑68                                TCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTGCTATGCTTAAGAAATTATTCATAGAC                                          <  1:914960/68‑1                                                   TGTATTGCATGCTGCCACTCCTGCTATGCTTAAGAAATTATTCATAGACTCTAAATAATTCGAGTTGC                       <  1:765096/68‑1                                                                     cCCTGCTATGCTTAAGAAATTATTCATAGACTCTAAATAATTCGAGTTGCAGGAAGGCGGCAAACGAG     <  2:129276/67‑1                                                                      CCTGCTATGCTTAAGAAATTATTCATAGACTCTAAATAATTCGAGTTGCAGGAAGGCGGCAAACGAGT    >  2:533786/1‑68                                                                        TGCTATGCTTAAGAAATTATTCATAGACTCTAAATAATTCGAGTTGCAGGAAGGCGGCAAACGAGTGA  >  1:332570/1‑68                                                                                                                                             ACAGCCGAGGCAAAGTCGGGAAACGCGAAGGTCGGTTTTTCCTAAGGGGTTGTATTGCATGCTGCCACTCCTGCTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/438076‑438152 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccactcctgctatGCTTAAGAAATTATTCATAGACTCTAAATAATTCGAGTTGCAGGAAGGCGGCAAACGAGTGA  <  NC_000913/438293‑438232

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