New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	1 (0.100)	8 (0.890)	8/232	NT	94.8%	coding (305/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	= 1209602	0 (0.000)					pseudogene (54/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

GGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGG  >  NC_000913/1207681‑1207946                                                                                                                             |                                                                                                                                              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)                                                                                                                       gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtg                                                                  >  2:90896‑M1/1‑7 (MQ=255)                                                                                                                       gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtg                                                                  <  1:90896‑M1/84‑78 (MQ=255)                                                                                                                        agaAGGgtctgcattacttgttggtgtcccggttccgtggtcttcagctactccgcgaacaggctggctgaaatgcaatggtgccgctttttctgctgaagcatacccgaaactggcaaaggcttatccgacaaataaattgc      <  2:51343‑M1/143‑138 (MQ=255)                                                                                                                            GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCGGCAGAAGCATCTGCCTCTGGGGgggacgg  <  1:155747‑M1/143‑4 (MQ=255)                                                                                                                             |                                                                                                                                              GGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGG  >  NC_000913/1207681‑1207946

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.