New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	3 (0.070)	8 (0.200)	7/460	1.6	73.4%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	3 (0.080)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

TAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAG  >  NC_000913/1207572‑1208035                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       tAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTAATTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATACGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                                                                                                                                                                                                       >  1:361709‑M1/1‑234 (MQ=255)                                                                  aaCTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGCCCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgct                                                                                                                                                      >  1:284036‑M1/1‑169 (MQ=255)                                                                       tGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttttt                                                                                                                                                  <  2:128028‑M1/250‑87 (MQ=255)                                                                                                               ggcggcCCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgaca                                                                                                         <  2:260786‑M1/251‑128 (MQ=255)                                                                                                               ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGcc                                                                                                         <  2:379309/251‑1 (MQ=255)                                                                                                                                    ccGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgt                                                                                    <  1:88578‑M1/251‑149 (MQ=255)                                                                                                                                                                                                  aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttgtctgctgaagaatacccggaactgtcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgagacaggtcgctctattttaat                      <  2:361709‑M1/251‑211 (MQ=255)                                                                                                                                                                                                      atttcttcatgttcataaaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCGGCAACGTCTACATTAGGACGCACCGTTTTCATAGCTATACCGTTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGTGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATGCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCTCATTAGGACAATGTCAgc                  <  2:257199/235‑1 (MQ=255)                                                                                                                                                                                                        cagacgtcactcgaaaaccttGGTTTGGGAGAAGGcgctgcattacctgtggtggtcccggttccgtggccgtcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgcttaagaatacccggaactggcaaaggcttatccgacaaataaatggcctgatttacgtggtgagtttattcgttgctgggatgacgggcgggggattgatacggttcgctcgattttaagtattca                <  2:17403‑M1/242‑217 (MQ=255)                                                                                                                                                                                                                 ttcataaaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGATTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGACTGTGTTGCCGTTGTTGAAAATTATGATAATGTCAGCGCACCTATCa        >  2:188506/6‑250 (MQ=255)                                                                                                                                                                                                                      aaaCCTTCGTTTGGGAGAAGGccctgaattccctgttggtgtcccggttcagtggccttcagccactccgccaaaaggctggctgcaatgacatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggagagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagaattaaggggagtgcaacgg  >  2:177886‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       TAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAG  >  NC_000913/1207572‑1208035

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