Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,602,723 G→A R302R (CGG→CGA)  hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,7230GA97.5% 144.9 / NA 40R302R (CGG→CGA) hyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/18);  minor base C (1/0);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTATTTTCGGCATCCTGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGGCTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTC  >  NC_000913/2602578‑2602870
                                                                                                                                                 |                                                                                                                                                   
gTATTTTCGGCATCCTGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGactgct                                                                                                                                                >  1:646485/1‑151 (MQ=255)
    tttCGGCATCCTGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTgg                                                                                                                                              >  2:841730/1‑149 (MQ=255)
     ttCGGCATCCTGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCt                                                                                                                                           >  2:362597/1‑151 (MQ=255)
       cGGCATCCTGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTgg                                                                                                                                         <  1:304461/151‑1 (MQ=255)
       cGGCATCCTGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTgg                                                                                                                                         >  1:85389/1‑151 (MQ=255)
               tGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGt                                                                                                                                 >  1:291375/1‑151 (MQ=255)
                  aaGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCg                                                                                                                               >  2:174720/1‑150 (MQ=255)
                     tAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGaaa                                                                                                                            >  1:175916/1‑150 (MQ=255)
                        cGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGt                                                                                                                        >  2:281379/1‑151 (MQ=255)
                            ggATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCgg                                                                                                                     <  1:614556/150‑1 (MQ=255)
                               tCTGCTGGCGCAAACGGGTTTGCCTCTGTGGCGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCatt                                                                                                                 >  2:57551/1‑151 (MQ=255)
                                  gctgGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCattatt                                                                                                              >  2:139870/1‑151 (MQ=255)
                                  gctgGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCa                                                                                                                   <  2:216624/146‑1 (MQ=255)
                                    tgGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTAtttt                                                                                                            >  2:357867/1‑151 (MQ=255)
                                     ggCGCAAACGGGTTTTCCTCTGTGGTGGGGCATTCTGGTGAGGGCGATCGGCGCCATCTCCGCGCTTATGGGCGTGCTATATGCGCAGGCATAACAGGATATTAAACAACTGCTGGCCTGGAGTACCGTCGAAAACGACGGCATTAttt                                                                                                             <  1:551751/149‑1 (MQ=255)
                                      gcgcAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGc                                                                                                          >  1:639794/1‑151 (MQ=255)
                                       cgcAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCt                                                                                                         <  2:856281/151‑1 (MQ=255)
                                        gcAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTg                                                                                                        >  1:127320/1‑151 (MQ=255)
                                          aaaCGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTg                                                                                                        <  2:150104/149‑1 (MQ=255)
                                           aaCGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAACCGCCTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCt                                                                                                         >  1:569330/1‑147 (MQ=255)
                                           aaCGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCa                                                                                                     <  2:744661/151‑1 (MQ=255)
                                            aCGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGc                                                                                                      >  1:510733/1‑149 (MQ=255)
                                                  ttGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGgt                                                                                               <  2:124651/150‑1 (MQ=255)
                                                    gCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGtgtg                                                                                            <  2:752687/151‑1 (MQ=255)
                                                     ccTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTgg                                                                                           >  2:456590/1‑151 (MQ=255)
                                                          gtggtgGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATg                                                                                      <  2:340000/151‑1 (MQ=255)
                                                                      tCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTc                                                                           <  2:639836/150‑1 (MQ=255)
                                                                          gTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGc                                                                      <  2:38658/151‑1 (MQ=255)
                                                                            gATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCa                                                                     <  1:840807/150‑1 (MQ=255)
                                                                               ggCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGAc                                                                 >  2:659958/1‑151 (MQ=255)
                                                                                     cGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCgctg                                                           <  2:256549/151‑1 (MQ=255)
                                                                                      ggCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCgctgc                                                          >  2:418298/1‑151 (MQ=255)
                                                                                                       ccTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACt                                          <  1:656444/150‑1 (MQ=255)
                                                                                                         tGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTg                                         >  1:209665/1‑149 (MQ=255)
                                                                                                             cGTGCTATATGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGgc                                   >  2:726655/1‑151 (MQ=255)
                                                                                                                      tGCGCTGGCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGATCACCGTGGTTGGACTGCTCGGCGCACTGtt                           >  2:639683/1‑150 (MQ=255)
                                                                                                                            ggCGGAACAGGATATCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCt                     <  1:496457/150‑1 (MQ=255)
                                                                                                                                      gATATCAAACGACTGCTGGCCTTGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACGGCTCGGCGCACTGTTTCATCTGCTCAACCATg          <  1:400608/151‑1 (MQ=255)
                                                                                                                                          tCAAACGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCt      >  1:430185/1‑151 (MQ=255)
                                                                                                                                            aaaCGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGt    >  2:354959/1‑151 (MQ=255)
                                                                                                                                              aCGACTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTc  <  2:405078/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTATTTTCGGCATCCTGAAAGTAGCGATGGATCTGCTGGCGCAAACGGGTTTGCCTCTGTGGTGGGGCATTCTGGTGATGGCGATCGGCGCAATCTCCGCGCTCCTGGGCGTGCTATATGCGCTGGCGGAACAGGATATCAAACGGCTGCTGGCCTGGAGTACCGTCGAAAACGTCGGCATTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTC  >  NC_000913/2602578‑2602870

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