New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	3 (0.200)	14 (1.070)	14/232	NT	85.8%	coding (290/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	1209619 =	2 (0.150)					pseudogene (37/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

TCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209478‑1209759                                                                                                                                              |                                                                                                                                             tCGGATAAGCCCTTGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                                                                                             <  1:134601/143‑1 (MQ=255)               gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                               <  1:170005‑M2/16‑1 (MQ=255)                              tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctaccaaaccaaCGTTTAGGAAAATGAAGAAAACACAAGCaaa                                                                                                                >  2:326‑M2/113‑143 (MQ=255)                                aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                              <  2:109996‑M2/33‑1 (MQ=255)                                           attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                   <  1:299667‑M2/44‑1 (MQ=255)                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctacattgaatattgctacttctaccaaacaaaCGTTTATGAAAAAAAAGAAATAACAAGCAAATGGAATCATTCATGAttt                                                                                              >  1:173268‑M2/95‑143 (MQ=255)                                                        accaatcataacctgtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      <  1:263839‑M2/57‑1 (MQ=255)                                                        accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATCCCTGCTCTTACCAggg                                                                                      >  1:203417‑M2/87‑143 (MQ=255)                                                             tcattacctgtatagctatgataccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAGTAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                                                                                 <  1:303774‑M2/62‑1 (MQ=255)                                                                 ccagcgTGGTGGCCGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                             <  2:254623/143‑1 (MQ=255)                                                                              atgaaagcggtgcgtccgaatgtagcagttgctggattgagttttgctccttcggccaaaccagCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTGTACCAGGGGGATTTAACATGCTTATTGGCt                                                                <  1:37189‑M2/79‑1 (MQ=255)                                                                               ttaaacctgtgcgtcctagtctaggcgttccttcattgagttttgctccttcctccaaaccaaCGTTTCTGATAATGAAGAAATAACAAGCTAATGGCATCATTCCTTCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                               <  2:300437‑M2/80‑1 (MQ=255)                                                                                           gacccGGAACCTGTACACCAACAGGTAATGCAGAGCCTTCTCGCCAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                   <  2:127619/139‑1 (MQ=255)                                                                                                aatgtagccgttgctgaagtgagttttgcgacttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaacaa                                              <  1:36072‑M2/97‑1 (MQ=255)                                                                                                   gtagctgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAGATAACAAGCAAATGGCATCATTCCTGCTTTTTCCAGGGGGATTTAACATGCTTATTGTCTATGTACGCGTATCAACAAATg                                           >  1:42146‑M2/44‑143 (MQ=255)                                                                                                               gcattgagttttgctccttctcacaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                               <  2:210877‑M2/112‑1 (MQ=255)                                                                                                                 attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGAATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                             <  2:288194‑M2/114‑1 (MQ=255)                                                                                                                               cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCCTCATTCCTGCTTTTAGCAGGGGGATTTAACATGCTTAGTGGCTATGTACGCGTAGCGACAAATGACCAGAAAACAGATCTACAAcggaagcg               >  2:271000‑M2/16‑140 (MQ=255)                                                                                                                                            aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  2:253934‑M2/141‑1 (MQ=255)                                                                                                                                              |                                                                                                                                             TCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209478‑1209759

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