Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,602,375 +G coding (558/2019 nt) hyfB → hydrogenase 4 component B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,3711.G97.4% 132.5 / ‑3.5 39coding (554/2019 nt)hyfBhydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (22/16);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAGCCGGGATGATGCCGTTGCACAG  >  NC_000913/2602232‑2602503
                                                                                                                                            |                                                                                                                                     
gTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGAt                                                                                                                                        <  2:748318/140‑1 (MQ=255)
                tgTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCgg                                                                                                                       >  2:542445/1‑140 (MQ=255)
                tgTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCgg                                                                                                                       <  1:49866/140‑1 (MQ=255)
                tgTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATg                                                                                                                                      >  1:712535/1‑125 (MQ=255)
                tgTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATg                                                                                                                                      <  2:712535/125‑1 (MQ=255)
                     cTGTCTTCCTGGTTTGTGGTCATTGTCAGCCAGGATACAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGcc                                                                                                                  <  1:420331/140‑1 (MQ=255)
                      tGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCt                                                                                                                 >  2:183156/1‑140 (MQ=255)
                           tCCTGGTTTCTGGTCATGGCCAGGCAGGATACAACTTCGGTCAACGCTCGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGAtt                                                                                                            <  2:342943/140‑1 (MQ=255)
                             cTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGAtttt                                                                                                          >  1:467241/1‑140 (MQ=255)
                                       gTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTcc                                                                                                >  1:5472/1‑140 (MQ=255)
                                        tCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCg                                                                                               <  2:5472/140‑1 (MQ=255)
                                        tCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCg                                                                                               <  2:467241/140‑1 (MQ=255)
                                         cATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGc                                                                                              <  1:542445/140‑1 (MQ=255)
                                            tGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCAc                                                                                            <  2:612541/139‑1 (MQ=255)
                                            tGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCAc                                                                                            >  1:612541/1‑139 (MQ=255)
                                              ccAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCt                                                                                         <  1:183156/140‑1 (MQ=255)
                                                 ggcaggATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTc                                                                                      >  2:490958/1‑140 (MQ=255)
                                                     ggATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCt                                                                                         <  1:92656/133‑1 (MQ=255)
                                                     ggATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCt                                                                                         >  2:92656/1‑133 (MQ=255)
                                                      gATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACttt                                                                                 >  2:19170/1‑140 (MQ=255)
                                                      gATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACttt                                                                                 >  1:508462/1‑140 (MQ=255)
                                                         aaaaCGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTtctc                                                                              <  1:490958/140‑1 (MQ=255)
                                                                cGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCt                                                                       >  2:819414/1‑140 (MQ=255)
                                                                   tCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGc                                                                    >  1:639138/1‑140 (MQ=255)
                                                                            gCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGgt                                                           >  1:811305/1‑140 (MQ=255)
                                                                             cATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCTCACGCTTGCACATTGTCCGGGACTGGCGTCGGCGgtg                                                          >  1:258569/1‑140 (MQ=255)
                                                                                 ctctACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTcc                                                      >  1:63711/1‑140 (MQ=255)
                                                                                      cTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCctgctg                                                 >  1:57730/1‑140 (MQ=255)
                                                                                         ttttATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGcc                                              <  2:63711/140‑1 (MQ=255)
                                                                                                           aTCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGcgcg                            <  2:508462/140‑1 (MQ=255)
                                                                                                           aTCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGcgcg                            <  2:57730/140‑1 (MQ=255)
                                                                                                             cGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGaa                          <  1:19170/140‑1 (MQ=255)
                                                                                                                tGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGTTGTTCCTGCTGGCCGTTTTCGATTTTGGCGCGAAAGc                       >  1:263997/1‑140 (MQ=255)
                                                                                                                tGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAg                        >  1:69304/1‑139 (MQ=255)
                                                                                                                tGCTGATAATGATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAg                        <  2:69304/139‑1 (MQ=255)
                                                                                                                        aTGATCGCCTTCTTGCTCATGGCGGCGCGAACGCGGCAGCCTCGATTTCGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAGCCGGgatga               <  2:263997/140‑1 (MQ=255)
                                                                                                                          gATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAGCCGGgatga               <  2:133185/138‑1 (MQ=255)
                                                                                                                          gATCGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAGCCGGgatga               >  1:133185/1‑138 (MQ=255)
                                                                                                                             cGCCTTCTTGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAGCCGGGATGATGCCg          >  2:866225/1‑140 (MQ=255)
                                                                                                                             cGCCTTCTTGCTGATCGGGGCGTGACAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAGCCGGGATGATGCCg          <  1:411230/140‑1 (MQ=255)
                                                                                                                                cttcttGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTACTGCTGGCCTTTTTCGTGTTTGGCGGGAAAACCGGGATGTTGCCGttt       >  2:411230/1‑139 (MQ=255)
                                                                                                                                   cttGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTATTCGGTTTTGGCGCGATAGCCACGATGATGGCGTTGcac    >  2:404742/1‑140 (MQ=255)
                                                                                                                                     tGCTGATGGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCTGTTTTGGGGCGAAAGCCGGGATGATGCCGTTGCACAg  >  1:292720/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                     
GTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGATCAACGCTGGCATGCTCTACTTTTTTATCGCCCACGCCGGATCGGTGCTGATAATGATCGCCTTCTTGCTGATGGGGCGCGAAAGCGGCAGCCTCGATTTTGCCAGTTTCCGCACGCTTTCACTTTCTCCGGGGCTGGCGTCGGCGGTGTTCCTGCTGGCCTTTTTCGGTTTTGGCGCGAAAGCCGGGATGATGCCGTTGCACAG  >  NC_000913/2602232‑2602503

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