New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	0 (0.000)	9 (0.370)	8/232	NT	90.0%	coding (290/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	1209619 =	2 (0.080)					pseudogene (37/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

TTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209490‑1209760                                                                                                                                  |                                                                                                                                              cggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                                                                                                  >  1:325327‑M2/130‑143 (MQ=255)                  ggagtgcttactgagggctcaagaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaaggtagccggtgctgcattgagttttgctcattctcccaaaccaaCGTTTATGAAAATGAAGAAATAACCAGCcaa                                                                                                                 <  2:222959‑M2/31‑1 (MQ=255)                    ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCAAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCACCAGGTAATGCAGAGCCTTCGACCACACCAACGTTCATGACAACGAAGAAATAGCAAGCaaacg                                                                                                               >  2:308829/1‑141 (MQ=255)                     atgattactgattgctcaataccaccaatcataaccgggatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                                                              <  1:243103‑M2/34‑1 (MQ=255)                                         accaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                          <  1:404778‑M2/54‑1 (MQ=255)                                           caccaatcataatcggtatagctatgaaatctgttcgtcctcatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCATATGGCATCATTCCTGCTTTTACCAgg                                                                                        <  2:148220‑M2/56‑1 (MQ=255)                                                                ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                   <  1:287449‑M2/77‑1 (MQ=255)                                                                                        tagccgttgctgcatttagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTCTTGGCTATGTACGCGTATCAACAAATGa                                           >  2:151392‑M2/43‑143 (MQ=255)                                                                                        tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTGAACATGCTGATTGGCTATGAACGCGTATGAACAAATGa                                           >  1:12901‑M2/43‑143 (MQ=255)                                                                                                gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                   >  1:455463‑M2/35‑143 (MQ=255)                                                                                                  aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATGTAACATGCTTATTGGCTATGTACGCGTAGCAACAAATGACCAGACcaca                                 >  2:307345/1‑143 (MQ=255)                                                                                                                             accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGGGATTTAACATGCTTATTGGCTATGTAAGCGTATCAACAAATGACCAAAACAAAGACCTACAACGTAATGAgctgaaatat      >  2:238525‑M2/6‑139 (MQ=255)                                                                                                                                 aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg  <  2:108846‑M2/142‑1 (MQ=255)                                                                                                                                  |                                                                                                                                              TTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209490‑1209760

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