New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	45 (0.630)	15 (0.230)	13/458	1.8	35.3%	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

TTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTC  >  NC_000913/1209376‑1209864                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      cattgataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                                                                                                                                                                                                                <  2:28430‑M2/8‑1 (MQ=255)                             gCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCCCCACGTAAATCCGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                                                                                                                                                                     >  1:356819/1‑250 (MQ=255)                                    gTATCAATCGCCCGCCCGCCACCCCAGCCACGAACCAACTCACCACGTAAATCAGACATTTTATTTGTCGCATAAGCCCTTGCCAGTTCCGGGTCTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCCTCATTcc                                                                                                                                                                                                             >  2:411447/1‑251 (MQ=255)                                    acaacggaaacccaggccttatgaaacgcacccgggaataccactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgcatgattagtgattgctcaataccaccaatcataaccggtatagctatgaaaacggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAATTAACAAGCAAATGGCATCATTcc                                                                                                                                                                                                             >  2:263083‑M2/209‑251 (MQ=255)                                                                  accgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgcttgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                                                                                                                                >  1:297310‑M2/179‑250 (MQ=255)                                                                     gggaatacaactgtatttccatccccgccagaggcagatgcttttgcggcattccacccccacggaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccagCGTGTATGAAAATGAAGAAATGAGAAGCAAATGGCATCATTCCGTCTGTTACCAGGGTGATTTAAGAGGCTTattc                                                                                                                                                                            <  2:229529‑M2/76‑2 (MQ=255)                                                                       aaaCTCTCCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                                                                                                                          >  1:564048/1‑251 (MQ=255)                                                                           tCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                                                                                                                                >  1:444898/1‑209 (MQ=255)                                                                           tCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                                                                                                                                <  2:444898/209‑1 (MQ=255)                                                                             aactgtatttccatcccccccagcggcagatgcctttgcggcattccacccccactgcatgattactgattgctcaataccaccaatcataaccggtatagatatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGcttattggaaatgt                                                                                                                                                                     >  2:94555‑M2/168‑244 (MQ=255)                                                                             aactgtatttccatcccacccagaggcagatgcttttgcggccttccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatcaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTATTGGCTATGt                                                                                                                                                                     >  2:257109‑M2/168‑250 (MQ=255)                                                                                 cGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                                                                                                                                <  1:172071/251‑1 (MQ=255)                                                                                                        agatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaacaatcataaccggtatagctctgaaacctgtgcgtcctaatgtagccgttgctgcattgagttttgctctttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCATCAAATGACCAGAacac                                                                                                                                          >  2:342517‑M2/141‑250 (MQ=255)                                                                                                                  gcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgcgttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                                                                                                                               >  2:583845‑M2/131‑251 (MQ=255)                                                                                                                       attccacccccactgaatgcttactgattgctcaataccaccaatcatcaccggtatagctattcatccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcacaaaccacCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTACTGCTTTTACCAGGGGGATTTAACATGCTTATTGGTTTTTTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                                                                                                                          >  1:390235‑M2/126‑251 (MQ=255)                                                                                                                         gTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgc                                                                                                                        >  1:180167/1‑251 (MQ=255)                                                                                                                                  agtcttcAGCAGAAAAAGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCTAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTTCAACGTAATGCGCTGAACt                                                                                                                <  1:214385/248‑1 (MQ=255)                                                                                                                                     gaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc                                                                                                            <  1:580747‑M2/140‑1 (MQ=255)                                                                                                                                              aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCg                                                                                                    >  1:522003/1‑250 (MQ=255)                                                                                                                                                                   ccataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAAAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGAt                                                                              <  2:336438‑M2/170‑1 (MQ=255)                                                                                                                                                                         agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCgg                                                                        >  1:240464/1‑251 (MQ=255)                                                                                                                                                                                       ggggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAATTGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGTTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGaaa                                                          <  2:356819/249‑1 (MQ=255)                                                                                                                                                                                               gcgtcctaatgtagctgctgctgcattgtgttttgctacttctccaaaaccaaCGTTTATGAACATGAAGAAATAACAAGCATCTGGCATCATTCCTGATTTTACCAGGGGGTTTTAACATGCTTATTGGCTATGAACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCggg                                                   >  2:413872‑M2/54‑250 (MQ=255)                                                                                                                                                                                                       gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGTTTATGAAAATGAAGCACTCAAACGCAAATGGCATAATTCCTGCTTTTACCAGGGGGCTTTACCATGCTTATTGGATCTGTACGCGTTTCAACAACTGCCCAGCACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTACGAAAGGCCGGGACTGaaac                                           >  1:267295/1‑249 (MQ=255)                                                                                                                                                                                                                           aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCgg                      <  1:177380/251‑1 (MQ=255)                                                                                                                                                                                                                            gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGc                     >  1:528212‑M2/25‑251 (MQ=255)                                                                                                                                                                                                                                 gctccgtctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTGCAACGTAATGCGCTGAACTGTGCAGGAGGCGAGCTGAGTTGTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGCACTGAAAAAACTGCTCAGGACATGAGCGGCAGGTg                <  2:250858‑M2/232‑1 (MQ=255)                                                                                                                                                                                                                                  ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGAAGGATGCGAGCTGATTTTTGACGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGcaggga                >  1:263218‑M2/19‑248 (MQ=255)                                                                                                                                                                                                                                       tctcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGTGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACAct           <  1:94555/250‑1 (MQ=255)                                                                                                                                                                                                                                          cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTg        <  1:508484/250‑1 (MQ=255)                                                                                                                                                                                                                                          cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAATCCTGCTTTCACCAGGTGGATTTAACATGCTTATTGGCTATGTACGCGGATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCGGATTGTTGAAGACAAGATAAGCGGCAGAAAGTCCGAAAAGCCGGGACTGAAAAAACTGCTCAAGACATTATCTGCAGGGGACACTCTgg       <  2:390235/251‑1 (MQ=255)                                                                                                                                                                                                                                           ccAAACCAACGTTTATGAAAATGACGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGacaacccggg      >  2:266022/1‑244 (MQ=255)                                                                                                                                                                                                                                               aCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCGGATTTTTGAAGACAAGAGAAGCTGCACAAAGTCCGAAAGGCCGGGACGGAAGAAACTGCTCAGGACATTATCGGCAGGTTACACTCAGGTTGTc  <  2:386269/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTC  >  NC_000913/1209376‑1209864

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