New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	12 (0.350)	8 (0.250)	8/454	1.6	51.3%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	4 (0.130)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAAC  >  NC_000913/1209357‑1209826                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 ccATGAGCATGATCCTCCGTTGCATACCCCTGAATACGTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCATAGCCtt                                                                                                                                                                                                                              <  2:350598/250‑1 (MQ=255)                                               agCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                                                                                                                               <  1:124079/251‑1 (MQ=255)                                                  gACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGACTAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAACATGAAGAAATAACAAGCAAATGGcatc                                                                                                                                                                            >  2:264445/1‑251 (MQ=255)                                                                                             cTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcga                                                                                                                                  <  1:318631‑M2/250‑97 (MQ=255)                                                                                              ccaccacgTAAATCAGGCAATTTACTTGGCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAAGGCAGAGccttctcccaaaccaagggttgcgagagccgtttgcacagtgccatccgatgtgacatcgccaatcggattcttgcggcgtaacagcagctcgcgaag                                                                                                                                >  2:122258‑M2/2‑153 (MQ=255)                                                                                                       aaaTCAGGCAATTTATTTGTCGGATAATCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAATCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCATCGTTTATGAAAAGGAAGAAATAACAAGCAAATGGCGTCATTCCGGCTTTGACCGGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                                                                                       <  1:333195/251‑1 (MQ=255)                                                                                                                                         ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggc                                                                                      <  1:249707‑M2/250‑141 (MQ=255)                                                                                                                                                                                                 gTTGGCGGAGTGGCTGAAGGCCACCGAACCGGGCCACCAACAGGTAATGCAGAGccttatcccaacccccggttttcgagagccgtttgcacagtgccctccgatttgatatcgccaaccggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgggccgccccttctccaggctggcaccggatgccaccaccacgctgcaaagttcatcctgcaacatgtcaaagaagtcatc                             <  2:351116‑M2/251‑198 (MQ=255)                                                                                                                                                                                                           tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatcc                         <  1:357925‑M2/245‑202 (MQ=255)                                                                                                                                                                                                                gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCg               <  1:58152/250‑1 (MQ=255)                                                                                                                                                                                                                gAAGGCCACAGAACCTGTACACCAATAGGTAATGCAGTGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtgactggtcgtgccgcgcctgctccaggctggcaccggatgcctccaccacgccgcaaaggtcctcctgcaacgtgtcaaagtaatcttgatccagatcggtg                >  2:235848‑M2/1‑39 (MQ=255)                                                                                                                                                                                                                       aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcttgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaatttcctcctgcaacatgtcaaagtagtcatcatccagagcggtggcaggcgtg       >  1:113601‑M2/1‑32 (MQ=255)                                                                                                                                                                                                                            aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccggt  <  1:21490‑M2/251‑225 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAAC  >  NC_000913/1209357‑1209826

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.