Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012967	3,745,391	A→C	F13V (TTT→GTT)	waaG ←	glycosyltransferase family 4 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	3,745,391	0	A	C	100.0%	103.2 / NA	26	F13V (TTT→GTT)	waaG	glycosyltransferase family 4 protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/10);  total (16/10)

ACCCGAACATGATGACCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAAAGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTACCTGTCGTGGAGGAGGGCAG  >  NC_012967/3745321‑3745456                                                                       |                                                                  aCCCGAACATGATGACCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACgg                                                                 <  1:3370129/73‑1 (MQ=255)    cGAACATGATGACCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGaa                                                              >  1:3873705/1‑73 (MQ=255)       aCATGATGACCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAat                                                           >  1:8269447/1‑73 (MQ=255)       aCATGATGACCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAat                                                           <  1:4621185/73‑1 (MQ=255)         atgatgACCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAatat                                                         >  1:12174890/1‑73 (MQ=255)               aCCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTtatat                                                   <  1:8878591/73‑1 (MQ=255)                 cTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATaa                                                 >  1:6271228/1‑73 (MQ=255)                 cTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATaa                                                 <  1:3876342/73‑1 (MQ=255)                     ggcggcGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACaa                                             >  1:7484250/1‑73 (MQ=255)                         gcgACTGTCTGAGCAATAAGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTACATAAACAAAAAg                                         >  1:6905644/1‑73 (MQ=255)                          cgACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGc                                        >  1:12276537/1‑73 (MQ=255)                           gACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCa                                       <  1:1408733/73‑1 (MQ=255)                            aCTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCaa                                      >  1:10510153/1‑73 (MQ=255)                                  tGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCa                                >  1:3266489/1‑73 (MQ=255)                                    aGCAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAACCAAAAAGCAACGAACATa                              >  1:7032974/1‑73 (MQ=255)                                      cAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAac                            >  1:10161815/1‑73 (MQ=255)                                      cAATACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGAAACAATAAAAac                            >  1:5055503/1‑73 (MQ=255)                                       aaTACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAacg                           >  1:1462309/1‑73 (MQ=255)                                         tACGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAacgac                         >  1:138941/1‑73 (MQ=255)                                          aCGCATAAAATCGCGCTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAacgacg                        <  1:10527678/73‑1 (MQ=255)                                                        gcTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTACCTGTCGTggag          <  1:3873321/73‑1 (MQ=255)                                                         cTGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTACCTGTCGTggagg         >  1:10017231/1‑73 (MQ=255)                                                          tGCAAACCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTACCTGTCGTggagga        <  1:6567081/73‑1 (MQ=255)                                                           gCAATCCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTTCCTGTCGTGgaggag       <  1:9703325/73‑1 (MQ=255)                                                             aaaCCGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTACCTGTCGTGGAGGAggg     <  1:1980555/73‑1 (MQ=255)                                                                ccGCCAACGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTACCTGTCGTGGAGGAGGGCAg  >  1:6407549/1‑73 (MQ=255)                                                                       |                                                                  ACCCGAACATGATGACCTCGGGCGGCGACTGTCTGAGCAATACGCATAAAATCGCGCTGCAAACCGCCAAAGGGAAAATATTTATATAAACAAAAAGCAACGATCATAACGACGTACCTGTCGTGGAGGAGGGCAG  >  NC_012967/3745321‑3745456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑