New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	2808944 =	2 (0.090)	22 (1.020)	15/164	0.3	91.7%	coding (855/1557 nt)	gshA	glutamate‑‑cysteine ligase
	?	NZ_CP009273	3576788 =	NA (NA)					coding (391/417 nt)	BW25113_RS17825	YrhA family protein

TTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCA  >  NZ_CP009273/3576712‑3576861                                                                             |                                                                          tttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCaa                                                              >  1:364360/1‑90 (MQ=255)   tCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACt                                                            >  2:351300/1‑90 (MQ=255)      tAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTAc                                                         <  1:634418/90‑1 (MQ=255)            cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAttt                                                   <  2:364360/90‑1 (MQ=255)                    cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                           >  2:215249/1‑90 (MQ=255)                    cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                           >  2:91880/1‑90 (MQ=255)                                      ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTAAGGTGc                         >  1:81351/1‑90 (MQ=255)                                          aTTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                     >  2:565589/1‑90 (MQ=255)                                                  aTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTTCTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc             <  2:422890/90‑1 (MQ=255)                                                   tttCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCt            >  1:387329/1‑90 (MQ=255)                                                     tCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGt          >  1:6630/1‑90 (MQ=255)                                                             tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa  >  2:469909/1‑90 (MQ=255)                                                             tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa  >  2:565935/1‑90 (MQ=255)                                                                             |                                                                          TTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCA  >  NZ_CP009273/3576712‑3576861

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