New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	16 (0.260)	9 (0.150)	9/458	2.4	36.0%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	17 (0.290)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGG  >  NC_000913/1209373‑1209867                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         ccGTTGCATACCCCTGAATTCTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGACTTTGCCAGTTCCTGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGTGGAGTGGCTGAAGGCCACGGAACCGGGTCACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                                                                                                                                                      <  2:510618/251‑1 (MQ=255)   gcattgataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                                                                                                                                                                                                                      <  1:290097‑M2/5‑1 (MQ=255)      gCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaa                                                                                                                                                                                                                                                 <  1:367558/251‑1 (MQ=255)       cATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAAGCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                                                                                                                                                                                <  2:313914/251‑1 (MQ=255)                        acaaTAGAGCGACCTGTATCAATCCCCCCCCCGTACTCCCAGCCACGAATAACCTCACCACGTAAATCAGGTAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                                                                                                                               <  2:879/249‑1 (MQ=255)                           tcatactttgcaacaacggcaccacaggccttattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggcatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                                                                                                                                                                                                             >  2:265019‑M2/221‑250 (MQ=255)                                           cAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACTATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                                                                                                                             <  1:188368/250‑1 (MQ=255)                                                  acaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTACTGCTTTTAcc                                                                                                                                                                                                      >  2:475396‑M2/198‑250 (MQ=255)                                                    aggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                                                                                                                                   <  1:251838‑M2/56‑1 (MQ=255)                                                           aTCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTCTTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCAAGGAACCGGGACACCAACAGGTAATGCAGAGCATTATCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAACAGGGGGAttt                                                                                                                                                                                            >  1:355132/1‑251 (MQ=255)                                                                   caaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccacggaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                                                                                                                                                                    <  1:117119‑M2/71‑1 (MQ=255)                                                                      ccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                                                                                                                                 >  2:113862‑M2/178‑251 (MQ=255)                                                                                           catcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                                                                                                                            >  1:355502‑M2/157‑251 (MQ=255)                                                                                                                            gTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgc                                                                                                                           >  1:410082/1‑251 (MQ=255)                                                                                                                               ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGATGGCGGAGTGGCTGATGGACACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCAAACCAACGTTTATGAAAATGAAGCACTAACAAGCAAATGGCATCATTCCTGCTTTTACCAGTGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTCATgcgat                                                                                                                        >  1:269417/1‑249 (MQ=255)                                                                                                                                     aTTCTTCAGCAGAAAAAGCGGCACCATTGCAGTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGTAAATGGCATCATTCCTGCTTTTACCAGGGTGATTGAACATGCTTATTGGCTATGTACGCGTATCATCAAATGACCAGAACACAGATCTACATCTTCATGCGCCGAACtg                                                                                                                  <  2:385139/251‑1 (MQ=255)                                                                                                                                         ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                                                                                                                                      >  2:133356/1‑195 (MQ=255)                                                                                                                                         ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                                                                                                                                      <  1:133356/195‑1 (MQ=255)                                                                                                                                                   gattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAg                                                                                                     <  2:40384‑M2/150‑1 (MQ=255)                                                                                                                                                     aacGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGACACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCGAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTCACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAATACAGATCTACAACGTAATGCGCTGACCTGTGCAGGATGCGAGCTg                                                                                                  >  2:850/3‑251 (MQ=255)                                                                                                                                                        agCACCATTGCATCTAAGCCAACCTGTTGGCGGAGTGGCTGGACGACACGGCACCGGAACACCCACAGGAAGTGCAGAGCCTTAACCCAAAACAACGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATACCTACTTTTACCAGGGGGATTTAACCTGCTTATTGGATATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGAtt                                                                                               >  2:625858/2‑251 (MQ=255)                                                                                                                                                                     atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACa                                                                                     >  2:221395‑M2/83‑248 (MQ=255)                                                                                                                                                                      tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGAGTTTTGAAGACAAGa                                                                                  <  2:410082/250‑1 (MQ=255)                                                                                                                                                                                                                    acCAACAGGTAATGCAGAGCCTTCTCCCAAACCGACGTTTAAGAAAATGAATAAATAACAAGCAAATGGCATCATTCCTGCTTTGACCAGGGGGATTTAACATGCTTATTGACTATGTACGCGTATCAACAAATGACCAGAACACAGATCTATAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGATAAGATAAGCGGCACAAAGCACGAAAGGCCGGGACTGAAAAAACTGCTCAgg                                   <  1:610665/251‑1 (MQ=255)                                                                                                                                                                                                                                  cAGAGCCTTCTCCCAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGa                                                                                                                                          >  1:185073/1‑134 (MQ=255)                                                                                                                                                                                                                                  cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGa                                                                                                                                          <  2:185073/134‑1 (MQ=255)                                                                                                                                                                                                                                             cccAAACCAACGTATATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGAACTGAAAAAACTGCTCAGGACATTATCGGCAGGCGTCActcttg          >  1:144099/1‑249 (MQ=255)                                                                                                                                                                                                                                                     aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGACCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTgg  >  2:433779/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         CCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGG  >  NC_000913/1209373‑1209867

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.