New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	15 (0.260)	20 (0.370)	19/450	1.1	61.5%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	11 (0.210)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATC  >  NC_000913/1209358‑1209842                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 cATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTCAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTATCGGATAAGCCTTTGCCAGTTCAGGGTATTCTTCAGCAGAAAAAGCGTCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCACAAGGTAATGCCGAGCCTtct                                                                                                                                                                                                                                            >  2:257336/1‑251 (MQ=255)              ctccGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGACAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGAAACCAACAGGTAATGCAAAGccttctccaaaaccaaggg                                                                                                                                                                                                                               <  2:98589‑M2/251‑20 (MQ=255)                ccGTTGCATACCCCTGAATACTTAAAAGAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACTCCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                                                                                                                             <  1:257336/251‑1 (MQ=255)                      cattccccgtaaTACTTAAAATAGAGCGACCCGTATCAATCCCCCGCCCGTCAGCCCAGCCACGAATAAACTCACCTCGTAAATCAGGCAATTTATTTGTCGGATAAGACTTTGCCAGTTCCGGGTATTCTTCGGCAGAAAAAGCGGCACCAGTGGATTTCAGCCAGCCTGTGGGCGGAGTGGCTGAACGCCGCGGACCCGGGGCACCAACAGGTAAGGCAGAGCCTTCTCCCAAACCaaggtggtcgaga                                                                                                                                                                                                                       <  2:247121/241‑12 (MQ=255)                       aTACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaatgttttcgagag                                                                                                                                                                                                                      >  2:67517‑M2/1‑223 (MQ=255)                                        aaaTAGAGCGACCTGTATCAATCCCCCGACAGTCATACCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTGGTCGGACAAGCCTTTGCCAGTTCCGGGTATGCTTCAGAAGAAAAAGCGGCACCATGTCATTGCAGCCATCCTGTTGGCGGAGTGGCTGAAGGCCACCGTACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGaataaata                                                                                                                                                                                                            <  1:576986/244‑1 (MQ=255)                                              agCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgct                                                                                                                                                                                               <  2:201629‑M2/251‑52 (MQ=255)                                                cGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                                                                                                                             >  1:201629‑M2/1‑198 (MQ=255)                                                     tGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                                                                                                        >  1:32675/1‑251 (MQ=255)                                                      gTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTAACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGACAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc                                                                                                                                                                                       <  1:475878‑M2/251‑60 (MQ=255)                                                             tCCCCCGCCCGTCATCCCAGCCACGAATAAACTTACCACGTAAATCAGGCAATATATTTGTCGGATAAGCCTTTGCCAGCTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGAAGTGGCTGAAGGCCACGGAACCGGCACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaac                                                                                                                                                                                <  1:79487‑M2/251‑67 (MQ=255)                                                                 ccgcccgTCATCCCATCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTTTCGGATTAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                                                                                                            <  1:65294/251‑1 (MQ=255)                                                                          aTCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCATCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                                                                                                   <  2:472764/251‑1 (MQ=255)                                                                                               accacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtg                                                                                                                                              <  1:564281‑M2/251‑101 (MQ=255)                                                                                                          cAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaac                                                                                                                                   <  1:281199‑M2/251‑112 (MQ=255)                                                                                                                    atttGTCGGGTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATGTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCGCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGTCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAGTGGCTATGTACGCGTATCATCAAATGACCAg                                                                                                                         <  2:248134/251‑1 (MQ=255)                                                                                                                       tGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcc                                                                                                                      <  1:24100‑M2/251‑125 (MQ=255)                                                                                                                       tGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                                                                                                      >  2:413438/1‑251 (MQ=255)                                                                                                                             ttaAGACTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCCCGGAACCGGGACACCACCAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaaagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttct                                                                                                                <  2:422865‑M2/250‑131 (MQ=255)                                                                                                                                  ccTTTGCAAGTTCCGCGTATTCTTCAGCAGACAAAGCGGCACCATTGCATTTCAGCCAGACTGTTAGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccagg                                                                                                           <  2:373817‑M2/251‑136 (MQ=255)                                                                                                                                         cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACATCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcac                                                                                                    >  2:270051‑M2/1‑109 (MQ=255)                                                                                                                                                gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg                                                                                              >  2:506446/1‑250 (MQ=255)                                                                                                                                                     ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCACCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaccagcagcgcgcgcagtgcgggaagcaactggtcgtgccgccccttctccaggctggcaccggatgcctccc                                                                                        <  2:188657‑M2/251‑155 (MQ=255)                                                                                                                                                          cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccccg                                                                                   <  1:102509‑M2/251‑160 (MQ=255)                                                                                                                                                                aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGAAGTGGCTGACGGACAAGGAACCGGGACACCAACAGGTCATGCAGAGccttctcccaaaccaaggatttcgagcgccgtttgcacagagccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactgggcgtgccgccccttctccaggcttgcaccggatgcctccaccacgctgcc                                                                              <  2:43071‑M2/250‑165 (MQ=255)                                                                                                                                                                  aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaag                                                                           >  1:342954‑M2/1‑84 (MQ=255)                                                                                                                                                                                gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGTTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgctgtaagtaacgggtcgtgccgccccttctccaggctggcaccggatgcctccacctcgctgcaaagttcctcctgcaaca                                                             >  2:496595‑M2/1‑70 (MQ=255)                                                                                                                                                                                 cATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAACTGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGAc                                                             >  2:494586/1‑250 (MQ=255)                                                                                                                                                                                       agccaaccTGTTGGCGGAGTGGCTGAAGGCCACGGCACCGGGACACCAACAGGTAATGCAGCGccttctcccaaaccaaggctttcgagagccgtttgcacagtgccatccgatttgatatagccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcacagttcctcctgcaacatgtcac                                                       <  2:387161‑M2/250‑188 (MQ=255)                                                                                                                                                                                        gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAg                                                     <  2:32675/251‑1 (MQ=255)                                                                                                                                                                                        gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGcctgctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgccgccaccacgctgcaaagttccttcggcaacatgtcaaag                                                     >  2:564281‑M2/1‑62 (MQ=255)                                                                                                                                                                                                gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGATccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggaggcctccaccacgctgcaaagtgcctcctgcaacatgtcaaagtagtcat                                              >  2:281199‑M2/1‑54 (MQ=255)                                                                                                                                                                                                                                           gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATc  >  1:29497/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 CATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATC  >  NC_000913/1209358‑1209842

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