Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,002,536 A→G 43.9% D63D (GAT→GAC hybG ← hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,002,5360AG43.9% 21.8 / 78.0 98D63D (GAT→GAChybGhydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base A (35/20);  new base G (32/11);  total (67/31)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATC  >  minE/2002485‑2002598
                                                   |                                                              
ggTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:2431461/69‑1 (MQ=255)
ggTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:181134/69‑1 (MQ=255)
ggTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:1478761/69‑1 (MQ=255)
ggTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:1453362/69‑1 (MQ=255)
ggTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:2143653/69‑1 (MQ=255)
ggTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:215782/69‑1 (MQ=255)
ggTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:2733878/69‑1 (MQ=255)
 gTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               <  1:1891519/68‑1 (MQ=255)
       tCGTAATCCATTTGGCGCATTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:3110899/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCgataa                                                   >  1:639376/1‑56 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGc                                         >  1:805312/1‑68 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGc                                         >  1:1007219/1‑68 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGc                                         >  1:806312/1‑68 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGc                                         >  1:2858902/1‑68 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGc                                         >  1:2977784/1‑68 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGc                                         >  1:2102391/1‑68 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGc                                         >  1:693202/1‑68 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:1661884/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:1217487/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:2962849/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:2765993/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:2638214/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:2448284/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:1053433/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:3344275/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:220349/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:3386742/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:3430917/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:1963835/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:1470773/1‑69 (MQ=255)
       tCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCGTCTTCATCGATAATGCTCATGGCa                                        >  1:1852859/1‑69 (MQ=255)
              ccATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATc                                                              <  1:666472/40‑1 (MQ=255)
              ccATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATc                                                              <  1:3245242/40‑1 (MQ=255)
              ccATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATc                                                              <  1:886074/40‑1 (MQ=255)
                   tGGCGCAGTGCGTCTAATGTGGCTTTTGCTTCATCTTCGTCGATGATGCt                                               <  1:1392458/50‑1 (MQ=255)
                            gCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAcc                   >  1:3083184/1‑69 (MQ=255)
                            gCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAcc                   <  1:980342/69‑1 (MQ=255)
                            gCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAcc                   <  1:1329656/69‑1 (MQ=255)
                            gCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAcc                   <  1:2338961/69‑1 (MQ=255)
                            gCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAcc                   <  1:2318728/69‑1 (MQ=255)
                            gCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAcc                   >  1:1868867/1‑69 (MQ=255)
                            gCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAcc                   <  1:2284250/69‑1 (MQ=255)
                             cGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAccc                  >  1:648577/1‑69 (MQ=255)
                             cGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCAccc                  >  1:2230974/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCt                                               >  1:2199543/1‑36 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACt               >  1:1692064/1‑68 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:2837271/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:1549436/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:666968/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:1423202/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:462654/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:3253688/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:3266060/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:2338806/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:3274768/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:80651/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:2564716/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:830087/1‑69 (MQ=255)
                                 tAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTg              >  1:2337516/1‑69 (MQ=255)
                                      tggctttggcttCATCTTCGTCGATGATGCTCGTGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:2922619/1‑68 (MQ=255)
                                      tggctttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:3308715/1‑68 (MQ=255)
                                      tggctttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:166315/1‑68 (MQ=255)
                                      tggctttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:714501/1‑68 (MQ=255)
                                      tggctttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:1787740/1‑68 (MQ=255)
                                      tggctttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:95305/1‑68 (MQ=255)
                                      tggctttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:124352/1‑68 (MQ=255)
                                          tttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:1663803/69‑1 (MQ=255)
                                          tttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:3017610/69‑1 (MQ=255)
                                          tttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:218259/69‑1 (MQ=255)
                                          tttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:2301115/69‑1 (MQ=255)
                                          tttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:693106/69‑1 (MQ=255)
                                          tttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:2852910/69‑1 (MQ=255)
                                          tttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:2908981/69‑1 (MQ=255)
                                           ttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:717200/68‑1 (MQ=255)
                                           ttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:1855410/68‑1 (MQ=255)
                                           ttggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:3145728/68‑1 (MQ=255)
                                           ttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:1330846/1‑63 (MQ=255)
                                           ttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:1408852/1‑63 (MQ=255)
                                           ttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:613511/1‑63 (MQ=255)
                                           ttggcttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGccc          >  1:3206768/1‑63 (MQ=255)
                                            tggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCcagcag     <  1:3230100/67‑1 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGAtg  >  1:929938/1‑68 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:994772/1‑69 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:3321076/1‑69 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:3316086/1‑69 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:2999074/1‑69 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:2458368/1‑69 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:1816961/1‑69 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCCGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:1010172/1‑69 (MQ=255)
                                             ggcttCGTCTTCATCGATAATGCTCATGGCAAATCAGACGTGTACCAGTACCCACTGGCCCAGCAGATc  >  1:68171/1‑69 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     >  1:361274/1‑63 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     >  1:2731376/1‑63 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     >  1:2323532/1‑63 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     >  1:1767798/1‑63 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     <  1:1451476/63‑1 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     <  1:1408976/63‑1 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     >  1:1227390/1‑63 (MQ=255)
                                                ttCATCTTCGTCGATGATGCTAATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAg     <  1:1781644/63‑1 (MQ=255)
                                                   |                                                              
GGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATC  >  minE/2002485‑2002598

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