New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	269666 =	NA (NA)	16 (0.870)	10/166	0.5	94.1%	intergenic (+2/+146)	ykfC/BW25113_RS01325	reverse transcriptase domain‑containing protein/IS5‑like element IS5 family transposase
	?	NZ_CP009273	= 671829	1 (0.050)					coding (175/978 nt)	ybeQ	SEL1‑like repeat protein

CAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAA  >  NZ_CP009273/269580‑269751                                                                                       |                                                                                      cAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGaa                                                                                    <  1:472614/90‑1 (MQ=255)             cGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAg                                                                        >  2:497259/1‑90 (MQ=255)               agaTGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCg                                                                      >  1:366108/1‑90 (MQ=255)                     gAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGaaa                                                                >  2:32108/1‑90 (MQ=255)                                aCGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                                                     <  2:12453/90‑1 (MQ=255)                                aCGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                                                     <  2:366108/90‑1 (MQ=255)                                    tCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                                                 >  1:486661/1‑90 (MQ=255)                                      ccGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTc                                               <  2:434218/90‑1 (MQ=255)                                          aGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtc                                           <  1:352738/90‑1 (MQ=255)                                                aaGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttca                                     >  2:260542/1‑90 (MQ=255)                                                 aGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat                                    <  1:32108/90‑1 (MQ=255)                                                                             gCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc        >  2:432923/1‑90 (MQ=38)                                                                                cTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCg     >  1:423069/1‑90 (MQ=35)                                                                                  gTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGTCGATTTTTTCTCCCGTa   >  1:311734/1‑90 (MQ=25)                                                                                   ttATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa  >  2:403174/1‑90 (MQ=35)                                                                                       |                                                                                      CAGTCGTTATGGCGAGATGCGAGTGGTCTCAACGATCCCGCAGGAGAAAGCCAGAAACTTCGCCGCATCGCTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAA  >  NZ_CP009273/269580‑269751

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