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_0009133,791,0040AG100.0% 156.6 / NA 42G114G (GGT→GGC) tdhL‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

GACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTACC  >  NC_000913/3790870‑3791139
                                                                                                                                      |                                                                                                                                       
gACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCaa                                                                                                                                     <  2:1099305/139‑1 (MQ=255)
   aGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGc                                                                                                                                  >  2:294944/1‑139 (MQ=255)
      gcgGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTa                                                                                                                               <  1:823656/139‑1 (MQ=255)
                   cGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGca                                                                                                                  <  1:748086/139‑1 (MQ=255)
                      tACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaa                                                                                                               >  1:468066/1‑139 (MQ=255)
                       aCCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGTGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                              <  1:682218/139‑1 (MQ=255)
                             gggATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTa                                                                                                        >  1:807027/1‑139 (MQ=255)
                                aTCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGt                                                                                                         <  1:930466/135‑1 (MQ=255)
                                aTCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGt                                                                                                         >  2:930466/1‑135 (MQ=255)
                                          gCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGAc                                                                                           >  1:100213/1‑139 (MQ=255)
                                                       cGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATg                                                                              >  1:567229/1‑139 (MQ=255)
                                                            atatTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACcac                                                                         >  1:549105/1‑139 (MQ=255)
                                                            atatTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACcac                                                                         <  2:100213/139‑1 (MQ=255)
                                                                    ggggATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                              >  1:928253/1‑94 (MQ=255)
                                                                    ggggATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                              <  2:928253/94‑1 (MQ=255)
                                                                    ggggATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGata                                                                 >  2:1038713/1‑139 (MQ=255)
                                                                    ggggATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGata                                                                 >  2:277336/1‑139 (MQ=255)
                                                                         ttttGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGAc                                                                            >  1:874477/1‑123 (MQ=255)
                                                                         ttttGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGcc                                                            <  2:567229/139‑1 (MQ=255)
                                                                         ttttGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGcc                                                            <  1:639970/139‑1 (MQ=255)
                                                                         ttttGAAGGCGTTGAATGCCGGGAACACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGAc                                                                            <  2:874477/123‑1 (MQ=255)
                                                                                    ttGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCa                                                                                 >  2:355642/1‑107 (MQ=255)
                                                                                    ttGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCa                                                                                 <  1:355642/107‑1 (MQ=255)
                                                                                            cGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGcc                                                            >  1:36243/1‑120 (MQ=255)
                                                                                            cGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGcc                                                            <  2:36243/120‑1 (MQ=255)
                                                                                                   aCCTGATATTCGGCAAAGCAGCCCGGGCGATTAACGcc                                                                                                                                       >  2:1240242/1‑38 (MQ=255)
                                                                                                   aCCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGcc                                                                                                                                       <  1:1240242/38‑1 (MQ=255)
                                                                                                    ccAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                              <  1:130737/62‑1 (MQ=255)
                                                                                                    ccAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                              >  2:130737/1‑62 (MQ=255)
                                                                                                        atatTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAg                             <  1:1117343/139‑1 (MQ=255)
                                                                                                         tatTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGc                            >  1:646280/1‑139 (MQ=255)
                                                                                                         tatTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGc                            <  2:468066/139‑1 (MQ=255)
                                                                                                          atTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACg                                                                                                      <  1:865044/64‑1 (MQ=255)
                                                                                                          atTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACg                                                                                                      >  2:865044/1‑64 (MQ=255)
                                                                                                          atTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGcc                           >  1:196759/1‑139 (MQ=255)
                                                                                                            tCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGacc                                                                                                   >  1:918359/1‑65 (MQ=255)
                                                                                                            tCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGacc                                                                                                   <  2:918359/65‑1 (MQ=255)
                                                                                                               gCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCa                      >  1:931260/1‑139 (MQ=255)
                                                                                                                   agcagcCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTc                  <  1:294944/139‑1 (MQ=255)
                                                                                                                        cccGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGAc             <  1:437749/139‑1 (MQ=255)
                                                                                                                           gggCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCaa          <  2:807027/139‑1 (MQ=255)
                                                                                                                               gATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAAtacc      <  1:277336/139‑1 (MQ=255)
                                                                                                                                   aacgccaacgccTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAAtacctacc  <  2:644420/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTACC  >  NC_000913/3790870‑3791139

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: ‑