New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	9 (0.310)	18 (0.670)	17/452	0.5	72.9%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	5 (0.190)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

GACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGT  >  NC_000913/1207581‑1208036                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        gACCGGCCCGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGGATGGCGATATCAAATCGGATGGCACTGGGCAAACGGCTTTCGAAAACCTTGGTGTGGGAGAAGGctctgcattacctgctggtgtcccgg                                                                                                                                                                                                               >  2:218042‑M1/1‑225 (MQ=255) gACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgg                                                                                                                                                                                                               <  1:165757‑M1/251‑27 (MQ=255)     ggCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGTTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                                                                                                                                                                                           <  1:48056/251‑1 (MQ=255)                 aCCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGAGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                                                                                                                                                                                               >  2:274158/1‑251 (MQ=255)                      tCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc                                                                                                                                                                                          >  1:33248‑M1/1‑204 (MQ=255)                         ggATGATGACTACTGTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCTCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcc                                                                                                                                                                                         <  2:289438‑M1/249‑49 (MQ=255)                                gaCTACTTTGCCATGTTGCAGGCGGAACTTTGCAGCGTGGTGGCGGCATCCGGTGCCAGCCTGGAGCAGGGGCGGCACGCCCAGTTACTTACCGCACTTCGCGCGCTGCTGTTACGCCGCAAGAATCCGTTTGGCGATCTCAAATCGGATGGCACTGTGCAACCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg                                                                                                                                                                                >  2:204576‑M1/1‑194 (MQ=255)                                              gTTGCAGGAGGAACTTGGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGTGCGGCGCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTTCAATCGGCTCTCGAAAACCTTGGTTTGGTAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgc                                                                                                                                                                   <  2:267377‑M1/250‑71 (MQ=255)                                                tGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTGCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgattgtgtcccggttccgtggccgtcagccactccgccaacaggctggctgaaatgcaa                                                                                                                                                                 <  1:204576‑M1/250‑73 (MQ=255)                                                  cAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                                                                                                                              >  1:156866/1‑251 (MQ=255)                                                                 aGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGAGGAAAGCCCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCACGACTCCGCTTGGCGATATCAAATCGGATGCCCCTGTGCAAACGCCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAt                                                                                                                                               >  1:26935/1‑251 (MQ=255)                                                                          gAGGCATCCGGTGCCTGCCTGGAGAAGGTGCGGCACGACCAGTTATTTACCGCCCTTCGCTCGCTGATGTTAAACCGCACGTATCCGGTTGGCTAGATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGTTTTGGGAGAAGGccctgcattacctgttggtgtcccggttccgtggccttcatccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaata                                                                                                                                      <  2:238835‑M1/251‑100 (MQ=255)                                                                                 cgcattacggtgtagatctgtggtctggtcatttgttgatacgcgtacttagacaataagcctgttaaatccccctggtaaaagcaggtatgatgcgagttgcgtgttatttcttcagtttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTgg                                                                                                                               >  2:112677‑M1/130‑251 (MQ=255)                                                                                     ttacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaatcagcatgttaaatccccctggtcaaagcaggaatgatgccctttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCaga                                                                                                                                       <  2:13000‑M1/114‑2 (MQ=255)                                                                                     ttacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt                                                                                                                                       >  1:13000‑M1/126‑239 (MQ=255)                                                                                          gCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGAGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactagcaaag                                                                                                                      <  1:91199‑M1/251‑116 (MQ=255)                                                                                                                     tacttacCGCACTGCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGAGGGCACTGTGCAAACGGCTCTCGAAACCCTAGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctga                                                                                            <  2:240217‑M1/250‑142 (MQ=255)                                                                                                                             gtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGt                                                                                   >  1:197400‑M1/86‑251 (MQ=255)                                                                                                                                      gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATTCAGCAACGGCTACATTAGGACGCTCCGGTTTTATAGCTATACCGGTTATGATTTGTGGTATTGAGCAATCAGTTATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAGCTGCCGCTGGGGGGGATGGAAATACAGTTGTATTccc                                                                          <  1:107447/251‑1 (MQ=38)                                                                                                                                       cgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttatt                                                                         <  1:18305‑M1/251‑161 (MQ=255)                                                                                                                                                    ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatg                                                            <  1:276950‑M1/251‑174 (MQ=255)                                                                                                                                                     cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagactacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatg                                                            >  1:150030‑M1/1‑77 (MQ=255)                                                                                                                                                         aGAATCTGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacggg                                                       >  1:92191‑M1/1‑73 (MQ=255)                                                                                                                                                            aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCtgt                                                    <  1:270211/251‑1 (MQ=255)                                                                                                                                                              ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGTTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGTATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTgtgt                                                  <  1:274158/251‑1 (MQ=255)                                                                                                                                                                aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGGTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGtt                                                 >  2:198801‑M1/51‑250 (MQ=255)                                                                                                                                                                    gcaggaatgatgccatttgcttgttatttcttccttttcctaaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGATGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTACTAATGCCTGTGTTGCCgt                                            <  2:181876‑M1/205‑1 (MQ=255)                                                                                                                                                                                 cGGATGGGACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAGTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTTGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATTGAAATACAGTTTTATTCCCGGTTGCGTTTAATAATTCCTGTGTTGCCGCTGTTGCAAATTAt                               <  2:48897/251‑1 (MQ=255)                                                                                                                                                                                             atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTACTAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                   <  1:265987‑M1/230‑1 (MQ=255)                                                                                                                                                                                              tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcg                  <  1:1314‑M1/231‑1 (MQ=255)                                                                                                                                                                                                 aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggcactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattcagg               >  1:213514‑M1/1‑33 (MQ=255)                                                                                                                                                                                                   cGGCTCTCGAAATCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacggggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattcagggg             <  2:213514‑M1/251‑221 (MQ=255)                                                                                                                                                                                                     attttcataaatgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAGTATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGGTGTATTCCCGTTTGCGTTTAATAATGTCAGTGTTGCCGGTGTTGCAAattattacaatgttatcccgcctat           >  2:46197‑M1/14‑239 (MQ=255)                                                                                                                                                                                                          cataaacgTTGGTTTGGGAGAAGACGCAAAAATCAACGCACCAACGGCAACATTAGGCCGCACCAGTTTCATAGCTCTACCGGCTATGATTGGCGGCATTGAGCAATCAGTAATCATTCACTGAGGGTGGAATGCCGCAACAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTCTTCCCGGTTGCGTTTCATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAACGCACCTatcactg      <  2:16353‑M1/243‑4 (MQ=255)                                                                                                                                                                                                          cGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattcaggggtatgca       <  1:244914‑M1/250‑227 (MQ=255)                                                                                                                                                                                                          cGAAAACCTTGGTTTGGGAGAAGGctctgcattaccggttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattcagtggtatgca       <  2:150030‑M1/250‑227 (MQ=255)                                                                                                                                                                                                              gaacttgggctGGGAGAAGGAGCAAAACTCAATGCAGTAACGGATACATTAGGACGCAACGGTTTAATAGCTATAACGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGAATCTGCCTATGGGGGGGATGGAAAAACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTCTGACAATGTCAGCTCACCTATCAATGCAgg  >  2:18820/11‑250 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        GACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGT  >  NC_000913/1207581‑1208036

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