New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	6 (0.190)	7 (0.270)	7/128	1.9	58.9%	coding (290/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	1209619 =	5 (0.200)					pseudogene (37/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

CCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATG  >  NC_000913/1209546‑1209687                                                                          |                                                                     ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATg                                                               <  1:237424/81‑1 (MQ=255)  gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                               >  1:237457‑M2/73‑80 (MQ=255)           tgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                      <  1:237209‑M2/17‑1 (MQ=255)           gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                     <  2:237428/81‑1 (MQ=255)                   gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                             <  2:237210‑M2/26‑1 (MQ=255)                         cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                       >  1:237465‑M2/50‑81 (MQ=255)                          ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                      >  1:237459‑M2/49‑81 (MQ=255)                           cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATgg                                     >  2:237463/1‑81 (MQ=255)                                     gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                           >  1:237462‑M2/38‑81 (MQ=255)                                             aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                   <  2:237432/81‑1 (MQ=255)                                             aGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                   <  2:237433/81‑1 (MQ=255)                                                       gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt         >  1:237464‑M2/20‑81 (MQ=255)                                                          ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa      >  2:237491‑M2/17‑81 (MQ=255)                                                                 ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg  <  1:237434‑M2/69‑1 (MQ=255)                                                                          |                                                                     CCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATG  >  NC_000913/1209546‑1209687

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