Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,971,307	G→A	L16L (CTG→CTA)	wecC →	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,971,307	0	G	A	100.0%	63.3 / NA	20	L16L (CTG→CTA)	wecC	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/5);  total (15/5)

AGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGC  >  NC_000913/3971173‑3971441                                                                                                                                       |                                                                                                                                       aGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCaa                                                                                                                                    >  2:597306/1‑139 (MQ=255)                          tATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCaa                                                                                                                                    >  2:523790/1‑114 (MQ=255)                          tATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCaa                                                                                                                                    <  1:523790/114‑1 (MQ=255)                          tATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCa                                                                                                           >  1:181776/1‑139 (MQ=255)                          tATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCa                                                                                                           >  1:323983/1‑139 (MQ=255)                                 aTGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCGGCGTTTGCCTCACGGCAAAAACAg                                                                                                    >  2:343323/1‑139 (MQ=255)                                 aTGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAg                                                                                                    >  2:329450/1‑139 (MQ=255)                                                 tcGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGata                                                                                    >  2:280380/1‑139 (MQ=255)                                                                                    cACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGc                                                 <  2:181776/139‑1 (MQ=255)                                                                                             ttttGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCat                                        >  1:166308/1‑139 (MQ=255)                                                                                              tttGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCg                                                                      >  1:537985/1‑108 (MQ=255)                                                                                              tttGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCg                                                                      <  2:537985/108‑1 (MQ=255)                                                                                                tGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATAtc                                     >  1:512269/1‑139 (MQ=255)                                                                                                     ccATTTCTGTTATCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGa                                >  1:587111/1‑139 (MQ=255)                                                                                                                 tCGGACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGc                    <  2:323983/139‑1 (MQ=255)                                                                                                                   ggACTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGa                  >  2:105791/1‑139 (MQ=255)                                                                                                                      cTGGGTTATATCGGGCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAgtg               <  1:329450/139‑1 (MQ=255)                                                                                                                        gggTTATATCGGGCTAACAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTgta             >  1:116101/1‑139 (MQ=255)                                                                                                                                   ggCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGc  >  2:164564/1‑139 (MQ=255)                                                                                                                                   ggCTACCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGc  >  1:182587/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       AGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCATTTCTGTTATCGGACTGGGTTATATCGGGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGC  >  NC_000913/3971173‑3971441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑