New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	17 (0.470)	9 (0.260)	8/456	2.1	36.8%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	15 (0.440)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

TTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGC  >  NC_000913/1207562‑1208033                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    ttACCCGTGGTAACACCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTAGCGCGCTGCTGTTAAGCCGCAAGAAACCGTTTGGCGATATCACATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                                                                                                                                                                                                                >  2:243433/1‑250 (MQ=255)    cccGTGGTAACCTCCAGACAGGCACGCCTGCCACCTAACTGGATGATGGCTACTTTGACATGTTGCAGGAGGACCTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGCAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAATCCGCACGACTCCGTTTGGCGATATCAACTCGGATGGCCCTGTGCAACCGGCTCTCGAAACCCTTGGTTTGTGAGCAgggtacaaacac                                                                                                                                                                                                                            >  2:453045/1‑241 (MQ=255)      cGTGGTAACCCCCAGACCGGCACGCCCGCCACCGCTCTGGAATATGACTCCTTTCACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGCGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCACATCGGCTGGCACTGTGCACACGGCCCTCGAAAACCTTGCTTTGGGAGAAGGctcagcattacc                                                                                                                                                                                                                          >  2:57934‑M1/1‑239 (MQ=255)        tGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGATGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCTGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGAAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGTAGCAGAAATCAATg                                                                                                                                                                                                                        >  2:9876/1‑251 (MQ=255)            aaCCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACGTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgg                                                                                                                                                                                                                    <  1:184914‑M1/251‑19 (MQ=255)             aCCCCCAGACCGGCACGTCGTCCAACGAGCTTGAGGATGACTTCTTTGTCCTGTTGCCGGTGGAACTTGTCAGCTTGGTGGAGGCTCCCGTTGCCAGCCTGGGGAAGGGGCGGCACTACCATTTCCTTTCCGGACTTCGCGCGCTGCTTTTAAGTCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTCAGCAAACGGCTCTCGTAAACCTTGGTTTGGGAGaattcgctgcattccctgttggt                                                                                                                                                                                                                   <  2:358786/251‑22 (MQ=255)                cccAGACCGGCCCGCCGTCCGCCTATCGGGATGATTAGTACTTTGACAGGTTGCAGGAGGAATTTTGCAGCTGGGTTGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACTACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTATGCCGCAAGAATCCGTTTGTCGATATCAAATCGGATGGCACTGTGCAAACGGCGCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattaccttttggtgta                                                                                                                                                                                                                <  2:95450‑M1/251‑23 (MQ=255)                        ggCACGACTGCTACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCAGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTACATGCAGCAACGGCTACagt                                                                                                                                                                                                        >  2:70298/1‑249 (MQ=255)                                         tCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAGTCGGATGGCACTGTGCAACCGGCTCTCGCAAACCTTGGTTTGGGAGAAGGAGCACAACTCAATGCAGCAACGGCTACCTTAGGACGCACAGGTTTc                                                                                                                                                                                        >  2:392007/1‑250 (MQ=255)                                                    actactTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc                                                                                                                                                                            <  1:248958‑M1/251‑59 (MQ=255)                                                                 gtggaagggagaactttgctgcGTGGGGGAGTCATCCGGTGCCAGGCTGGATAAGGGTCGGCACGTCCAGTTACTGTCCTCACTTCGCGCGCTGCTGGTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACTGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattaccggttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgca                                                                                                                                                               <  2:343330‑M1/241‑72 (MQ=255)                                                                     caaaaatcagctcgcatcctgcacaggtcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaagccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGTTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTATATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                                                                                                                           >  2:131436‑M1/161‑251 (MQ=255)                                                                             agctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                                                                                                                            <  1:386901‑M1/74‑1 (MQ=255)                                                                             agctcgcaacctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttattgcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGTTTTCATAGCTATACCg                                                                                                                                                                            >  2:386901‑M1/153‑226 (MQ=255)                                                                                        gggtggATGCATACGGTTCCAGCCGGGAGAAGGTGCGGCTAGACCCGTTACGTACCGCACTTCTCTCGCTGCTGTGAAGCCGCAAGAAGCCGTTTGGCGATATCAAAGCGGATGGCACTGTGCAAACGGCTCTCGTAAACCTTGGTTTTGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACTTTAGTACTCATCTGTGTCATAGCTATACCGGTTATGATTTGTGGTATTGAGTAATCAGTAATCAtt                                                                                                                                        <  2:196349/250‑1 (MQ=255)                                                                                          cacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTc                                                                                                                                       <  1:23281‑M1/111‑1 (MQ=255)                                                                                              aGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTCGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgg                                                                                                                                  >  1:443481/1‑251 (MQ=255)                                                                                                            aGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTGGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccgctccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaa                                                                                                                    <  1:57934‑M1/251‑115 (MQ=255)                                                                                                                       gtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcagtaatgatgccagttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTGCATTAGGACGCATCGTTTGCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGTTGTAATGGCTCAAAAGCATCt                                                                                                          >  2:287976‑M1/111‑250 (MQ=255)                                                                                                                                    cAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGc                                                                                                                       <  1:110622/224‑1 (MQ=255)                                                                                                                                            tacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                                                                     <  1:311606/250‑1 (MQ=255)                                                                                                                                                 cACTTCGCGCGCTGCTGTTAAGACGCAAGAATCCGGTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTCGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATTCCGGTTATGATTGGTGGTATTGATCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCCGCCTCTGGGGGGGATGGAAATACAGtt                                                                               <  2:443481/251‑1 (MQ=255)                                                                                                                                                         caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTAGACCGGTTATGAGTGGTGGGATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTccc                                                                       >  1:289551‑M1/77‑251 (MQ=255)                                                                                                                                                                tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAACAGCATCTGCCTCTGGGGGGGATGGAAATACAGTGGTATTCCCGGTTGCg                                                                >  2:290366/1‑251 (MQ=255)                                                                                                                                                                 gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCg                                                                <  1:243433/250‑1 (MQ=255)                                                                                                                                                                   tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGtt                                                              <  1:70298/250‑1 (MQ=255)                                                                                                                                                                         tccccctggtaaaagcaggaatgatgccatttgcttgatatttcttccttttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGAATTCTCGGTTGCGTTTattaat                                                       <  2:452806‑M1/191‑1 (MQ=255)                                                                                                                                                                                ggtaaaatcaggaatgatgccatgtgcgtgttatttcttcattttcataaacgTTGGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTTGGGGTGGAATGCCGCAAAAGCATCGGCTTCTGGGGGGGATGGAAATACAGTTGTATTCCCTGTTTCGTTTAATAatgcttg                                                  >  2:23281‑M1/54‑246 (MQ=255)                                                                                                                                                                                     ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaa                                                                                                                 >  2:159531/1‑181 (MQ=255)                                                                                                                                                                                     ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaa                                                                                                                 <  1:159531/181‑1 (MQ=255)                                                                                                                                                                                      tGGCGATATCAAATCGGAGGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttagtcgtggctgggatgacgggcgggggatt                                           <  2:89241‑M1/250‑188 (MQ=255)                                                                                                                                                                                              tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggt                                  >  1:57118‑M1/1‑55 (MQ=255)                                                                                                                                                                                                                  cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattc               >  1:229503‑M1/1‑35 (MQ=255)                                                                                                                                                                                                                              gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc  <  1:392007/251‑1 (MQ=255)                                                                                                                                                                                                                               taaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGCCAATGTCAGCGCACCTATCAATGc  <  2:139639‑M1/244‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    TTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGC  >  NC_000913/1207562‑1208033

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