Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,399,171	Δ3 bp	7.1%	coding (806‑808/1104 nt)	mreC ←	cell wall structural complex MreBCD transmembrane component MreC

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3399170	52 (0.710)	4 (0.050)	3/410	NT	7.1%	coding (809/1104 nt)	mreC	cell wall structural complex MreBCD transmembrane component MreC
	?	NC_000913	3399174 =	53 (0.730)					coding (805/1104 nt)	mreC	cell wall structural complex MreBCD transmembrane component MreC

TGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3398936‑3399170 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGT  >  NC_000913/3399174‑3399345                                                                                                                                                                                                                                                                                                                                                                                                                         TGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGGTCATTAGCAACACGATGCACCCCTTCCAGCGTCATCGGGTTAGCGCCGTTACAATCTGCCCCCCCCA                                                                                                                                                                           >  2:948544/1‑238        AGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGCAGTTCATTAGCAACACGATACACCTCATCCGGCGTCATCGGGTGAGCCCCGTTACGGACTGCCCCCCCCACCAaccgc                                                                                                                                                                   >  2:770377/1‑234                                                                                                                                              AATACCGGCGGCATCATCTGCACCAGACGTTCATTAGCAACACGATGCCCTCCTGCCGGCGTCATCGGGTTAGCGCCGTTACGATCTG‑CCCCCACAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACC                            <  2:875514/239‑1                                                                                                                                                               CTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCA                        <  1:390151/227‑1                                                                                                                                                                         CGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCGCAACAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCGGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGCAAGCCTTCCGGGAAACAACCGCCCAAACCGGAAGGCACCAGACAACCACCAACACagaccactaa  >  2:390151/1‑229                                                                                                                                                                                                                                                                                                                                                                                                                         TGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3398936‑3399170 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGT  >  NC_000913/3399174‑3399345

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