Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,073,554 A→C F6L (TTT→TTG)  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,5540AC100.0% 176.6 / NA 50F6L (TTT→TTG) yihU3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (24/26);  total (24/26)

CGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGAAAT  >  NC_000913/4073409‑4073701
                                                                                                                                                 |                                                                                                                                                   
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)
                                                                                             ctggctgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAAt                                                   >  1:219456/4‑151 (MQ=255)
                                                                                             ctggctgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAAt                                                   >  2:214247/4‑151 (MQ=255)
                                                                                                   gcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATc                                              <  1:13088/150‑1 (MQ=255)
                                                                                                      aTAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTc                                           >  1:666625/1‑150 (MQ=255)
                                                                                                         aaTTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaaa                                       <  2:60110/151‑1 (MQ=255)
                                                                                                         aaTTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaa                                        >  2:698133/1‑150 (MQ=255)
                                                                                                            tGCTCGTCATTGGCGAACCCATTTGTCCTAAACCGATCAACTCGATTGCTGCCATAACCCTCTCCTGAATTCAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaaaaa                                     <  1:306638/150‑1 (MQ=255)
                                                                                                              cTCGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTg                                  >  1:714448/1‑151 (MQ=255)
                                                                                                                cGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTg                                  <  2:323870/149‑1 (MQ=255)
                                                                                                                cGCCATTGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGAc                                >  1:214986/1‑151 (MQ=255)
                                                                                                                     ttGGCGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCg                           >  1:538803/1‑151 (MQ=255)
                                                                                                                     ttGGCGAACCCAGTTGTCCTAAGCCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCTCTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGa                                 <  2:173767/145‑1 (MQ=255)
                                                                                                                         cGAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCAc                       >  1:245842/1‑151 (MQ=255)
                                                                                                                          gAACCCATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACt                      >  1:538444/1‑151 (MQ=255)
                                                                                                                              ccATTTGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACttttt                  >  1:485155/1‑151 (MQ=255)
                                                                                                                                  ttGTCCTAAACCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTaaa               <  1:458358/150‑1 (MQ=255)
                                                                                                                                          aaCCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTg      <  2:498674/151‑1 (MQ=255)
                                                                                                                                          aaCCGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTg      <  2:147217/151‑1 (MQ=255)
                                                                                                                                             cGATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTGACTTTTTAAACAATTGGTGaaa   <  1:530389/151‑1 (MQ=255)
                                                                                                                                              gATCAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGAAAt  >  2:167771/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGCCTGCGCCGGGTTGGCGGCGGGAGTCGCACCTTTGTCTACCAGATGCCGCACAGCTTCGGCGTTCACATCAAAGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGCCGTCACTTTTTAAACAATTGGTGAAAT  >  NC_000913/4073409‑4073701

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