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

ACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTA  >  NC_000913/3790871‑3791137
                                                                                                                                     |                                                                                                                                     
aCAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAAc                                                                                                                                  <  1:299585/139‑1 (MQ=255)
    cgcgGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCt                                                                                                                              >  1:461097/1‑139 (MQ=255)
     gcgGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTa                                                                                                                             <  1:65855/139‑1 (MQ=255)
     gcgGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTa                                                                                                                             <  2:1171090/139‑1 (MQ=255)
                      aCCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                            <  2:461097/139‑1 (MQ=255)
                            gggATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTa                                                                                                      >  1:824602/1‑139 (MQ=255)
                            gggATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTa                                                                                                      >  2:317505/1‑139 (MQ=255)
                            gggATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTa                                                                                                      <  2:910447/139‑1 (MQ=255)
                               aTCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaa                                                                                                             <  2:162647/129‑1 (MQ=255)
                               aTCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaa                                                                                                             >  1:162647/1‑129 (MQ=255)
                               aTCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGa                                                                                                   >  2:1112109/1‑139 (MQ=255)
                                        tGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGa                                                                                          >  2:1075338/1‑139 (MQ=255)
                                                   cATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCa                                                                               <  2:824602/139‑1 (MQ=255)
                                                     tCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAAt                                                                             <  2:28282/139‑1 (MQ=255)
                                                            tatTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                            >  2:484197/1‑101 (MQ=255)
                                                            tatTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACaaa                                                                                                            <  1:484197/101‑1 (MQ=255)
                                                            tatTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACcaca                                                                      >  2:960059/1‑139 (MQ=255)
                                                             atTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAg                                                                     <  1:960059/139‑1 (MQ=255)
                                                             atTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAg                                                                     <  1:1075338/139‑1 (MQ=255)
                                                                   ggggATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGata                                                               >  2:507342/1‑139 (MQ=255)
                                                                   cgggATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGcc                                                                                                                                     >  1:736271/2‑69 (MQ=255)
                                                                                gCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAg                                                  >  2:572269/1‑139 (MQ=255)
                                                                                  gTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCgtgt                                                                                                                       <  2:233257/68‑1 (MQ=255)
                                                                                  gTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCgtgt                                                                                                                       >  1:233257/1‑68 (MQ=255)
                                                                                           cGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACACATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCg                                       >  2:732656/1‑139 (MQ=255)
                                                                                                  aCCAGATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCt                                >  1:564831/1‑139 (MQ=255)
                                                                                                      gATATTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGaa                            <  1:732656/139‑1 (MQ=255)
                                                                                                         atTCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGcc                         >  2:19352/1‑139 (MQ=255)
                                                                                                           tCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGcc                                     <  1:91481/125‑1 (MQ=255)
                                                                                                           tCGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGcc                                     >  2:91481/1‑125 (MQ=255)
                                                                                                            cGGTAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCg                                                                                  <  1:756680/79‑1 (MQ=255)
                                                                                                            cGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCg                                                                                  >  2:756680/1‑79 (MQ=255)
                                                                                                            cGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCtt                                                        <  1:1063966/105‑1 (MQ=255)
                                                                                                            cGGCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCtt                                                        >  2:1063966/1‑105 (MQ=255)
                                                                                                              gCAAAGCAGCCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCa                    >  1:324957/1‑139 (MQ=255)
                                                                                                                   gcagcCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTcc               <  1:19352/139‑1 (MQ=255)
                                                                                                                    cagcCCGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCt              >  2:486480/1‑139 (MQ=255)
                                                                                                                       cccGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGa                                                                                          >  2:183423/1‑60 (MQ=255)
                                                                                                                       cccGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGa                                                                                          <  1:183423/60‑1 (MQ=255)
                                                                                                                         cGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCa         <  1:317505/139‑1 (MQ=255)
                                                                                                                         cGGGCGATTAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCa         <  1:486480/139‑1 (MQ=255)
                                                                                                                                ttAACGCCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAAtaccta  >  1:560322/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                     
ACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTA  >  NC_000913/3790871‑3791137

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