Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,785,967	G→A	V236V (GTG→GTA)	gpmM →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,785,967	0	G	A	100.0%	144.2 / NA	39	V236V (GTG→GTA)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/14);  total (25/14)

CAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACT  >  NC_000913/3785890‑3786043                                                                              |                                                                             cAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTaaaa                                                                           <  2:731800/81‑1 (MQ=255)     gCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGa                                                                       >  2:731854/1‑81 (MQ=255)        aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCg                                                                    <  1:731801/81‑1 (MQ=255)            ccAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTAt                                                                >  2:731893/1‑81 (MQ=255)             cAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGtatcc                                                               >  2:731858/1‑77 (MQ=255)             cAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTAt                                                                <  2:731803/80‑1 (MQ=255)             cAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATc                                                               <  1:731802/81‑1 (MQ=255)             cAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATc                                                               >  2:731856/1‑81 (MQ=255)             cAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATc                                                               >  1:731868/1‑81 (MQ=255)             cAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATc                                                               <  2:731804/81‑1 (MQ=255)               ggCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCg                                                             <  1:731805/81‑1 (MQ=255)                        cgccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAgg                                                    >  1:731872/1‑81 (MQ=255)                         gccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGt                                                   <  2:731806/81‑1 (MQ=255)                             ttGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                               >  1:731874/1‑81 (MQ=255)                              tGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                               >  2:731890/1‑80 (MQ=255)                              tGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                               >  2:731853/1‑80 (MQ=255)                               gCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                               >  2:731879/1‑79 (MQ=255)                                  ggTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAc                                          >  1:731876/1‑81 (MQ=255)                                   gTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcg                                         >  1:731870/1‑81 (MQ=255)                                      tGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                       <  1:731807/80‑1 (MQ=255)                                        cAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGc                                     >  1:731882/1‑80 (MQ=255)                                            ctgctTATGCT‑‑CGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                              <  1:731808/81‑1 (MQ=255)                                              gctTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATgg                                >  2:731889/1‑79 (MQ=255)                                               ctTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATgg                                >  1:731910/1‑78 (MQ=255)                                               ctTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATgg                                >  2:731881/1‑78 (MQ=255)                                               ctTATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                             <  2:731809/81‑1 (MQ=255)                                                ttATGCTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                            >  1:731900/1‑81 (MQ=255)                                                     cTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCg                       <  2:731810/81‑1 (MQ=255)                                                     cTCGCGACGAAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCg                       <  1:731811/81‑1 (MQ=255)                                                         cgacgaAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg                    >  1:731901/1‑80 (MQ=255)                                                         cgacgaAAATGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg                    >  1:731905/1‑80 (MQ=255)                                                              aaaaTGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGa              >  1:731886/1‑81 (MQ=255)                                                               aaaTGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAt             <  1:731812/81‑1 (MQ=255)                                                                aaTGATGAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtt            >  1:731883/1‑81 (MQ=255)                                                                 atgatgAATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAttt           >  1:731915/1‑81 (MQ=255)                                                                   gatgaATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt          >  1:731897/1‑80 (MQ=255)                                                                    atgaATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTc         <  2:731813/80‑1 (MQ=255)                                                                    atgaATTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa        >  1:731911/1‑81 (MQ=255)                                                                        aTTCGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa    >  2:731907/1‑81 (MQ=255)                                                                           cGTAAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACt  <  2:731814/80‑1 (MQ=255)                                                                              |                                                                             CAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACT  >  NC_000913/3785890‑3786043

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