New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	376748 =	NA (NA)	3 (0.120)	3/152	2.0	NA	intergenic (‑5/‑59)	BW25113_RS01855/BW25113_RS01870	transferase/IS3‑like element IS2 family transposase
	?	NZ_CP009273	376768 =	NA (NA)					intergenic (‑25/‑39)	BW25113_RS01855/BW25113_RS01870	transferase/IS3‑like element IS2 family transposase

ATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAG  >  NZ_CP009273/376685‑376801                                                                                    |                                  attGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                             >  1:16446/1‑90 (MQ=255)    gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAc                          >  2:113805/1‑90 (MQ=255)      gCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                        <  1:464441/90‑1 (MQ=255)      gCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                        <  1:47802/90‑1 (MQ=255)      gCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                        <  2:16446/90‑1 (MQ=255)       ccAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAg                       >  1:314101/1‑90 (MQ=255)         aGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                     >  2:546722/1‑90 (MQ=255)         aGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                     >  2:699613/1‑90 (MQ=255)              gAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                >  1:204147/1‑90 (MQ=255)                aTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc              <  2:468087/90‑1 (MQ=255)                aTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc              <  2:86369/90‑1 (MQ=255)                     cGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc         <  1:413633/90‑1 (MQ=255)                     cGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc         <  1:699613/90‑1 (MQ=255)                            ttACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTCAGACTAGTCATTGgag  <  1:398666/90‑1 (MQ=17)                                                                                    |                                  ATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAG  >  NZ_CP009273/376685‑376801

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