New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	2 (0.050)	13 (0.340)	13/464	1.1	90.1%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	1 (0.030)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAT  >  NC_000913/1207565‑1208050                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      cccGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTGTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatta                                                                                                                                                                                                                                             <  1:259567‑M1/251‑11 (MQ=255) cccGTGGAAACCCCCAGCCCGGCAAGCCTGCCACCGATCTGGCTGAGGAATACTTGGACATGCTGCAGGAGGCACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGACGGGGCGGCACGAACAGATACATACCGCACTTCGAGCGCTGCTATTAAGCCGCAAGAATCCGTTTGGCGATATACAATCGGATGGCACTGTCCAAACGGCTATCGAAAACCGTGGTTTGGGAGAAGGCGCAAAACTc                                                                                                                                                                                                                                             >  2:325970/1‑251 (MQ=255)                      tttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcgtggtatttcttcattttcataaacgTTGGTTTGGGAGAATGAGCAAAACTCAATGCAGCAACGGCTACAtt                                                                                                                                                                                                                         >  1:257408‑M1/205‑250 (MQ=255)                                             gatgaCTACTGTGACATGTTGCAGGAGGAACTGTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTGTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaac                                                                                                                                                                                                 <  1:85046‑M1/251‑55 (MQ=255)                                              tgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCGACGGCTACATTAGGACGCACCGGTTTCATAGCtat                                                                                                                                                                                                 >  1:382363‑M1/181‑250 (MQ=255)                                                                    ggaggaACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGcgctgcatttcctgttggtgtcccggttccgtggccgtcagccattccgccaacaggctggctgaaatgcaatggtg                                                                                                                                                                          <  1:26195‑M1/251‑78 (MQ=255)                                                                    aaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTg                                                                                                                                                                          <  1:56335‑M1/93‑1 (MQ=255)                                                                      aggaACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgcc                                                                                                                                                                        <  1:217479‑M1/251‑80 (MQ=255)                                                                        tcagctcgcatcctgcacagttcagcgcattacgttgtcgatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGc                                                                                                                                                                       <  1:286507‑M1/96‑1 (MQ=255)                                                                        gaACGTTGCAGCGTGGGGGAGGCATCCGGGGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgc                                                                                                                                                                      <  2:14060‑M1/251‑82 (MQ=255)                                                                                   cGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGAACTTCGCGCGATGCTGTTAAGCCGCACGACTCCGTTTGGCGATATCAAATCGGATGGCACTGCGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacatgttggggtcccggttacgtggccttcagccactccgccaacaggctggctgaaatgctatggttccgctttttctgctg                                                                                                                                                           >  2:296758‑M1/1‑159 (MQ=255)                                                                                    gtggtgGAGGCCTCCGGTGCCAGCCTGGAGAAGGGGCTGCACGACCAGTTACTTACCGCACGTCGCGCGCTGCTGTTAAGCTGCAAGAATCCGTTTGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattaccggttggtgtcccgttgccgtggccttcatccacttcgccaacaggctggctgaaatgcaatggtggctctttttctgctga                                                                                                                                                          <  2:355692‑M1/251‑94 (MQ=255)                                                                                       cacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCa                                                                                                                                                       >  2:56335‑M1/140‑251 (MQ=255)                                                                                         ggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaat                                                                                                                                                     >  1:352573‑M1/1‑153 (MQ=255)                                                                                           gttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttctgcattgtcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACTCGCCGTTTTCATAGCTATGCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgg                                                                                                                                                   >  2:432206‑M1/136‑251 (MQ=255)                                                                                                                agaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATTCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaactgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttat                                                                                                                              >  1:59710‑M1/1‑130 (MQ=255)                                                                                                                               atttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAATCATCTGCCTCTgggggg                                                                                                               >  2:186546‑M1/100‑251 (MQ=255)                                                                                                                                         tacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgattta                                                                                                      <  1:441005‑M1/250‑146 (MQ=255)                                                                                                                                             gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcgtatccgacaaataaattgcctgatttacgtgg                                                                                                 <  1:152695‑M1/251‑151 (MQ=255)                                                                                                                                                                                  ttGGAGATATCAAATAGGAAGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcacggtgccgtgtccttcagccactcagccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggaaaaggcttatccgacaattacatttcctgatttcagtggtgagtttattcgtggctgggatgacgggggggggatt                                                            >  2:12864‑M1/1‑64 (MQ=255)                                                                                                                                                                                                             atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                                 <  1:444149‑M1/230‑1 (MQ=255)                                                                                                                                                                                                               cAAACGGCTCTCGAAACCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccacgccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctatttgaagtattca                               <  2:425060‑M1/251‑217 (MQ=255)                                                                                                                                                                                                                                            agaAGGAGCAAAACTCCATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTAGGACAATGTCAGCGCACCTATCAATGCATTGGCAACGGGGGGat  <  1:136830/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      CCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAT  >  NC_000913/1207565‑1208050

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