Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,791,004	A→G	G114G (GGT→GGC)	tdh ←	L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,791,004	0	A	G	100.0%	124.8 / NA	33	G114G (GGT→GGC)	tdh	L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/16);  total (17/16)

ACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTC  >  NC_000913/3790868‑3791123                                                                                                                                         |                                                                                                                        aCGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGcc                                                                                                                       <  1:100532/139‑1 (MQ=255)    aCAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAAc                                                                                                                    <  2:601692/139‑1 (MQ=255)            gTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCg                                                                                                            <  2:194384/139‑1 (MQ=255)                    ggCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCgg                                                                                                    >  2:6477/1‑139 (MQ=255)                       gTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGcac                                                                                                 >  2:600114/1‑139 (MQ=255)                       gTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGcac                                                                                                 >  2:27906/1‑139 (MQ=255)                         tACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaa                                                                                               >  2:68195/1‑139 (MQ=255)                           ccGAAGGGATCAAAAATTGCATCTAAGTCATCGGAAATATAGTCGGGGATTTTGAAGCCGTTGCATGCCGGGATCACCAGCTATTCGGCAAA‑CAGCCCGGGCGCTTAACGCCAACGCCTATCGTGTAGCGGCACAAATg                                                                                            >  1:33321/1‑139 (MQ=255)                                gggATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTa                                                                                        >  2:578829/1‑139 (MQ=255)                                  gATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACg                                                                                      <  1:578829/139‑1 (MQ=255)                                   aTCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGa                                                                                     >  1:38276/1‑139 (MQ=255)                                    tCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGac                                                                                    <  1:27906/139‑1 (MQ=255)                                                            gAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGAc                                                            <  2:449679/139‑1 (MQ=255)                                                                 atTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAg                                                       <  2:567078/139‑1 (MQ=255)                                                                       ggggATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGata                                                 <  1:6477/139‑1 (MQ=255)                                                                                       ttGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaa                                                                                               >  2:511071/1‑77 (MQ=255)                                                                                       ttGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaa                                                                                               <  1:511071/77‑1 (MQ=255)                                                                                       ttGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGAc                                                            >  1:61466/1‑112 (MQ=255)                                                                                       ttGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGAc                                                            <  2:61466/112‑1 (MQ=255)                                                                                               cGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCg                         >  2:424127/1‑139 (MQ=255)                                                                                                 tgATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGcc                       <  2:38276/138‑1 (MQ=255)                                                                                                        cAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTg                <  1:235991/139‑1 (MQ=255)                                                                                                             atTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACcac                                                         >  2:459125/1‑93 (MQ=255)                                                                                                             atTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACcac                                                         <  1:459125/93‑1 (MQ=255)                                                                                                               tCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCtt         >  2:79343/1‑139 (MQ=255)                                                                                                                 ggCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTc       >  2:195307/1‑139 (MQ=255)                                                                                                                  gCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTg                >  2:59847/1‑129 (MQ=255)                                                                                                                  gCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTg                <  1:59847/129‑1 (MQ=255)                                                                                                                  gCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCa      >  2:466754/1‑139 (MQ=255)                                                                                                                     aaGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACtt   >  2:553455/1‑139 (MQ=255)                                                                                                                      agcagcCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGATCACCACGACAGTTGCGGCAATGACcac                                                         <  2:395275/84‑1 (MQ=255)                                                                                                                      agcagcCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTc  <  1:466754/139‑1 (MQ=255)                                                                                                                      agcagcCCGGGCGATTAACGCCAACGCCAATCGTGTTGCGGCACAAATGCGTACGACCACCACGACAGTTGCGGCAATGACcac                                                         >  1:395275/1‑84 (MQ=255)                                                                                                                                         |                                                                                                                        ACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTC  >  NC_000913/3790868‑3791123

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