New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	4 (0.140)	13 (0.520)	10/232	NT	85.2%	coding (290/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	1209619 =	1 (0.040)					pseudogene (37/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

AGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTG  >  NC_000913/1209485‑1209757                                                                                                                                       |                                                                                                                                           ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttctgctccttctcccaaaccaaAGTTTGTGa                                                                                                                                    >  2:334824‑M2/135‑143 (MQ=255)                                 ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                    >  2:110565‑M2/103‑143 (MQ=255)                                   gattcctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgatttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                                  <  2:38225‑M2/43‑1 (MQ=255)                                               ctaccaatcatcaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccacCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                      <  1:479527‑M2/55‑1 (MQ=255)                                                             cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAacagg                                                                        >  1:376839‑M2/75‑141 (MQ=255)                                                               gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctcctgctcccaaaacaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCCGCTTTGACCAGGGGGAGTGAACatgga                                                                      >  1:69007‑M2/73‑141 (MQ=255)                                                                               ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                      <  1:269856‑M2/87‑1 (MQ=255)                                                                                        taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcacca                                             >  1:174330‑M2/48‑143 (MQ=255)                                                                                        taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca                                             >  1:454786‑M2/48‑143 (MQ=255)                                                                                                        gcattgagttttgctcctgcccccaaaccgaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAaaatactc                             <  2:232791‑M2/112‑3 (MQ=255)                                                                                                          attgagttttgctccttctcccaaaccaaCGTGTATGAAAATGAAGAAATAACAACCAAATGGCATCATTCCTGGTTTTACCAGGGGAATTTAACATGCTTATTCGCTATGTACGCGTATCAACAAAGGACcagaacacggct                           >  2:236032‑M2/30‑139 (MQ=255)                                                                                                          aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                           <  2:458703/143‑1 (MQ=255)                                                                                                            tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                         >  1:230679‑M2/28‑143 (MQ=255)                                                                                                             gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTa                        >  2:349462‑M2/27‑143 (MQ=255)                                                                                                                              ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa       >  2:288779‑M2/10‑143 (MQ=255)                                                                                                                              ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa       >  1:301866‑M2/10‑143 (MQ=255)                                                                                                                                  accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAAAAAACATGATTATTGGCTATGTACGAATAATAACAAAGGACAAGAACATAGATTTACAACGAAATGAGCTGAACtgt   >  2:197017‑M2/6‑143 (MQ=255)                                                                                                                                   ccaaCGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgtg  <  1:207192‑M2/139‑1 (MQ=255)                                                                                                                                   ccaaCGTTCATGAGAATTAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgtg  <  1:502913‑M2/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                           AGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTG  >  NC_000913/1209485‑1209757

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.