Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	2,121,524	Δ4 bp	9.6%	coding (28‑31/1395 nt)	wcaJ ←	colanic biosynthesis UDP‑glucose lipid carrier transferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2121528 =	78 (0.970)	8 (0.100)	8/484	NT	9.6%	coding (27/1395 nt)	wcaJ	colanic biosynthesis UDP‑glucose lipid carrier transferase
	?	NC_000913	= 2121523	76 (0.980)					coding (32/1395 nt)	wcaJ	colanic biosynthesis UDP‑glucose lipid carrier transferase

CCGCGTGGAACAGCATTTTAGCCGTGAGGCGCTGGCGGTGGATCGCACCGATGGCATCAGCATGACCTTTGCCGACTGGCGCTTTAACCTGCGCACCTCCAATACCGAACCGGTGGTGCGCCTGAATGTGGAATCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2121807‑2121528 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgagcgaAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCATTCCTCTACATGCACCTGTTGGTGGCGCTGATTACGCTGGTGGTGTTCCAGATGCTGGGCGGCATCACCGATTTTTATCGCTCATGGCGCGGTGTTCGGGCAGCGACAGAATTTGCCCTGTTGCTACAAAACTGGACCTTAAGCGTGATTTTCAGCGCCGGACTGGTGGCGTTCAACAATGATTTCGACA  <  NC_000913/2121523‑2121236                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          CCGCGTGGAACAGCATTTTAGCCGTGAGGCGCTGGCGGTGGATCGCACCGATGGCATCAGCATGACCTTTGCCGACTGGCGCTTTAACCTGCGCACCTCCAATACCGAACCGGTGGTGCGCCTGAATGTGGAATCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATG                                                                                                                                                                                                                                                                                         <  2:545955/289‑1                                                                                    TTAACCTGCGCACCTCCAATACCGAACCGGTGGTGCGCCTGAATGTGGAATCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGA                                                                                                                                                                                                      <  2:641848/289‑1                                                                                                                         CCTGAATGGGGAATCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCATTCCTCTACATGCACCTGTTGGTGGC                                                                                                                                                                      <  1:429432/284‑1                                                                                                                                   GAATCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACAATCATGTTTGCCGGACTAGGGCTGGTTTGCGAAGTAAGCGGACTGTCATTCCTCAACATGCACCTGTTGGTGGGGTTGATTACGCcgt                                                                                                                                                        >  2:429432/1‑285                                                                                                                                      TCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCA                                                                                                                                                                                                >  1:313584/1‑245                                                                                                                                      TCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCA                                                                                                                                                                                                <  2:313584/245‑1                                                                                                                                                        CTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCATTCCTCTACATGCACCTGTTGGTGGCGCTGATTACGCTGGTGGTGTTCCAGATGCTGGGCGG                                                                                                                                  >  1:296795/1‑289                                                                                                                                                                                                     taagtgcgcATCTGCCCGTGCCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGAAAACCAATGCATCTTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCAGCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCATTCCTCTACATGCACCTGTTGGTGGCGCTGATTACGCTGGTGGTGTTCCAGATGCTGGGCGGCATCACCGATTTTTATCGCTCATGGCGCGGTGTTCGGGCAGCGAC                                                                                     <  2:296795/280‑1                                                                                                                                                                                                                                                                                        AAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCATTCCTCTACATGCACCTGTTGGTGGCGCTGATTACGCTGGTGGTGTTCCAGATGCTGGGCGGCATCACCGATTTTTATCGCTCATGGCGCGGTGTTCGGGCAGCGACAGAATTTGCCCTGTTGCTACAAAACTGGACCTTAAGCGTGATTTTCAGCGCCGGACTGGTGGCGTTCAACAATGATTTCGACA  <  1:280623/289‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          CCGCGTGGAACAGCATTTTAGCCGTGAGGCGCTGGCGGTGGATCGCACCGATGGCATCAGCATGACCTTTGCCGACTGGCGCTTTAACCTGCGCACCTCCAATACCGAACCGGTGGTGCGCCTGAATGTGGAATCGCGCGGTGATGTGCCGCTGATGGAAGCGCGAACGCGAACTCTGCTGACGTTGCTGAACGAGTAATGTCGGATCTTCCCTTACCCCACTGCGGGTAAGGGGCTAATAACAGGAACAACGATGACAAATCTAAAAAAGCGCGAGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2121807‑2121528 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgagcgaAAACCAATGCATCGTTAATCTCTATGGTGCAACGCTTTTCAGATATCACCATCATGTTTGCCGGACTATGGCTGGTTTGCGAAGTCAGCGGACTGTCATTCCTCTACATGCACCTGTTGGTGGCGCTGATTACGCTGGTGGTGTTCCAGATGCTGGGCGGCATCACCGATTTTTATCGCTCATGGCGCGGTGTTCGGGCAGCGACAGAATTTGCCCTGTTGCTACAAAACTGGACCTTAAGCGTGATTTTCAGCGCCGGACTGGTGGCGTTCAACAATGATTTCGACA  <  NC_000913/2121523‑2121236

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