Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,968,671	Δ9 bp	5.5%	coding (756‑764/1104 nt)	wecA →	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3968670	72 (0.910)	4 (0.050)	4/514	NT	5.5%	coding (755/1104 nt)	wecA	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
	?	NC_000913	3968680 =	65 (0.820)					coding (765/1104 nt)	wecA	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

TTGATGGCTTGCTGGGCGGGTTGTCCTGCGTCTCGTTTGCAGCAATCGGTATGATTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3968388‑3968670 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCC  >  NC_000913/3968680‑3968870                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            TTGATGGCTTGCTGGGCGGGTTGTCCTGCGTCTCGTTTGCAGCAATCGGTATGATTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCACGGTTACCGGTTTGGGGATAATCGCAATGGA                                                                                                                                                                                           >  2:140935/1‑289                                            AATCGGTATGATTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGC                                                                                                                                                <  1:140935/289‑1                                                                                                                                  TGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCT                                                          >  1:633380/1‑289                                                                                                                                                                                             ACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGCAATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCC  <  2:633380/286‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            TTGATGGCTTGCTGGGCGGGTTGTCCTGCGTCTCGTTTGCAGCAATCGGTATGATTTTGTGGTTCGACGGGCAAACCAGCCTCGCAATCTGGTGCTTTGCGATGATCGCCGCCATCCTGCCATACATCATGCTTAACCTTGGTATCCTGGGTCGCCGCTACAAAGTCTTTATGGGTGATGCGGGCAGTACGCTGATTGGTTTTACCGTTATCTGGATCCTGCTCGAAACGACCCAGGGCAAAACCCATCCCATCAGCCCGGTTACCGCTTTGTGGATAATCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3968388‑3968670 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AATGGATATGGTGGCGATTATGTACCGTCGCCTGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCC  >  NC_000913/3968680‑3968870

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