New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 1391494	NA (NA)	9 (0.680)	5/166	0.8	100%	intergenic (+145/‑128)	BW25113_RS06980/ynaJ	IS5‑like element IS5 family transposase/DUF2534 family protein
	?	NZ_CP009273	1972747 =	0 (0.000)					intergenic (‑293/‑487)	flhD/uspC	flagellar transcriptional regulator FlhD/universal stress protein UspC

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCT  >  NZ_CP009273/1972660‑1972814                                                                                        |                                                                    atttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcctt                                                                   <  1:113484‑M2/3‑1 (MQ=255) atttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcctt                                                                   <  1:203245‑M2/3‑1 (MQ=255)       tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                                                             <  1:91268‑M1/90‑5 (MQ=255)       gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                             <  1:89494‑M2/9‑1 (MQ=255)       gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                             <  2:96647‑M2/9‑1 (MQ=255)        cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                            >  1:252217‑M1/1‑85 (MQ=255)         gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                           <  1:170145‑M2/11‑1 (MQ=255)                       aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                             <  1:68885‑M1/90‑21 (MQ=255)                        tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct                                            <  2:252217‑M1/90‑22 (MQ=255)                            aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                        >  1:374418‑M1/1‑65 (MQ=255)                               gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatga                                     >  1:289345‑M1/1‑62 (MQ=255)                                 ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                   >  1:200637‑M1/1‑60 (MQ=255)                                 ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                   >  1:89772‑M1/1‑60 (MQ=255)                                 agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                   <  2:311974‑M2/35‑1 (MQ=255)                                              tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcat                      >  1:125191‑M1/1‑47 (MQ=255)                                               tgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                     <  2:285017‑M2/49‑1 (MQ=255)                                               tgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                     <  2:285016‑M2/49‑1 (MQ=255)                                                             aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac       <  1:241545‑M1/90‑59 (MQ=255)                                                                 tttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc   <  1:118385‑M2/67‑1 (MQ=255)                                                                  gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt  <  2:125191‑M1/90‑64 (MQ=255)                                                                  gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt  <  2:103289‑M1/90‑64 (MQ=255)                                                                  gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgaagacatagaacagggt  >  2:176188‑M1/1‑27 (MQ=255)                                                                                        |                                                                    CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCT  >  NZ_CP009273/1972660‑1972814

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