New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	17 (0.470)	9 (0.260)	8/456	2.1	36.8%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	15 (0.440)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

GTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTG  >  NC_000913/1209375‑1209866                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         gcattgataggtgcgctgacattggcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                                                                                                                                                                     <  2:139639‑M2/6‑1 (MQ=255)          cccTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATCAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTATCCCAAACCAACGTTTATGAAAATga                                                                                                                                                                                                                                            >  1:46148/1‑250 (MQ=255)              gAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCGTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                                                                                                                                                                                                                       <  1:118264/251‑1 (MQ=255)               aaTACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCAGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTGAAGAAAt                                                                                                                                                                                                                                      >  1:353326/1‑251 (MQ=255)                                                 cgccctcactCCCAGCCACGAATAAACACACCACGCAAATCAGGCAATTTATTTTTCGGATAAGCCGTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCATGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAGCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                                                                                                                                    <  2:243734/242‑1 (MQ=255)                                                 caagcattattaaacgaaacagggaatacaactgtatttccatcccccccagaagccgatgcttttgcggcattccaccccaactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccacaccaaCGTTTATGAAAATGAAGAAATAACACGCACATGGCATCATTCCTGATTTTAcc                                                                                                                                                                                                      >  2:23281‑M2/197‑249 (MQ=255)                                                      attaataaacgcaaccgagaattcaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAGGAAGAAATATCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                                                                                                                                               <  2:452806‑M2/60‑1 (MQ=255)                                                                      gggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaatcccaccactcataaccggtctagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                                                                                                                               >  1:289551‑M2/176‑251 (MQ=255)                                                                             caccacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGAATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCTAATGGCATCATTCCTGCTTTTACCAGGTGGAGTTAACATGCTTATTGGCTAt                                                                                                                                                                          <  1:392054/249‑1 (MQ=255)                                                                                                         agatgcttttgagccattacaaccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgcaaacgatgcgtcctaatgcagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAACTGGCATCATTACTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                                                                                                                                             >  2:287976‑M2/141‑250 (MQ=255)                                                                                                             taaaccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGAAAATGGCATCATTCCt                                                                                                                                                                                                               >  2:272286/5‑180 (MQ=255)                                                                                                             tAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                                                                                                                               <  1:272286/180‑1 (MQ=255)                                                                                                                 ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc                                                                                                                                     >  1:80209/1‑250 (MQ=255)                                                                                                                                      gaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACtgtg                                                                                                               <  1:23281‑M2/139‑1 (MQ=255)                                                                                                                                             agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAAACGGGACACCAACAGGTAATGCAGAGCCTTCTCCCACACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGCATGc                                                                                                       >  2:106592/1‑251 (MQ=255)                                                                                                                                                          ataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatatagccgttgctgcattgagttttgctccttctcccaaaacaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑ACCTGTGCAGGATGCGAGCTGATTTTTg                                                                                          >  2:131436‑M2/92‑251 (MQ=255)                                                                                                                                                             ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGa                                                                                         >  2:153034/1‑249 (MQ=255)                                                                                                                                                                       gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATa                                                                               >  2:75396/1‑249 (MQ=255)                                                                                                                                                                           cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCt                                                                                                  <  1:386901‑M2/152‑1 (MQ=255)                                                                                                                                                                           cggtatagctatgaaaacggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGCAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGTTGCGAGCt                                                                                                  >  2:386901‑M2/75‑226 (MQ=255)                                                                                                                                                                            ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGca                                                                        <  1:154105/251‑1 (MQ=255)                                                                                                                                                                                       aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGaa                                                             <  1:153034/251‑1 (MQ=255)                                                                                                                                                                                             tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGAGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCg                                                       >  2:342209/1‑251 (MQ=255)                                                                                                                                                                                              gAAGGCCATGGACACGGGACCACAACAGGTAATGCATAGTCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGTTTGTTCCATGTGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCATAACTCAGATCGTCAACGTACTGC‑TTGTAACTGTGCAGGAGGCGAGCTGTTGTTTGAAGACAAGATAATCTTCTCAAAGTCCGCAAGGcc                                                        <  2:439824/249‑1 (MQ=255)                                                                                                                                                                                                                                          tcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                                                                                                                                           <  1:267908/91‑1 (MQ=255)                                                                                                                                                                                                                                          tcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTg          >  1:394645/1‑251 (MQ=255)                                                                                                                                                                                                                                          tcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGAATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                                                                                                                                           >  2:267908/1‑91 (MQ=255)                                                                                                                                                                                                                                                 ccAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCt   <  1:287976/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTg  >  1:331187/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         GTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTG‑AACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTG  >  NC_000913/1209375‑1209866

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.