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

TCACTCGTTTGCCGCCTTCCTGCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATGCCAGCCAGGCCTTGATTTTGGCTTCCATACCAGCGGCATC  >  NC_000913/438232‑438359                                                              |                                                                   tCACTCGTTTGCCGCCTTCCTGCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCatagca                                                              >  1:332570‑M2/1‑62 (MQ=255)   aCTCGTTTGCCGCCTTCCTGCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCatagcagg                                                            >  2:533786‑M2/1‑60 (MQ=255)    cTCGTTTGCCGCCTTCCTGCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCatagcaggg                                                           <  2:129276‑M2/68‑10 (MQ=255)         ttGCCGCCTTCCTGCAACTCGAATTATTTCGAGTCTATGAATAATTTCTTAAGCATAGCAggagtgga                                                      <  2:858404/68‑1 (MQ=255)           gacgccATCCTGCAACTAGAATAATATAGAGTCTATGAATAATTTCTTAAGCAGAGCAggagtggagt                                                    >  2:419202/3‑68 (MQ=255)                      gCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCatagcaggagtggcagcatgcaataca                                         <  1:765096‑M2/68‑28 (MQ=255)                      gCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATgc                                         >  2:266806/1‑68 (MQ=255)                      gCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATgc                                         >  2:125795/1‑68 (MQ=255)                      gCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATgc                                         >  1:574336/1‑68 (MQ=255)                      gCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATgc                                         >  1:445545/1‑68 (MQ=255)                      gCAACTCGAATTAATTAGAGTCTATGAATAATTTCTTAGGCATAGCAGGAGTGGAGGAGGGATTATgc                                         >  2:901344/1‑68 (MQ=38)                         aCTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATgccgg                                      >  2:698634/1‑66 (MQ=255)                                    ttAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATGAAAGCCAGTCCTTGa                           >  1:12558/1‑68 (MQ=255)                                        agTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATGCCAGCCAGGCCTTGAtttt                       >  2:48088/1‑68 (MQ=255)                                         gTCTATGAATAATTTCTTAAGCatagcaggagtggcagcatgcaatacaaccccttaggaaaaaccga                      <  1:914960‑M2/68‑47 (MQ=255)                                              tGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATGCCAGCCAGGCCTTGATTTTGGCTTc                 >  2:497442/1‑68 (MQ=255)                                                aataatTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATGCCAGCCAGGCCTTGATTTTGGCTTCCa               >  2:306567/1‑68 (MQ=255)                                                   aatTTCTTAAGCatagcaggagtggcagcatgcaatacaaccccttaggaaaaaccgaccttcgcgtt            >  2:914957‑M2/1‑12 (MQ=255)                                                     tttCTTAAGCATAGCAGGAGTGGAGTAGGGATTATGCCAGCCAGGCCTTGATTTTGGCTTCCATACCa          <  1:48088/68‑1 (MQ=255)                                                           aaGCatagcaggagtggcagcatgcaatacaaccccttaggaaaaaccgaccttcgcgtttcccgact    >  2:765093‑M2/1‑4 (MQ=255)                                                             gCATAGCAGGAGTGGAGTAGGGATTATGCCAGCCAGGCCTTGATTTTGGCTTCCATACCAGCGGCATc  >  1:788666/1‑68 (MQ=255)                                                              |                                                                   TCACTCGTTTGCCGCCTTCCTGCAACTCGAATTATTTAGAGTCTATGAATAATTTCTTAAGCATAGCAGGAGTGGAGTAGGGATTATGCCAGCCAGGCCTTGATTTTGGCTTCCATACCAGCGGCATC  >  NC_000913/438232‑438359

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