New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	22 (0.310)	21 (0.320)	21/458	1.5	55.6%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	13 (0.200)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

GTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTG  >  NC_000913/1209375‑1209866                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        gTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                                                                                                                                   >  1:607647/1‑251 (MQ=255)    gCATACCCCTCAATACTTCAAATAGAGCAACCTGTATCCATCCCCCGCCCGTCATCCCATCCACGCATAAACTCACCACGTCCATCCGGCCATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGACAACGCGGCACCATTGCATTTCCGCCAGCCTGTTGGCGGAGTGGCTGAACGCCACGGCACCGTGACCCCCACATGTAATGCAGCGCCTTCTCCCAACCCAAGGttttctcg                                                                                                                                                                                                                                                >  2:141221/1‑246 (MQ=255)             gcgttgacattgccataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatccccaccagaggcagaggctgttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                                                                                                                                                                                                                        <  2:13267‑M2/18‑1 (MQ=255)                   acattgtcataatttgcaacaacggcaacacaggcattattaaacgcaacagggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                                                                                                                                  >  1:321975‑M2/227‑250 (MQ=255)                        gtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                                                                                                                                                                                                             <  1:89821‑M2/29‑1 (MQ=255)                         aTAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                                                                                                                                           <  1:677304/251‑1 (MQ=255)                             agCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                                                                                                                                                                       <  1:573393/251‑1 (MQ=255)                             aatttgcaacaacggcatcacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgcatgattactgattgctcaataccaccaatcataacctttatagctatgaaaacggtgcgtcctaatgtagccgttgctgcattgctttttgctccttctcaccaaccaaCGCTTATGAAAATGACGAAATAACAAGCAAATGGc                                                                                                                                                                                                                       >  2:627696‑M2/217‑251 (MQ=255)                                           aTCCCCCGCCCGTCATCCCAGCCTCGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                                                                                                                                         <  2:6990/251‑1 (MQ=255)                                             aacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                                                                                                                                        >  1:334361‑M2/201‑250 (MQ=255)                                                             aacgcaaccgggaatccaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcggcctaatgtagccgttgctgcattgaggtttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                                                                                                                       <  2:334361‑M2/67‑1 (MQ=255)                                                                  aaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccacttaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTGAACATTCtt                                                                                                                                                                                  <  2:35936‑M2/72‑1 (MQ=255)                                                                   accgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                                                                                                                                  >  1:484126‑M2/179‑250 (MQ=255)                                                                           tacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                                                                                                                                        <  1:23882‑M2/50‑1 (MQ=255)                                                                           tacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                                                                                                                                        >  2:23882‑M2/171‑220 (MQ=255)                                                                                    atttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaatacaaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTCCGCGt                                                                                                                                                                 >  1:133658‑M2/162‑250 (MQ=255)                                                                                                    gaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                                                                                                                                          <  1:233006‑M2/32‑1 (MQ=255)                                                                                                    gaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                                                                                                                                          >  2:233006‑M2/146‑177 (MQ=255)                                                                                                               ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                                                                                                                      >  1:207931‑M2/135‑250 (MQ=255)                                                                                                                tttgcggcattccacccccactgaatgattactgattgctcaatacccccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                                                                                                                      <  1:492841‑M2/116‑1 (MQ=255)                                                                                                                                      gaatgattactgattgctcaaacccaccaatcataaccggtatagctattaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt                                                                                                                <  2:484126‑M2/138‑1 (MQ=255)                                                                                                                                       ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGTAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCTGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTCATGCGCTGAACTGTGc                                                                                                              >  1:644040/1‑250 (MQ=255)                                                                                                                                                        caataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccatCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAATACAGATCTACAACGTTATGCGCTGAACTGTGCAGGATGCGAGCTGAtt                                                                                              <  1:50616‑M2/156‑1 (MQ=255)                                                                                                                                                               gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAATCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                                                                                                                      <  2:500140/202‑1 (MQ=255)                                                                                                                                                               gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                                                                                                                      >  1:500140/1‑202 (MQ=255)                                                                                                                                                               accaatcataaccggtatatctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttcttccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATTGCATCATTCCTGCTTTTACCAGTGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAGGACCAGAACACAGATCTACAAAGTAATGCGCTGAACTGTGCAGGATTCGAGCTTATTTTTGAAGAc                                                                                     <  2:132534‑M2/165‑1 (MQ=255)                                                                                                                                                                         accggcatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgcgttttgctccttctcccaaaccaaCGTTTATGACAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACAAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCg                                                                           >  2:105639‑M2/77‑251 (MQ=255)                                                                                                                                                                                     ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCg                                                               <  2:79569/251‑1 (MQ=255)                                                                                                                                                                                     ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGAGGCGAGCTGATTTTTGAAGACAAGATAATCGGCTCAAAGTCCg                                                               <  2:474666/251‑1 (MQ=255)                                                                                                                                                                                     atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGTCCAGAACACAGATCTACAACGTAATGCGCTGAACTTTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTc                                                                 <  1:35168‑M2/185‑1 (MQ=255)                                                                                                                                                                                        gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTGTGAAGACAAGATAAGCGTCACAAAGTCCGaaa                                                            <  2:644040/251‑1 (MQ=255)                                                                                                                                                                                         gggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAg                                                           <  1:300317/250‑1 (MQ=255)                                                                                                                                                                                                      cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAGCCAACGTTTATGCAAATGAAGAAATAACAAGCAAATGGCCTCATTACTGCTCTTACCAGGGGGCTTTAACATGCTTATTGGCTATGTACGCGTCTCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGTTAAGCGGCACAAAGTACGAAAGGCAGGTACTGaa                                               >  2:64762/1‑250 (MQ=255)                                                                                                                                                                                                                       cattgagttttgctccttctcccaaaccaaCGTTTATGCAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTGCAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGAATGAAAAAACTGCTCAGGACAtt                             >  2:85117‑M2/31‑251 (MQ=255)                                                                                                                                                                                                                                          tcccaccccAACGTTTATGAAAATGCAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGAATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGAAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGGGACAATCTg          >  1:582515/8‑251 (MQ=255)                                                                                                                                                                                                                                           cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTgg         <  1:205072/251‑1 (MQ=255)                                                                                                                                                                                                                                              aaaCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACAGTATCGGCAGGTGACACTCTGGTTg      <  2:321975/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ccAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCt   <  1:396338/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cAACGTTTATGAAAATTACGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTg  >  2:174467/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        GTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTG  >  NC_000913/1209375‑1209866

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