New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 1391494	NA (NA)	7 (0.510)	6/164	0.7	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

GAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGG  >  NZ_CP009273/1391416‑1391572                                                                               |                                                                               gagaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATc                                                                     >  1:343341/1‑90 (MQ=37)                           aaatgaaatgACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTaa                                           <  2:195773/90‑1 (MQ=255)                                      cTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTa                                            >  1:169932/1‑78 (MQ=255)                                      cTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTa                                            <  2:169932/78‑1 (MQ=255)                                         aGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTc                             <  2:148512/90‑1 (MQ=255)                                           tCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCag                           <  1:101651/90‑1 (MQ=255)                                                  agaagaATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATcc                                                 <  1:281621/61‑1 (MQ=255)                                                  agaagaATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATcc                                                 >  2:281621/1‑61 (MQ=255)                                                   gaagaaTTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACa                   >  1:121410/1‑90 (MQ=255)                                                            ccccGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGaaa          <  1:170957/90‑1 (MQ=255)                                                                  ttATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCAt    <  2:343341/90‑1 (MQ=255)                                                                    aTTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATgg  <  1:275026/90‑1 (MQ=255)                                                                               |                                                                               GAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGG  >  NZ_CP009273/1391416‑1391572

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