Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,843,708 T→A 7.6% H168L (CAC→CTC)  hycG ← formate hydrogenlyase subunit HycG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,843,7080TA 7.6% 134.5 / 4.0 54H168L (CAC→CTC) hycGformate hydrogenlyase subunit HycG
Reads supporting (aligned to +/- strand):  ref base T (20/29);  new base A (2/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

AATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGAA  >  NC_000913/2843568‑2843847
                                                                                                                                            |                                                                                                                                           
aaTCGTCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGa                                                                                                                                           <  2:547989/143‑1 (MQ=255)
   catcGGCTATCGCACGACCGCAACGACACCCCGCCAGGCGACGTGCTGCGCGATCCACTTTCACGCGCAGCGGCTGCACCATAACCCCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAAtt                                                                                                                                        <  2:653578/143‑1 (MQ=255)
   catcGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAAtt                                                                                                                                        >  1:252309/1‑143 (MQ=255)
    atcGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCACGCAGGCTCTCCGCCGGCTGTTCATCCAGTCGACCCGGCCCGCGGGCGTGaacac                                                                                                                                       >  2:894026/1‑140 (MQ=255)
      cGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTc                                                                                                                                     <  1:318920/143‑1 (MQ=255)
         cAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGCGAATTTTCTGc                                                                                                                                  >  2:426468/1‑143 (MQ=255)
           aTCTGACGCTCGTCACGAGAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTc                                                                                                                                <  2:135931/143‑1 (MQ=255)
              tGACGACCGTAACGATAACCCGCCAGGCGCCGTGCTTCGCGATCCACTTTCACGCGCAGCAGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGTATTTTCTGCTCcag                                                                                                                             >  1:191204/1‑143 (MQ=255)
              tGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACCGGCGTGAATTTTCTGCCCcag                                                                                                                             >  2:662239/1‑143 (MQ=255)
                  gacCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCCgg                                                                                                                         >  1:810722/1‑143 (MQ=255)
                       gaccGATAACCCGCCAGGCGACGTGCTTCGCGATCCGCTTGCACGCGCAGCGGCGGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATGTTCTGCTCCAGCAGGCCGAg                                                                                                                    <  1:209561/140‑1 (MQ=255)
                        aaCGATAACCCGCCAGGCGACGTGCTTCGCGCTCCACGTCCACGCGGTGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAgc                                                                                                                   <  2:822754/143‑1 (MQ=255)
                        aaCGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCgcgc                                                                                                                   >  1:489661/1‑143 (MQ=255)
                           gATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCCCCGGTTGTTCATCCAGTTGCCACCGCCCACGCGCGTGAACTTTCTGCTCCAGCAGGCCggccacc                                                                                                                >  2:185307/1‑137 (MQ=255)
                                 ccGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCt                                                                                                                                    <  2:433165/117‑1 (MQ=255)
                                 ccGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCt                                                                                                                                    >  1:433165/1‑117 (MQ=255)
                                 ccGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCAACATATCAACATGCAGGATCTCCCCAGCCCATTCATCCGGTACACCAGGCCCACGGGCGTGAATTCTCTGCTCCAACAGGCCGGGCGccgtcggc                                                                                                          >  2:670235/1‑137 (MQ=255)
                                   gCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCaaa                                                                                                        <  1:760221/143‑1 (MQ=255)
                                   cccAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCaaa                                                                                                        <  1:760283/142‑1 (MQ=255)
                                      aGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGTAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGAGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGcc                                                                                                     >  1:324652/1‑143 (MQ=255)
                                      aGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGTAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGAGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGcc                                                                                                     >  1:834923/1‑143 (MQ=255)
                                       ggCGACGTGCTTCGCGATCCACTTTCACGCGCAACATGCGCTGCATATCGCCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCg                                                                                                    >  1:378798/1‑143 (MQ=255)
                                              ggCTGCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCg                                                                                             <  1:717736/142‑1 (MQ=255)
                                                     cgATCCGCCTTCCCGCGCAGCGGCTGCACCATACCACCATGCAGGATCTCCGCCTGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTggcggc                                                                                      <  2:582881/143‑1 (MQ=255)
                                                      gATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCa                                                                                     <  1:717547/143‑1 (MQ=255)
                                                        tCCACTTTCACGCGCAGCGGCTGCACCAGAGCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCACACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAgg                                                                                   <  1:89732/143‑1 (MQ=255)
                                                           aCTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATCTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTCCAGCGTGGCGACacccgc                                                                                >  2:613065/1‑138 (MQ=255)
                                                                 aCGCGCCGCGGCTGCACCATCTCCCCTTGCGGGATCTCCGCCGGTTGTTCACCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCgg                                                                          <  2:711524/143‑1 (MQ=255)
                                                                  cgcgCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCggg                                                                         <  1:718739/143‑1 (MQ=255)
                                                                  cgcgCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCAATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCggg                                                                         >  1:77600/1‑143 (MQ=255)
                                                                   gcgcAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGGACAGCGGGGCGGCAGGCGTTGGCGGGc                                                                        >  1:785427/1‑143 (MQ=255)
                                                                       aGCGTCTCCACCATACCACCATCCAGGATCGCCGCCGGTTGTTCATCCGGTTCACCCGGCGCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGATCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGcc                                                                    <  2:640936/143‑1 (MQ=255)
                                                                            cTGCACGATAGCACCATGCAGGATCTCCGCCGGTTGTTCCTCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGa                                                               <  2:773055/143‑1 (MQ=255)
                                                                              gCACCATATCACCATGCAGGATCTCCGCCTGTTGTTCATCCAGTTTACCCGGCTCACGGGCGTGAATTTTCTGCGCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGata                                                             <  1:17651/143‑1 (MQ=255)
                                                                                  cATCTCACCCTGCGGGATCGCCGCCGGTCGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATaaa                                                         <  2:896334/143‑1 (MQ=255)
                                                                                    tatCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACa                                                       >  1:570894/1‑143 (MQ=255)
                                                                                     aacACCATGCAGGATCTCACCCGGTGTGTCATCCAAGTCACCCGGCCCACGTGCGTGAATCTTCTGCTCCAGCAGGCCGAGCGCCAGTGGAAAGCTATACATCATCTCTGCAGGCGTTGGCGGGCAGCCAGGGATATAAACAt                                                      <  2:191204/141‑1 (MQ=255)
                                                                                                 gATCCCCCCCGGCTCTCCCTCCAGTCCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCTCCTAAGACTTCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACaa                                          <  1:126529/143‑1 (MQ=255)
                                                                                                            gttgttCATCCATTTCACCCGGCCCACGGGCGTGAATTTTGTGCTCCAGCAGGCCGAGTGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTa                               <  1:601823/143‑1 (MQ=255)
                                                                                                             ttgttCATCCAGTTCACCCGGCCCACGGGCGTGGATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTAc                              <  2:489661/143‑1 (MQ=255)
                                                                                                                  cATCCAGTTCACCCGGCCCACGGGCGAGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTCACGccc                         <  2:834923/143‑1 (MQ=255)
                                                                                                                  cATCCAGTTCACCCGGCCCACGGGCGAGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGccc                         <  2:324652/143‑1 (MQ=255)
                                                                                                                   aTCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGcc                                                                                                                <  1:574328/55‑1 (MQ=255)
                                                                                                                   aTCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGcc                                                                                                                >  2:574328/1‑55 (MQ=255)
                                                                                                                    tCCAGTTCACCCGGCCCACGGGCGTGAATTATCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACCACGGAGCGGCATGCCTGGGCGGCCCGCCAAGGACATGAACATCCACAGGCGCCACTTTACCCGTACCGccccc                       >  2:58802/1‑142 (MQ=255)
                                                                                                                     ccAGTTCACCCGGCCCACGGGCGTGAAAATTTTATTCCATCAGGCCGCGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGGTGGCGGGCAGCCAGGGATACAGACATCCACAGGGACAATGCTATCCGGACCGCCCcag                      >  2:542793/1‑142 (MQ=255)
                                                                                                                        gTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCCCAGGGACAAGTTTATCCGTCCCGCCCCACACg                   >  2:486521/1‑143 (MQ=255)
                                                                                                                             cccGG‑CCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTag             >  1:125013/1‑143 (MQ=255)
                                                                                                                                ggCCCACGGGCGTGAATTGGCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTagag           <  1:266482/143‑1 (MQ=255)
                                                                                                                                     aCGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTagagcacgt      >  2:372795/1‑138 (MQ=255)
                                                                                                                                     aCGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGt      <  1:426468/143‑1 (MQ=255)
                                                                                                                                       gggCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGTACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTgg    <  1:37049/143‑1 (MQ=255)
                                                                                                                                        ggCGTGTATTTTCTGCTCCTGCAGCCCGAGCGCCATTGCAAAGCCGTACTGCGTCTCGGCAGGCGTTCGCGGGCAGCCAGCGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACCCGCAGTAGAGACCGTTGa   <  2:773824/143‑1 (MQ=255)
                                                                                                                                         gCGTGACTTTTCTGCTCCAGCCGGCCGAGCGCCATCGCAAAGCCGAACGCCTTGGCGCCAGGCCTCGGCTGTCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGaa  <  1:243656/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
AATCATCGGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGAA  >  NC_000913/2843568‑2843847

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