New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	1203246 =	1 (0.040)	10 (0.490)	10/138	0.5	84.0%	coding (290/630 nt)	BW25113_RS05990	hypothetical protein
	?	NZ_CP009273	1205075 =	3 (0.150)					coding (40/540 nt)	BW25113_RS06005	phage tail protein

GCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGT  >  NZ_CP009273/1204989‑1205163                                                                                       |                                                                                         accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                       >  2:192143‑M2/87‑90 (MQ=255)      tcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                                                  >  1:628663‑M2/82‑90 (MQ=255)              ggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                          >  2:690183‑M2/74‑90 (MQ=255)                         gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                               >  2:585702‑M2/63‑90 (MQ=255)                           tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                                             >  1:207657/1‑90 (MQ=255)                               ggtgcgtcctaatgtagccgttgctgcattgagtttttctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                         >  2:527950‑M2/57‑90 (MQ=255)                                       aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                 >  1:56834/1‑90 (MQ=255)                                         aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                               >  2:121841‑M2/47‑90 (MQ=255)                                              agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                          >  2:326366‑M2/42‑90 (MQ=255)                                               gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                         >  1:660288‑M2/41‑90 (MQ=255)                                                             tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa                           <  1:86928/90‑1 (MQ=255)                                                                 tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                       <  1:648375‑M2/68‑1 (MQ=255)                                                                  ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                      >  1:374806‑M2/22‑90 (MQ=255)                                                                             cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa           <  1:527950‑M2/80‑1 (MQ=255)                                                                              ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt          <  2:679637‑M2/81‑1 (MQ=255)                                                                                      aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt  >  1:407122‑M2/2‑90 (MQ=255)                                                                                       |                                                                                         GCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGT  >  NZ_CP009273/1204989‑1205163

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