Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,399,169	Δ3 bp	6.3%	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	= 3399168	59 (0.750)	4 (0.050)	3/410	NT	6.3%	coding (809/1104 nt)	mreC	cell wall structural complex MreBCD transmembrane component MreC
	?	NC_000913	3399172 =	60 (0.760)					coding (805/1104 nt)	mreC	cell wall structural complex MreBCD transmembrane component MreC

TGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3398934‑3399168 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGT  >  NC_000913/3399172‑3399343                                                                                                                                                                                                                                                                                                                                                                                                                         TGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGGTCATTAGCAACACGATGCACCCCTTCCAGCGTCATCGGGTTAGCGCCGTTACAATCTGCCCCCCCCA                                                                                                                                                                           >  1:948544/1‑238        AGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGCAGTTCATTAGCAACACGATACACCTCATCCGGCGTCATCGGGTGAGCCCCGTTACGGACTGCCCCCCCCACCAaccgc                                                                                                                                                                   >  1:770377/1‑234                                                                                                                                              AATACCGGCGGCATCATCTGCACCAGACGTTCATTAGCAACACGATGCCCTCCTGCCGGCGTCATCGGGTTAGCGCCGTTACGATCTG‑CCCCCACAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACC                            <  1:875514/239‑1                                                                                                                                                               CTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCA                        <  2:390151/227‑1                                                                                                                                                                         CGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCGCAACAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCGGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGCAAGCCTTCCGGGAAACAACCGCCCAAACCGGAAGGCACCAGACAACCACCAACACagaccactaa  >  1:390151/1‑229                                                                                                                                                                                                                                                                                                                                                                                                                         TGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3398934‑3399168 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGT  >  NC_000913/3399172‑3399343

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.