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

GGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAT  >  NC_000913/1207650‑1207927                                                                                                                                             |                                                                                                                                          gcacagttcagcggattacgttgtagatctgtgttctggtcatttgttgagacgcgtaaatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTg                                                                                                                                         >  1:103455‑M1/141‑143 (MQ=255) gcacagtccagcgcattacgttggagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTg                                                                                                                                         >  1:80142‑M1/141‑143 (MQ=255)           tCCGGTGCCGGCTGGGAGAGGGGGCGTCATGACCATTTAGTTGCCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                                                                               <  2:37057‑M1/143‑1 (MQ=255)           gcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttctttcttcattttcataaacaTTGGtttgggcga                                                                                                                               <  2:149138‑M1/13‑4 (MQ=255)                                tgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                                                                          <  2:18652‑M1/34‑1 (MQ=255)                                 ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                                                                                         <  2:118170‑M1/143‑20 (MQ=255)                                      cACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                                                                                    <  1:59719‑M1/143‑25 (MQ=255)                                       tccaccaCTTACTTACCGCACTGCGGGTGCTGCTGTTAAGCCGCAAGATTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                                                                                   <  2:196705‑M1/140‑26 (MQ=255)                                               ttacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgctgtcccggttcagtg                                                                                           >  2:103047‑M1/1‑110 (MQ=255)                                                      ccGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttca                                                                                    <  1:96594‑M1/143‑41 (MQ=255)                                                                            ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACGTTAGGACGCACCGGTTTCATAGCTATACCGGTTa                                                              <  1:40919‑M1/78‑1 (MQ=255)                                                                              gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctg                                                            <  1:24997‑M1/143‑65 (MQ=255)                                                                                             ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgcc                                             <  2:188687‑M1/143‑80 (MQ=255)                                                                                             tagcaggaatgaagccatttgcatgttagttcttcactttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                             >  1:150763‑M1/49‑143 (MQ=255)                                                                                                     tggatgccagttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt                                     >  2:18275‑M1/41‑143 (MQ=255)                                                                                                       gatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACCCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTaa                                   >  1:18652‑M1/39‑143 (MQ=255)                                                                                                                            ctgcattttcgtaaaccTTCGCTTGGGAGAAGGAGCAAACTTCAATGCAGCAACGTCTACATTAGGACGCACCGGTTCCATAGCTATACCGGTCATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATg              >  2:133395‑M1/18‑143 (MQ=255)                                                                                                                                    tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcgaaggct      >  1:182017‑M1/1‑25 (MQ=255)                                                                                                                                        aaaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACGTGAGGACGCACCGGCTTCGTAGTCATCCCCGCCACGATCGCTGTCACTGAGCAATCACAAACCATTCAGCGGCCCTGGACTGCCGCACAAGCAt  >  2:153747/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                          GGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAT  >  NC_000913/1207650‑1207927

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