New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 1972751	0 (0.000)	16 (1.260)	9/168	0.3	100%	intergenic (‑297/‑483)	flhD/uspC	flagellar transcriptional regulator FlhD/universal stress protein UspC
	?	NZ_CP009273	= 2283592	NA (NA)					intergenic (‑68/‑387)	BW25113_RS11435/narP	IS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP

CTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972672‑1972840                                                                                |                                                                                          cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 >  2:145063‑M1/1‑80 (MQ=255)           aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       >  1:233571‑M1/1‑70 (MQ=255)                  ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg                                                                >  2:142357‑M1/1‑63 (MQ=255)                    aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                              >  2:275625‑M1/1‑61 (MQ=255)                     ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                             >  1:74948‑M1/1‑60 (MQ=255)                                       ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  1:128246‑M1/90‑49 (MQ=255)                                       ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  2:74948‑M1/90‑49 (MQ=255)                                            gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                      <  2:233571‑M1/90‑54 (MQ=255)                                                 aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                 >  1:110738‑M1/1‑32 (MQ=255)                                                                   acTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtca               <  2:110738‑M1/90‑77 (MQ=255)                                                                               aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccagagaacagggttcatcatgagtcatcaacttacctt   >  2:247724‑M1/1‑2 (MQ=255)                                                                                aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  1:28807‑M1/1‑1 (MQ=255)                                                                                aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  <  2:347926‑M1/90‑90 (MQ=255)                                                                                aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                            <  1:21367‑M1/64‑64 (MQ=255)                                                                                aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                            >  2:21367‑M1/1‑1 (MQ=255)                                                                                aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggacccatagaacaggtttcatcatgagtcatcaacttaccttc  >  1:94947‑M1/1‑1 (MQ=255)                                                                                |                                                                                          CTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972672‑1972840

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