Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,073,504 T→G Q23P (CAA→CCA)  yihU ← 3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,073,5040TG100.0% 199.0 / NA 55Q23P (CAA→CCA) yihU3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (23/32);  total (23/32)

CGCACCAGATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGT  >  NC_000913/4073357‑4073649
                                                                                                                                                   |                                                                                                                                                 
cgcACAAGATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTggct                                                                                                                                                <  2:341147/151‑1 (MQ=255)
   aCCAGATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCa                                                                                                                                             <  1:188890/151‑1 (MQ=255)
   aCCAGATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCa                                                                                                                                             <  2:70950/151‑1 (MQ=255)
    ccAGATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGATGATGCCCTGGCTGCa                                                                                                                                             >  1:624080/1‑150 (MQ=255)
     cAGATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAAt                                                                                                                                           <  1:221104/151‑1 (MQ=255)
       gATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATaa                                                                                                                                         <  1:280276/151‑1 (MQ=255)
        aTCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATaaa                                                                                                                                        <  2:365404/151‑1 (MQ=255)
              aTTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTc                                                                                                                                  <  2:722287/151‑1 (MQ=255)
               ttCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTc                                                                                                                                  <  1:565153/150‑1 (MQ=255)
                tCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCg                                                                                                                                 <  2:346625/150‑1 (MQ=255)
                tCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGc                                                                                                                                >  2:431304/1‑151 (MQ=255)
                  ggCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCGCCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCa                                                                                                                              <  2:433566/151‑1 (MQ=255)
                  ggCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCa                                                                                                                              <  1:318258/151‑1 (MQ=255)
                     aaCATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTg                                                                                                                           <  2:828461/151‑1 (MQ=255)
                     aaCATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTg                                                                                                                           <  1:325138/151‑1 (MQ=255)
                         tGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCg                                                                                                                        >  2:911833/1‑150 (MQ=255)
                            tAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAAcc                                                                                                                    >  1:207164/1‑151 (MQ=255)
                              aTGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAAccc                                                                                                                   >  2:242828/1‑150 (MQ=255)
                                gATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCAtt                                                                                                                <  1:271825/151‑1 (MQ=255)
                                   aaaTTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCAttt                                                                                                               >  2:34325/1‑149 (MQ=255)
                                     aTTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTcc                                                                                                           <  2:422815/151‑1 (MQ=255)
                                        cGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTa                                                                                                         <  1:547608/150‑1 (MQ=255)
                                             tCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCg                                                                                                    >  1:861851/1‑150 (MQ=255)
                                               tttAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGAtc                                                                                                 >  1:498658/1‑150 (MQ=255)
                                                ttAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGAtc                                                                                                 >  2:209281/1‑149 (MQ=255)
                                                 tAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGAtca                                                                                                >  2:176995/1‑148 (MQ=255)
                                                    cGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAAcgc                                                                                            >  2:493265/1‑151 (MQ=255)
                                                          gcgcCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATtgc                                                                                      <  2:196302/151‑1 (MQ=255)
                                                             ccGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATtgctgc                                                                                   >  2:868049/1‑151 (MQ=255)
                                                                 gTTGGCGGCGGGAGGCGCACCTTTGGATAACAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCCAAGCGATTGCTGCCa                                                                                 >  1:129958/1‑149 (MQ=255)
                                                                   tGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCATAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAAc                                                                             >  2:750272/1‑151 (MQ=255)
                                                                   tGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATaa                                                                              <  1:394414/150‑1 (MQ=255)
                                                                       ggcggGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCctc                                                                         >  1:594021/1‑151 (MQ=255)
                                                                        gcggGAGTCGCACCTTTGTCTACCAGATGCCGCACCGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCctc                                                                         >  1:579236/1‑150 (MQ=255)
                                                                              gTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGaa                                                                  <  1:85958/151‑1 (MQ=255)
                                                                               tCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAAt                                                                 >  2:255908/1‑151 (MQ=255)
                                                                                cGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAAt                                                                 >  2:536692/1‑150 (MQ=255)
                                                                                                aGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACtttt                                                <  1:835128/151‑1 (MQ=255)
                                                                                                  aTGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTg                                              <  1:783063/151‑1 (MQ=255)
                                                                                                      cGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCAt                                          >  2:729072/1‑151 (MQ=255)
                                                                                                           aGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATg                                     <  1:614331/151‑1 (MQ=255)
                                                                                                              ttCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGAc                                   >  2:924178/1‑150 (MQ=255)
                                                                                                              ttCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGAc                                   >  1:427215/1‑150 (MQ=255)
                                                                                                               tCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACa                                  <  2:152064/150‑1 (MQ=255)
                                                                                                                  gCGTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgct                              <  1:66969/151‑1 (MQ=255)
                                                                                                                    gTTCACATCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgcttt                            <  1:541998/151‑1 (MQ=255)
                                                                                                                       cacaTCAAAGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctttgc                          <  1:535658/150‑1 (MQ=255)
                                                                                                                               aGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGtt                 <  2:715341/151‑1 (MQ=255)
                                                                                                                               aGACGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGtt                 <  1:568717/151‑1 (MQ=255)
                                                                                                                                 aCGCGAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGtttt               >  1:918533/1‑151 (MQ=255)
                                                                                                                                   gcgAAGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTg             <  1:911812/151‑1 (MQ=255)
                                                                                                                                      aaGTTGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTg             <  1:682059/148‑1 (MQ=255)
                                                                                                                                          tGATGCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTa       <  1:652028/150‑1 (MQ=255)
                                                                                                                                              gCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTg   <  2:871196/150‑1 (MQ=255)
                                                                                                                                              gCCCTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGt  <  2:624100/151‑1 (MQ=255)
                                                                                                                                               cccTGGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGt  >  1:314004/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGCACCAGATCGCCATTCGGCAACATGGTAATGATAAATTCGGCATCTTTAGCGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGT  >  NC_000913/4073357‑4073649

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