Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 952,625 G→T 8.7% V310V (GTC→GTA)  pflB ← formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913952,6250GT 8.7% 132.5 / 3.1 55V310V (GTC→GTA) pflBformate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base T (3/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

TGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCC  >  NC_000913/952477‑952771
                                                                                                                                                    |                                                                                                                                                  
tgtTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATg                                                                                                                                                    <  1:801536/149‑1 (MQ=255)
      ggAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCGGATCGGGTGGGCAGAGAATCGTTCATGGTATTCCGGGGTACGCAGGGGGTGGAAGATACGCGGTTTCGTTACCggg                                                                                                                                              >  2:80795/1‑149 (MQ=255)
         aaCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtgg                                                                                                                                           <  1:59476/149‑1 (MQ=255)
           cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTc                                                                                                                                         <  1:721359/149‑1 (MQ=255)
           cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCCAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTc                                                                                                                                         <  1:140627/149‑1 (MQ=255)
            ggAAGCTGTTTTTGGTAACCAGGGTACGGCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGCGGTCa                                                                                                                                        >  1:427335/1‑149 (MQ=255)
                gCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa                                                                                                                                    >  1:759826/1‑149 (MQ=255)
                 cTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCACAAACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAt                                                                                                                                   >  1:792679/1‑149 (MQ=255)
                  tGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                                                                                  >  1:513554/1‑149 (MQ=255)
                         ggTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgcg                                                                                                                           <  1:87745/149‑1 (MQ=255)
                          gTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATCCGCAGTTTCATGACCAGGTGGTCAACCATTTCATgcag                                                                                                                          >  2:799685/1‑147 (MQ=255)
                                  ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                                                                                  >  1:411008/1‑149 (MQ=255)
                                  ggTACGACCGTCGAGGCCCATACCAACGATAGATTCGGTTGCCCAGATCGGGACGCCAGAGAAGAGTTCATCGTATTCCGGAGTACGCAGAAAAATTATCATTGTTTGATTCAGGACCAGGTGGTCAACCATTTCCTGCGCTACTTGtc                                                                                                                  >  1:350600/1‑148 (MQ=255)
                                     aCGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGACGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCgg                                                                                                               >  2:535113/1‑149 (MQ=255)
                                       gACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAGCAGTTCCTCGTATTCCGGAGTAAGCAGGAAGCAAAGCATACGCTGTTTCATGACCAGGTGGCTACCCTTTTCCTGCGCTTCTTGTTCGGTg                                                                                                             >  2:784327/1‑149 (MQ=255)
                                       gACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTACATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCCACCATTTCCTGCGCTTCTTGTTCGGTg                                                                                                             >  2:458851/1‑149 (MQ=255)
                                         ccGTCGAGGCCCATACCACCGATAGATTCGGTTGGCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGAt                                                                                                           <  1:792484/149‑1 (MQ=255)
                                                    cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCtt                                                                                                >  2:505697/1‑149 (MQ=255)
                                                     aggccaccGATAGATGCGGTTGCCCAGCTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                                                                                               <  2:302831/146‑1 (MQ=255)
                                                      tACCACCGATAGATTCGGTTGCCCAGGTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTTTTGTTCGGTGATCTTGCCAGCTTTc                                                                                              <  1:400402/149‑1 (MQ=255)
                                                        ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTcgg                                                                                            <  1:472933/149‑3 (MQ=255)
                                                         caccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATTACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAgg                                                                                           >  2:263637/1‑149 (MQ=255)
                                                         caccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATTACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAgg                                                                                           >  2:263432/1‑149 (MQ=255)
                                                            cGATAGATTCGGGTGCCCAGATCGCGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCa                                                                                        <  2:29252/149‑1 (MQ=255)
                                                                 gattgcgTTGCCCATGTTGGTGCGCCAGAGAGCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGTTGGTCAACCATTTCCTGCGCTTCTTGTTCGTTGATCTTGCCAGCTTTCAGGTCACGTTc                                                                                   <  2:356012/143‑1 (MQ=255)
                                                                      ggTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAACATCCGCAGTATCATGACCAGGGGGTCAACCATTTCCTGCGCTTCTTGTTCAGTGATCTTGCACTCTCTCACGTCACGTTCGATGt                                                                              >  2:18822/1‑149 (MQ=255)
                                                                       gTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTa                                                                             >  2:682836/1‑149 (MQ=255)
                                                                           cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCGGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTAcaca                                                                         >  1:144871/1‑149 (MQ=255)
                                                                            ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGCACGAAAACCAAACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGCCACGTTCGATGTacacct                                                                        >  2:330814/1‑147 (MQ=255)
                                                                              aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTGAACATcc                                                                      >  1:169076/1‑149 (MQ=255)
                                                                              aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTAAACATcc                                                                      >  2:351609/1‑149 (MQ=255)
                                                                               gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATTACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCa                                                                     <  1:263432/149‑1 (MQ=255)
                                                                               gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATTACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCa                                                                     <  1:263637/149‑1 (MQ=255)
                                                                                      tCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCGTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGt                                                              <  2:591384/149‑1 (MQ=255)
                                                                                        gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACTATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCACCTTTCAGGACACGTTCGATGTACACATCCAGGGAGGTgg                                                            >  1:555807/1‑149 (MQ=255)
                                                                                        gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTgg                                                            <  1:799685/149‑1 (MQ=255)
                                                                                                aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACAAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGa                                                    <  2:587472/149‑1 (MQ=255)
                                                                                                       aTCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAgg                                             >  2:429765/1‑149 (MQ=255)
                                                                                                       aTCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGGGGGGGTACGACCGagag                                             >  2:544872/1‑146 (MQ=255)
                                                                                                             ttccacaGTACGCAGGGAGCGTACCATACGGTGTTTCATGACCAGGTCGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTGGCCAGCTTTCAGGTCACGTTCGATGTACACACTCAGGAAGGTGGAGGTACGACCGAAGGACATTg                                       <  2:347943/143‑1 (MQ=255)
                                                                                                              tgcGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCTAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTgc                                      <  2:593849/147‑1 (MQ=255)
                                                                                                              tCCGGAGTACCCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGTTCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTgc                                      <  2:345859/149‑1 (MQ=255)
                                                                                                               ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAGGGACATTgca                                     >  1:703731/1‑149 (MQ=255)
                                                                                                                  gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCCGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTgcagca                                  <  2:759826/149‑1 (MQ=255)
                                                                                                                  gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTgcagca                                  >  1:700971/1‑149 (MQ=255)
                                                                                                                     tACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCg                               >  2:136138/1‑149 (MQ=255)
                                                                                                                          aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTg                          >  1:517289/1‑149 (MQ=255)
                                                                                                                             aaGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCGgaga                       >  2:582338/1‑149 (MQ=255)
                                                                                                                             aaGCGAACCATACGCAGTTTCATGACCAGGTAGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACCTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTgaga                       <  1:535113/149‑1 (MQ=255)
                                                                                                                                     cATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAAcag               >  2:181221/1‑149 (MQ=255)
                                                                                                                                         ggCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACtgcagc           <  2:336224/148‑1 (MQ=255)
                                                                                                                                          gCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCGTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGcc          <  1:299033/149‑1 (MQ=255)
                                                                                                                                           cAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCa         >  1:697717/1‑149 (MQ=255)
                                                                                                                                            aGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAg        <  1:429765/149‑1 (MQ=255)
                                                                                                                                            aGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAg        <  1:503732/149‑1 (MQ=255)
                                                                                                                                                  aTGCCCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGcc  <  2:513554/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCC  >  NC_000913/952477‑952771

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