Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,851,603 C→T 7.4% S201S (AGC→AGT hypB → GTP hydrolase involved in nickel liganding into hydrogenases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,851,6030CT 7.4% 147.8 / 3.2 56S201S (AGC→AGThypBGTP hydrolase involved in nickel liganding into hydrogenases
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base T (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

CGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAAC  >  NC_000913/2851459‑2851741
                                                                                                                                                |                                                                                                                                          
cGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCg                                                                                                                                        <  2:661579/149‑1 (MQ=255)
  gCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGAt                                                                                                                                      >  1:924018/1‑149 (MQ=255)
   cacaCCAGCGATGCCGGTGACCACCGGTAAAGGCTGCCATCTTGACTCAAAGATGATTGCCGACGACGCACCGCGTCTGCCACTGGTCGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGAtc                                                                                                                                     <  2:323216/149‑1 (MQ=255)
      acCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGCCGATATCGGTATTATGTTTATCGAACACGTTGGCAACCTCGTATCCCCGGCCAGCTTCGATCTCg                                                                                                                                  >  1:568089/1‑149 (MQ=255)
        cAGCGCTTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGt                                                                                                                                <  1:397298/149‑1 (MQ=255)
        cAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGCCGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACGGGACGATAACGGTATTCTGTTTATCTAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCGCGGt                                                                                                                                >  2:786950/1‑149 (MQ=255)
        cAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGt                                                                                                                                >  1:111520/1‑149 (MQ=255)
        cAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGt                                                                                                                                >  1:115295/1‑149 (MQ=255)
         aGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTg                                                                                                                               >  2:242778/1‑149 (MQ=255)
          gCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTCTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCATCTTCGATATCCGTGa                                                                                                                              >  2:740740/1‑149 (MQ=255)
           cGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGaa                                                                                                                             <  2:598036/149‑1 (MQ=255)
               tCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAac                                                                                                                         <  1:176069/149‑1 (MQ=255)
                   gTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAACCGTTGGCAACCTCGTATGCCCGGCCAGCTTCGGACTAGGTGAAAAACACaa                                                                                                                     >  2:578250/1‑149 (MQ=255)
                   gTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACaa                                                                                                                     <  1:560773/149‑1 (MQ=255)
                            ggTAAAGGCTCCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAAACTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGt                                                                                                            >  2:328335/1‑149 (MQ=255)
                             gTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTg                                                                                                           >  1:24216/1‑149 (MQ=255)
                              tAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGc                                                                                                          >  2:595134/1‑149 (MQ=255)
                                aaGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGACAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCtt                                                                                                        >  2:17494/1‑149 (MQ=255)
                                  ggCTGCCATCTTTACGCACATATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTAGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTc                                                                                                      <  2:707376/149‑1 (MQ=255)
                                  ggCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTc                                                                                                      <  2:709029/149‑1 (MQ=255)
                                     tGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGt                                                                                                   >  1:497817/1‑149 (MQ=255)
                                        cATCTTGACGCACAGATGATTGCCGACGCCGCACCTCGTCTGCCACTGGACGATAACGGTATTATGTTTATCGAAAACGTTGGCAACCTCGTATGCACGGCCGGCGTCGATCTCGCGGAAAAACCCGAAGTGGCGGTGCTTTCCGTTAc                                                                                                >  2:204245/1‑149 (MQ=255)
                                             tGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAg                                                                                           <  1:976907/149‑1 (MQ=255)
                                             tGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAg                                                                                           <  1:764283/149‑1 (MQ=255)
                                                cGCACAGAGGATAGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTg                                                                                        <  1:184112/149‑1 (MQ=255)
                                                  cacaGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGaa                                                                                      <  1:28592/149‑1 (MQ=255)
                                                           ttGCCTACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGAATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCa                                                                             <  1:928170/149‑1 (MQ=255)
                                                             gCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACt                                                                           <  2:111520/149‑1 (MQ=255)
                                                             gCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACt                                                                           <  2:115295/149‑1 (MQ=255)
                                                              ccGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCGAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAACCCACTg                                                                          >  2:849824/1‑149 (MQ=255)
                                                              ccGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGTTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTg                                                                          >  1:591998/1‑149 (MQ=255)
                                                              ccGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGTTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTg                                                                          >  1:591996/1‑149 (MQ=255)
                                                               cGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGa                                                                         <  2:367920/149‑1 (MQ=255)
                                                               cGACGCCGCACCCCGGCTGCCACTGGACGATTACGGTAATCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCTATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGAATAACCACTGa                                                                         <  2:196661/149‑1 (MQ=255)
                                                                    ccgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAatat                                                                    <  2:734477/149‑1 (MQ=255)
                                                                            cgTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCatat                                                            <  2:330260/149‑1 (MQ=255)
                                                                                 gCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTTGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACAAAAAGTGGCGGTGCTTTCCGTTACCTAAGGTGAAGACAAAACACTGAAATATCCTCATATGttta                                                       >  1:274100/1‑148 (MQ=255)
                                                                                 gCCACTGGACGATAACGGTATTCTCTTTATCGAAAACGATGGCAACCTCGTATGCACCGCCAGCTTCGATCTCGCTGAAAAACACAACGTGGCGGTGCTTTCCGTTTCCTAAGGTGAAGAAAAACCACTGAAATATCCGCATATGtttt                                                       >  2:891326/1‑148 (MQ=255)
                                                                                     cTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGCCAAGCCACTGAAATCTCCGCATATGTTTgccgc                                                   >  2:768819/1‑149 (MQ=255)
                                                                                             tAAAGGTATTCTGTTTAACGAAAACGTTGGCAACCTCTTATGCCCGGACATCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGc                                           <  1:578250/149‑1 (MQ=255)
                                                                                                 ggTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGAt                                       >  2:51285/1‑149 (MQ=255)
                                                                                                  gTATTCTGTTAATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAACCACAACGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGAtg                                      <  2:8718/149‑1 (MQ=255)
                                                                                                     ttCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGGTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCGGCCGACTCGctgctgctg                                   >  2:791140/1‑149 (MQ=255)
                                                                                                               cGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTGTCCGTTACCGAAGGGGAAGACCAACCACTGAAATATCCGCATATGTTTGCCGCCGCCGCGCTGATGCTGCTCACCAAAg                         >  1:865027/1‑149 (MQ=255)
                                                                                                                gAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAGACAAACCACTGAAAGATCCGCATATGTTTGCCGCCGCCGCGCTGATGCTGCTCAACAAAGt                        >  2:599358/1‑149 (MQ=255)
                                                                                                                 aaaaCTTTGGGAACAGCGTATTACCGGCCAGCTTCGAGCTCGGTGAAAAACACAACGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGtt                       <  1:943943/149‑1 (MQ=255)
                                                                                                                 aaaaCGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGtt                       >  2:423803/1‑149 (MQ=255)
                                                                                                                 aaaaCGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGtt                       <  1:511517/149‑1 (MQ=255)
                                                                                                                 aaaaCGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGtt                       >  1:476411/1‑149 (MQ=255)
                                                                                                                         ggCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGtt               <  2:24216/149‑1 (MQ=255)
                                                                                                                          gCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTg              >  1:1016242/1‑149 (MQ=255)
                                                                                                                              ccgcGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCGTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGt          <  2:924018/146‑1 (MQ=255)
                                                                                                                              ccTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAATACACAAAGTGGAGGTGCTTTCCGTTACAGAAGGTGAAGACAAACCACTGCTATATCCGCATATGTTTGCCGCCGCGTCGCTGATGCTGCTCAACAAAGTTGACATGGTGCCGt          >  2:87544/1‑149 (MQ=255)
                                                                                                                                 cGTATGCCCGGCCAGTTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTAtc       <  2:591998/149‑1 (MQ=255)
                                                                                                                                 cGTATGCCCGGCCAGTTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTAtc       <  2:591996/149‑1 (MQ=255)
                                                                                                                                      gCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAAc  >  2:652470/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                          
CGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAAC  >  NC_000913/2851459‑2851741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: