New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	25 (0.480)	8 (0.170)	8/456	2.6	28.4%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	17 (0.350)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

AATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGA  >  NC_000913/1209389‑1209868                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          aaTACGTAAAATAGGACGTCCTGGATCAATCCCCCCCCCGGCCTCCCATCCACGAATAAACTCACCACGTAAATCAGGCAATTTGTTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAATGCGGCTCCATGGCAATTCACCCAGCCTGTTCGCTAAGTGGCTGAAGTCCACGGATCCGGGACATCAACAGGTAAGGCAGAGCCTTCTCCCATACCAACGTTTATGaagaggtagaaat                                                                                                                                                                                                                                       <  2:467000/251‑12 (MQ=255)         tgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                                                                                                                                >  1:1152‑M2/223‑250 (MQ=255)         tgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                                                                                                                                <  1:495235‑M2/28‑1 (MQ=255)                           cAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCt                                                                                                                                                                                                             >  2:85253/1‑251 (MQ=255)                             gcaacaaaggcattattaaacggaaccgggtatacaactgtagttacatcccccgcagaagcagatgcttatgcggcattcaacccccactgaatgattactgattgatcaataccaccaatcataaccggtatagctatgacaccggtgcgtcctaatgtcaccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                                                                                                                                           >  2:489864‑M2/203‑251 (MQ=255)                                             cccAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                                                                                                                           <  1:341731/251‑1 (MQ=255)                                                       cgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                                                                                                                                 <  2:418310‑M2/75‑1 (MQ=255)                                                       aaTAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCCGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGAAACTGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                                                                                                                                 >  1:466097/1‑251 (MQ=255)                                                       aaTAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCAGTCCTGCTTTTACCAGGGGGATTGAACATGCTTAtt                                                                                                                                                                                 <  2:445308/251‑1 (MQ=255)                                                                             ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATAGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGACAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca                                                                                                                                                            >  1:215221/1‑250 (MQ=255)                                                                               cAATTTAGTTGTCTGATAAGCCTGTGCCAGTTCCGTGTATTCTTCAGCAGAAAAAGCGCCACCAGTGCATTTCAGCCAGCCTCTTGGCGGAGTGGCTGAAGGCCAAGGAACCGGGACACCAACAGGTAATGCAGAGCCGTCTCCCAAACCAACGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATGCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                                                                                                                         <  2:466097/251‑1 (MQ=255)                                                                                       ttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa                                                                                                                                                  >  1:535204/1‑250 (MQ=255)                                                                                         gTCGGAAAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGGAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                                                                                                                                <  2:139847/250‑1 (MQ=255)                                                                                            ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGTATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                                                                                                                                            <  1:426168/251‑1 (MQ=255)                                                                                            ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACTGAACCGGGACACCAACAGGTAATTCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTGATTGGCTATTTACGCGTATCAACAAATGACCAGAACACAg                                                                                                                                            <  1:315570/251‑1 (MQ=255)                                                                                               gcttttgcagcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                                                                                                                                              >  2:436255‑M2/137‑246 (MQ=255)                                                                                                               acccccactgaatgattactgattgctcaataccacctattataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagtttggctccgtctcccaaaccaaCGTTGATGAAAATGAATAAATAACATGCAAATGGCATCATTCCTTTTTTGACCAGGGGGATTTAACATGTTTATTGGCTATGTACGCGTATCATCAAATGACCAGAACACAGATCTACAACGTAATGCGCt                                                                                                                         <  1:242321‑M2/131‑1 (MQ=255)                                                                                                                            agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAgg                                                                                                             >  1:245929/1‑250 (MQ=255)                                                                                                                                 aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa                                                                                                       <  1:56051/251‑1 (MQ=255)                                                                                                                                                                          gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGaaa                                                              <  2:245929/251‑1 (MQ=255)                                                                                                                                                                                gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCgg                                                        >  1:260429‑M2/56‑251 (MQ=255)                                                                                                                                                                                           aaCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGaaaaa                                              >  1:458889/1‑250 (MQ=255)                                                                                                                                                                                                          attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTACTCATACAAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACAtt                               >  1:502463‑M2/30‑250 (MQ=255)                                                                                                                                                                                                            gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATc                            <  2:458889/251‑1 (MQ=255)                                                                                                                                                                                                                   agagCCTTCTACCCAACCAACCTTTATGAAAATGAAGAAATAAAAAGCAAATGGCATAATTCCTGCTTTTACCAGGGGGATCTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACAAGAACAAAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGAGGCACAAATTCCGAAAGGCAGTGACTGAAAACAATGCTCAGGCCATTAtaggaaggg                     >  2:363310/1‑243 (MQ=255)                                                                                                                                                                                                                             cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATTGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTACGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAAATAAGCGGCACAAAGCCCGAAAGGCCGGGGCGGAATAAACTGCTCAGGTCATTATCGGCAGGTGACActct             <  2:215221/249‑1 (MQ=255)                                                                                                                                                                                                                                      aCGTTTATGAAAATGAAGAAATAACAAGCAAATTGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGa  <  2:289196/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          AATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGA  >  NC_000913/1209389‑1209868

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.