New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	10 (0.260)	12 (0.340)	10/410	1.9	49.7%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	15 (0.430)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACT  >  NC_000913/1207545‑1208061                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                 cgcgGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                                                                                                                                                                                                                           >  1:168534‑M1/1‑261 (MQ=255) cgcgGGTAAGAACGGGTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGATTACTTTGACATGTTGCAGGAGGAAGTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGATCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCCATGCAGCAACGGc                                                                                                                                                                                                                                          <  1:562527/285‑1 (MQ=255)                 ctgccgataatgtcctgagcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                                                                                                                                                                                                                >  1:266874‑M1/230‑279 (MQ=255)                                                                                   ttGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccgga                                                                                                                                                         >  1:157865‑M1/1‑179 (MQ=255)                                                                                             aaCTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATAGCACTGTGCAAACGGCTCTCGAAAACCTTGGTCTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggacctggcaaagg                                                                                                                                              >  1:513811‑M1/1‑169 (MQ=255)                                                                                              gcatggcgcgggggtggGGGCTCCGGGTGCCGGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGAACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGaa                                                                                                                                                      <  2:446060/264‑1 (MQ=37)                                                                                                         tggtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAGCAGTAATCATTCAGTGGGGGAGGAAGGCAGCACAAGCATCTGcca                                                                                                                                  >  1:454727/1‑284 (MQ=255)                                                                                                                          tacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGt                                                                                                                                                          <  1:374395‑M1/120‑1 (MQ=255)                                                                                                                                 tatatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgacgccatttgattgttatttcttcattttcataaacgTTGGTTTGGGAGAAAGAGCAAAACTCAATGCAGCAACGTCTACATTATGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt                                                                                                                                                                >  1:546853‑M1/118‑231 (MQ=255)                                                                                                                                       tgtgttctggtcatttgttgatacgcgtacatagccactaagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCCGTATTCATTAAGTGGGGGt                                                                                                                                                          <  2:546853‑M1/120‑1 (MQ=255)                                                                                                                                       tggggCAGCACGACAGGTTACTTACCGCCCGTGCACCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc                                                                                                                                    <  2:454727/252‑1 (MQ=38)                                                                                                                                                                  tacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCATGTGTTg                                                                         <  1:469397‑M1/201‑1 (MQ=255)                                                                                                                                                                             taagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGG‑ATGGAAATACAg                                                                                                             <  1:190355‑M1/164‑1 (MQ=255)                                                                                                                                                                             taagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGG‑ATGGAAATACAg                                                                                                             >  2:190355‑M1/74‑237 (MQ=255)                                                                                                                                                                             ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGACTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcgactccgccaacagtctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgccggatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacagg                                                               >  1:577131‑M1/1‑89 (MQ=255)                                                                                                                                                                                   tgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                                                                                                                                             <  2:518069‑M1/85‑1 (MQ=255)                                                                                                                                                                                   tgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                                                                                                                                             >  1:518069‑M1/68‑152 (MQ=255)                                                                                                                                                                                      aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                           >  1:510058/1‑103 (MQ=255)                                                                                                                                                                                      aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                           <  2:510058/103‑1 (MQ=255)                                                                                                                                                                                          gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggatt                                                                       >  1:472089‑M1/1‑76 (MQ=255)                                                                                                                                                                                                            gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg                                                                                                                                                                            <  1:324078‑M1/102‑1 (MQ=255)                                                                                                                                                                                                            gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg                                                                                                                                                                            >  2:324078‑M1/43‑144 (MQ=255)                                                                                                                                                                                                                      atttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTGGCTTTTAATAATGCCTGTGTTGCCGTTGTTGCAACTTATGAAAATGTCagcccacctatcaagccatgggcaac                     <  1:730‑M1/253‑13 (MQ=255)                                                                                                                                                                                                                       gtggttttttttttcttcattttcataaccgTGAGTTTGGGAGAAGGAGCAATACTAAATGAAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAAAGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACg                    >  2:469397‑M1/32‑285 (MQ=255)                                                                                                                                                                                                                           gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGTGGGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCCTTCAGTAATCATTCAGTGGGGGTGGAAGGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATGCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                                         <  1:25825/260‑1 (MQ=255)                                                                                                                                                                                                                           ccctgcagaaaactgctcgcgAAACCCGGGGTTGGGGAGAAAGgtcatcatcaccgttggtggtcccggttccgtggccttctgccactccgccaaaaggctggctgaaaggcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaaaaaattgcctgatttacgtggtgagttgattcgtggctgggatgacgagcgtgggattgatacaggtcgctctattttaagtattcaggggtatgcaacggaggatcatgctc                <  2:513811‑M1/276‑243 (MQ=255)                                                                                                                                                                                                                               ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                           <  2:219585‑M1/39‑1 (MQ=255)                                                                                                                                                                                                                               ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                           >  1:219585‑M1/24‑62 (MQ=255)                                                                                                                                                                                                                                 tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattc                                           >  1:108210‑M1/1‑37 (MQ=255)                                                                                                                                                                                                                                 tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                                                                                                                                                      <  2:216705/49‑1 (MQ=255)                                                                                                                                                                                                                                 tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                                                                                                                                                      >  1:216705/1‑49 (MQ=255)                                                                                                                                                                                                                                      aCGCCTCTCGAAAACCTTGGCTTAGGAGCAGGAGCAACACTCCATGCTGCAAGGGCTACATTAGGCCGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTCTTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACa                                                    >  1:457561/1‑238 (MQ=255)                                                                                                                                                                                                                                         gCTCTCGAAGACCTTGGTTTGGGAGAAGGctctgcagtacctgttggtgtcccggttccgtggccttccgccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattcaggggtatgcaacggaggatcatgctcatggattaccgtca  <  1:573555‑M1/285‑257 (MQ=255)                                                                                                                                                                                                                                           tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc                                                                                                                                                                                                        <  1:571645‑M1/85‑59 (MQ=255)                                                                                                                                                                                                                                           tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc                                                                                                                                                                                                        >  2:571645‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                           tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccg                                                                                                                                                                                                                                   <  2:463125‑M1/58‑32 (MQ=255)                                                                                                                                                                                                                                           tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccg                                                                                                                                                                                                                                   >  1:463125‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                                                   AGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATc                                   <  1:541496‑M1/226‑1 (MQ=255)                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                 CGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACT  >  NC_000913/1207545‑1208061

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