New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	13 (0.410)	10 (0.340)	10/454	1.1	49.8%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	8 (0.270)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGA  >  NC_000913/1209357‑1209835                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          ccATGAGCATGAGCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCtt                                                                                                                                                                                                                                       <  1:195048/250‑1 (MQ=255)   aTGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGTCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATGGCAGCCAGCCTGTTGGCGGAGTGGCTGATGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTtctc                                                                                                                                                                                                                                    <  1:284983/251‑1 (MQ=255)     gAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTcc                                                                                                                                                                                                                                   <  1:311865/250‑1 (MQ=255)                     tGCATACCCCTGAATACTTCGAATAGAGCCACCTGTATCAATACCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGAATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCAGAGAGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTAATg                                                                                                                                                                                                                   >  1:98282/1‑250 (MQ=255)                           ccccTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCTTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                                                                                                                                                            <  2:300837/251‑1 (MQ=255)                                                  gACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATAAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGAAATGCAGAGccttctcccaaaccacgattttcgagcgccgtttgaacagtgccatccgatttg                                                                                                                                                                                     <  1:179735‑M2/251‑55 (MQ=255)                                                            aaTCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                                                                                                           <  1:87791/251‑1 (MQ=255)                                                                 cccgcccgTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgg                                                                                                                                                                       <  1:277097‑M2/250‑69 (MQ=255)                                                                    gcccgTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattc                                                                                                                                                                   <  1:40070‑M2/251‑73 (MQ=255)                                                                     cccgTCACCCCAGCCACGAATAAACTCACCACGTAAAACAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTGCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                                                                                                                                  >  2:28476‑M2/1‑178 (MQ=255)                                                                                           accTCACCACGTAAATAAGGCAATTTAGTTGTCGTATACGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGCAAAAGCGGAACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGTACCGGGACACCAACAGGAAATGGAGAGCCTTATCCCAAACCAACGTTTATGAAAATACAGAACTAACAAGCAACTGGCATCATTCCTGCTTTTACCAGGGGGATGTAACAAGCTTATTgga                                                                                                                                             >  2:34885/3‑249 (MQ=255)                                                                                           aaCTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcc                                                                                                                                                                             <  2:218224‑M2/218‑63 (MQ=255)                                                                                           aaCTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcc                                                                                                                                                                             >  1:218224‑M2/1‑156 (MQ=255)                                                                                                                aaTTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtc                                                                                                                       <  2:245330‑M2/251‑117 (MQ=255)                                                                                                                                                                           aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGTGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGAttt                                                              <  2:290563/249‑1 (MQ=255)                                                                                                                                                                           aCCATTGCATTTCAGCCAGCATGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGCAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGAttt                                                              >  1:290563/1‑249 (MQ=255)                                                                                                                                                                                       cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtca                                                 <  1:15088‑M2/250‑187 (MQ=255)                                                                                                                                                                                                      cGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTcc                                 >  1:316732/1‑251 (MQ=255)                                                                                                                                                                                                        gAGTGGCTGAATGCCACGGAACCGAGACACCAACAGGTAATGCAGAGccttctaccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagaaactggtcgtgcctccccttctccaggctggcaccggctgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcctcctccat                                <  2:326580‑M2/250‑204 (MQ=255)                                                                                                                                                                                                             gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAgg                           <  1:37475/250‑1 (MQ=255)                                                                                                                                                                                                                                     ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccggtctgggggtt  <  2:31354‑M2/251‑234 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGA  >  NC_000913/1209357‑1209835

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