New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	9 (0.230)	15 (0.410)	15/388	1.3	52.3%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	19 (0.520)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

AGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGC  >  NC_000913/1209401‑1209898                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                        acaccGACCTGTAACAATCCCCCGCCCGTCATCCCAGCCACGAAAAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTAAGAAGAAAAAGCGGCACGATTGCATTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGAAACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                                                                                                                                                                        <  1:66857/280‑1 (MQ=255)             aacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaatcggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                                                                                                                                                            <  1:17434‑M2/78‑1 (MQ=255)                  tCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGCTGGCGGAGTGGCTGAAGGCCTCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                                                                                                                                                                                                                            >  1:359139/1‑263 (MQ=255)                                                   cacatgtctctccctcccccccacgggcagatgcattttctgcattccacccccaacgaatgattactgattgctcaataccaccaatcataacgggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                                                                                                                                                                     <  2:88100‑M2/101‑1 (MQ=255)                                                     ccacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTGCGGGTATTCTTCAGCAGAAAAAGCAGCACCATTGCATTTCAGCCAGCCTGTTTGCGGACTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCAAACCAACGt                                                                                                                                                                                                                                                                                       >  1:627319/1‑169 (MQ=255)                                                     ccacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc                                                                                                                                                                    <  1:455362/284‑1 (MQ=255)                                                          tAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTGAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                                                                                                                                                               <  1:634212/284‑1 (MQ=255)                                                          tAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                                                                                                                                                               <  1:199032/284‑1 (MQ=255)                                                          atttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                                                                                                                                                              >  1:88100‑M2/162‑269 (MQ=255)                                                          attaccatccccccccgaggcagatgcttttgccgcactccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgacaccggtgcgtcctaatgtagccgttgctgcattgagtgttgctccttctcccaaaccaaCGTTTATGAAAATGAAAAAACAACAAGCAAATTGCATTATTCCTGCTTTTACCAGGAGGATTTAACCTGCTTAATGGCTATGTACGCGTAACAACAAATGACCAGAACACAGATCTACAACGt                                                                                                                                                               <  1:395394‑M2/123‑1 (MQ=255)                                                                   cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCAGGGTATTCTTCAGCAGAAAAAGCGGCGCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCCCTCAAACCAACGTGTATGAAAATGAAGAAATAACCAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTTTTGTCTCTGTACGCGTCACAACAAATGACCAGAac                                                                                                                                                                              >  1:283798/1‑260 (MQ=255)                                                                                  tgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAAACAGAACTAAAACGTAATGCGCTGAACTGTTCAGGATGc                                                                                                                                       >  1:288553‑M2/138‑284 (MQ=255)                                                                                           tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                                                                                                                                                               <  1:63720/251‑1 (MQ=255)                                                                                                  tCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGTTTTCACATGCTTATTGGCTATGTACGCCTATCAACAAATGACAAGAATCCAGATCGAAAAcaggaagtgcnn                                                                                                                                                      >  2:63720/1‑242 (MQ=255)                                                                                                          ctgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATATCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATTAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATa                                                                                                               >  1:491598‑M2/114‑284 (MQ=255)                                                                                                                                ataccaccactcataaccggtatagctaggaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaacccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAACAGGGGGATTTAACATGCTTATTGGCTATGTACGAGTATCAACAAATGACCAGACCACAGAGCTACAACGTAATgc                                                                                                                                                          <  1:621602‑M2/128‑1 (MQ=255)                                                                                                                                 taccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                                                                                                                                                                   <  2:292374‑M2/103‑1 (MQ=255)                                                                                                                                 taccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                                                                                                                                                                   >  1:292374‑M2/91‑193 (MQ=255)                                                                                                                                     gCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCAACGGAACCGGGACACCAACAGGTAATGCAAAGCCTTCTCCCAAACCAACGTTTATGAAAAGGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCGGGGGGATTTAACATGTTTATTGGTTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                                                                                                                                                      <  1:300256/202‑1 (MQ=255)                                                                                                                                                                      aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                                                                                                                                                               <  1:323198/176‑1 (MQ=255)                                                                                                                                                                             ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                                                                                                                                                              >  1:425513/1‑122 (MQ=255)                                                                                                                                                                             ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                                                                                                                                                              <  2:425513/122‑1 (MQ=255)                                                                                                                                                                             ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCACCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGCCCAGCACACAGACCTACACCGTAATccgc                                                                                                                                                        >  2:323198/1‑176 (MQ=255)                                                                                                                                                                                aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                                                                                                                                                       <  1:629/126‑1 (MQ=255)                                                                                                                                                                                aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                                                                                                                                                       >  2:629/1‑126 (MQ=255)                                                                                                                                                                                aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGCCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                                                                                                                                                             >  1:276831/1‑168 (MQ=255)                                                                                                                                                                                     cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                                                                                                                                                                 <  2:392311‑M2/41‑1 (MQ=255)                                                                                                                                                                                     cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                                                                                                                                                                 >  1:392311‑M2/39‑79 (MQ=255)                                                                                                                                                                                         gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                                                                                                                                                                      <  2:351574‑M2/36‑1 (MQ=255)                                                                                                                                                                                         gctgcattgagttttgctccttctcccaaaacaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                                                                                                                                                                      >  1:351574‑M2/35‑70 (MQ=255)                                                                                                                                                                                          cggtattgagttttgctccttctcccacaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTAGCTATGTACGCATATCAACAAATGACCAGAACACAGATCTACAACGTAATTCGCTGATCTGTGCAGGATGCGCGCTGATTTTTGAAGACAAGATAAGCTGCACAAAGTCCGATAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGaa                               <  1:380021‑M2/251‑1 (MQ=255)                                                                                                                                                                                            gcattgagttatgctactgctcccaaaccaaTGTTTATGAAAATGAAGAAATAACAAGAATATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAACGACCAGAACAAAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAg                                                                   <  1:270313‑M2/215‑1 (MQ=255)                                                                                                                                                                                                     gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                                                                                                                                                               >  1:474206/1‑145 (MQ=255)                                                                                                                                                                                                     gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGt                                                                                                                                                               <  2:474206/145‑1 (MQ=255)                                                                                                                                                                                                              tttcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTCGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTCAACTGTGCATGATGCGAGCTGATTTTTGAAGACAAGATAATCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGc              <  1:573988‑M2/268‑1 (MQ=255)                                                                                                                                                                                                               ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGCTTTTACCAGGGGGATTTAACTTGCTTATTGGCTATGTACGCGTAACAACAAATGAACAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGCTGCGAGCTGATTTTTGAAg                                                                                                                       <  1:622844‑M2/163‑1 (MQ=255)                                                                                                                                                                                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGC‑TTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCTGAACGTACGAACGGCCGGGACTGAAAAAGCTGATAAGGACATTAAGGGAAAGTGACCGCCTGGTGTTAtggaagctgtaacgggggggggggggtggggggg  >  2:457787‑M2/6‑271 (MQ=255)                                                                                                                                                                                                                           cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                                                                                                                                                                >  2:328182/1‑42 (MQ=255)                                                                                                                                                                                                                           cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                                                                                                                                                                  <  2:387517/72‑1 (MQ=255)                                                                                                                                                                                                                           cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                                                                                                                                                                  >  1:387517/1‑72 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                        AGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGC  >  NC_000913/1209401‑1209898

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