New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	19797 =	NA (NA)	15 (0.700)	11/162	0.5	83.5%	intergenic (‑2/+14)	BW25113_RS25960/BW25113_RS00100	fimbrial family protein/IS1‑like element IS1A family transposase
	?	NZ_CP009273	= 2808953	3 (0.140)					coding (846/1557 nt)	gshA	glutamate‑‑cysteine ligase

TTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGA  >  NZ_CP009273/19713‑19882                                                                                     |                                                                                      ttGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATg                                                                                  >  1:595209/1‑90 (MQ=255)  tGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGa                                                                                 >  2:99872/1‑90 (MQ=255)  tGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGa                                                                                 >  2:365864/1‑90 (MQ=255)       ccATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCa                                                                            <  1:266217/90‑1 (MQ=255)               aCGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTg                                                                    >  1:339686/1‑90 (MQ=255)               aCGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTg                                                                    >  2:475757/1‑90 (MQ=255)                   gggCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAg                                                                <  2:563858/90‑1 (MQ=255)                   gggCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAg                                                                <  2:512295/90‑1 (MQ=255)                          taaataGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTa                                                         <  2:112414/90‑1 (MQ=255)                               agagGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTc                                                    <  2:86670/90‑1 (MQ=255)                                  ggAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGa                                                 >  1:302666/1‑90 (MQ=255)                                  ggAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGa                                                 >  1:558216/1‑90 (MQ=255)                                    aaTCTGATTACTTCCTTCATTGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATa                                               >  2:475066/1‑90 (MQ=255)                                            tACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGa                                       <  1:618518/90‑1 (MQ=255)                                                tccttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGAc                                   <  1:365864/90‑1 (MQ=255)                                                  cttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACtt                                 <  1:545899/90‑1 (MQ=255)                                                   ttcATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACttt                                <  2:595209/90‑1 (MQ=255)                                                                   aaaGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCAcc                >  2:71003/1‑90 (MQ=255)                                                                      gagTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTc                                                    <  2:636142/51‑1 (MQ=255)                                                                      gagTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTc                                                    >  1:636142/1‑51 (MQ=255)                                                                        gtagtaATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGAttt           >  2:249338/1‑90 (MQ=255)                                                                             aaTTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTgaga      <  2:395629/90‑1 (MQ=39)                                                                                tGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACg   >  2:600748/1‑90 (MQ=35)                                                                                tGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACg   >  1:128845/1‑90 (MQ=35)                                                                                tGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACg   >  1:552833/1‑90 (MQ=35)                                                                                 gCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGa  <  1:244874/90‑1 (MQ=33)                                                                                     |                                                                                      TTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGA  >  NZ_CP009273/19713‑19882

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