Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,983,191 C→T G23D (GGT→GAT)  kduD ← 2‑dehydro‑3‑deoxy‑D‑gluconate 5‑dehydrogenase, KDG oxidoreductase, 20‑ketosteroid reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,983,1910CT100.0% 134.4 / NA 43G23D (GGT→GAT) kduD2‑dehydro‑3‑deoxy‑D‑gluconate 5‑dehydrogenase, KDG oxidoreductase, 20‑ketosteroid reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/28);  total (15/28)

TTCGCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAACTAGCGCAAATCTTTAACGGCCACATGGTCCATAT  >  NC_000913/2983057‑2983322
                                                                                                                                      |                                                                                                                                   
ttCGCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATccag                                                                                                                                 <  2:363491/139‑1 (MQ=255)
  cGCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATccagtc                                                                                                                               <  1:192538/139‑1 (MQ=255)
  cGCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATccagtc                                                                                                                               <  2:801788/139‑1 (MQ=255)
   gCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATccagtcc                                                                                                                              <  2:355188/139‑1 (MQ=255)
   gCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATccagtcc                                                                                                                              <  1:491419/139‑1 (MQ=255)
    cAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATccagtcca                                                                                                                             <  1:70503/139‑1 (MQ=255)
    cAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATccagtcca                                                                                                                             <  2:372093/139‑1 (MQ=255)
                 aGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAg                                                                                                                >  1:599573/1‑139 (MQ=255)
                         aaaaCGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGAcc                                                                                                        <  2:245806/139‑1 (MQ=255)
                                               aCCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGaaa                                                                                  <  1:133429/139‑1 (MQ=255)
                                                 cTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAAt                                                                                <  1:621893/139‑1 (MQ=255)
                                                     tCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCAc                                                                            >  2:318273/1‑139 (MQ=255)
                                                          ggTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTaaa                                                                       <  1:718688/139‑1 (MQ=255)
                                                            tttCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAAt                                                                     <  2:127005/139‑1 (MQ=255)
                                                             ttCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATc                                                                    <  1:711930/139‑1 (MQ=255)
                                                                    ggTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACaa                                                             <  2:599573/139‑1 (MQ=255)
                                                                          aCGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCtt                                                       >  1:586554/1‑139 (MQ=255)
                                                                                ttAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTgac                                                                                                            >  2:563491/1‑80 (MQ=255)
                                                                                ttAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTgac                                                                                                            <  1:563491/80‑1 (MQ=255)
                                                                                     gCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATg                                            <  1:344661/139‑1 (MQ=255)
                                                                                         acaaTGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCca                                        >  2:594393/1‑139 (MQ=255)
                                                                                             tGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACaat                                    >  1:622738/1‑138 (MQ=255)
                                                                                               tCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCa                                                                   <  1:613582/106‑1 (MQ=255)
                                                                                               tCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCa                                                                   >  2:613582/1‑106 (MQ=255)
                                                                                                       cGCTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAAt                          >  1:44555/1‑139 (MQ=255)
                                                                                                         cTTGCGCCAGCCCCAACGCCATCCCCTGATccagtcc                                                                                                                              >  1:54639/1‑37 (MQ=255)
                                                                                                         cTTGCGCCAGCCCCAACGCCATCCCCTGATccagtcc                                                                                                                              <  2:54639/37‑1 (MQ=255)
                                                                                                         cTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCt                        >  1:740877/1‑139 (MQ=255)
                                                                                                         cTTGCGCCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCt                        >  1:380968/1‑139 (MQ=255)
                                                                                                             cgcCAGCCCCAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTaa                    >  2:624747/1‑139 (MQ=255)
                                                                                                                   ccccAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCg                                                                                         <  1:535342/64‑1 (MQ=255)
                                                                                                                   ccccAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCg                                                                                         >  2:535342/1‑64 (MQ=255)
                                                                                                                   ccccAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCca              >  1:536563/1‑139 (MQ=255)
                                                                                                                    cccAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCcac             <  2:380968/139‑1 (MQ=255)
                                                                                                                    cccAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCcac             <  1:274892/139‑1 (MQ=255)
                                                                                                                     ccAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCcaca            <  1:555300/139‑1 (MQ=255)
                                                                                                                     ccAACGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCcaca            <  2:740877/139‑1 (MQ=255)
                                                                                                                       aaCGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATg          <  1:594393/139‑1 (MQ=255)
                                                                                                                         cGCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGt        <  1:624747/139‑1 (MQ=255)
                                                                                                                          gCCATCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTc       <  2:536563/139‑1 (MQ=255)
                                                                                                                              tCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAAt                                                                     >  1:767981/1‑73 (MQ=255)
                                                                                                                              tCCCCTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAAt                                                                     <  2:767981/73‑1 (MQ=255)
                                                                                                                               ccccTGATCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAATTAGCGCAAATCTTTAACGGCCACATGGTCCatat  <  1:318273/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
TTCGCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACAGCCCGCTTGCGCCAGCCCCAACGCCATCCCCTGACCCAGTCCAGTATCACAACCAGTGACGACCGCAACTTTACCTTCGAGAGAAAATGCACTTAAAATCATAACAATACCTTATTCGTTTATGCCCACAACTAGCGCAAATCTTTAACGGCCACATGGTCCATAT  >  NC_000913/2983057‑2983322

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