New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 3512449	2 (0.220)	12 (1.320)	10/166	0.1	85.7%	intergenic (‑26/+375)	aroK/hofQ	shikimate kinase AroK/DNA uptake porin HofQ
	?	NZ_CP009273	3576788 =	NA (NA)					coding (391/417 nt)	BW25113_RS17825	YrhA family protein

TAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGT  >  NZ_CP009273/3576703‑3576844                                                                                      |                                                         taatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGa                                                      <  1:198577/90‑1 (MQ=255)                     cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAttt                                  >  1:241931/1‑90 (MQ=255)                       gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAAc                                            >  1:137272/1‑78 (MQ=255)                       gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAAc                                            <  2:137272/78‑1 (MQ=255)                                 aaCGGTTAGTGTTCGGGGTTTAAAGCATTTCAATATTTACTTATGCCATTAAAGGTGATGCTGCCAACTTCCTGATTTAGTGTATGATgg                      >  2:213428/1‑90 (MQ=37)                                  aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgt                     <  2:15359/90‑1 (MQ=255)                                   cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAgg                                                         <  1:81768/53‑1 (MQ=255)                                   cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAgg                                                         >  2:81768/1‑53 (MQ=255)                                             tCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt          >  2:184728/1‑90 (MQ=255)                                               ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGc        >  1:255121/1‑90 (MQ=255)                                                     tataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAAGTGCTCCAGt  >  1:226291/1‑90 (MQ=255)                                                                                      |                                                         TAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGT  >  NZ_CP009273/3576703‑3576844

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