Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 694,020 G→T 24.0% T191T (ACG→ACT hyaF → protein involved in nickel incorporation into hydrogenase‑1 proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE694,0200GT24.0% 137.0 / 32.2 79T191T (ACG→ACThyaFprotein involved in nickel incorporation into hydrogenase‑1 proteins
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base T (5/14);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 3.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-02

GTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCG  >  minE/693954‑694074
                                                                  |                                                      
gTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTTACGc                                                       <  1:47364/68‑1 (MQ=255)
gTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:1152764/68‑1 (MQ=255)
gTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:664998/68‑1 (MQ=255)
gTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:2350345/68‑1 (MQ=255)
gTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:1505152/68‑1 (MQ=255)
gTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:2949983/68‑1 (MQ=255)
   ggagcATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:2983252/62‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:886221/65‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:2015469/65‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:1764768/65‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:2191243/65‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:2855356/65‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:2939093/65‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGc                                                       <  1:1346896/65‑1 (MQ=255)
   ggCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAAtt                                                   >  1:464029/1‑69 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:3324717/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:250439/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:2362975/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:586782/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:719847/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:868518/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:1139349/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:1910631/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:1861363/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:1860381/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:1334089/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:1211086/69‑1 (MQ=255)
                         gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:1785991/69‑1 (MQ=255)
                          tACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                             <  1:2811832/68‑1 (MQ=255)
                               aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                       >  1:1970939/1‑69 (MQ=255)
                               aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                       >  1:3364214/1‑69 (MQ=255)
                               aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                       >  1:3314201/1‑69 (MQ=255)
                               aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                       >  1:3092119/1‑69 (MQ=255)
                               aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                       >  1:2993736/1‑69 (MQ=255)
                               aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                       >  1:1330084/1‑69 (MQ=255)
                               aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                       >  1:2877106/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:3274832/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:3427225/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:3004169/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:1388050/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:2826724/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:130072/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:890571/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:1902930/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:2376249/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:1088841/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:2328235/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:2151829/1‑69 (MQ=255)
                                   ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                   >  1:2155610/1‑69 (MQ=255)
                                     gACGCGCAGCCGCACATCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:360955/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:848753/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:400596/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:3331295/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:1144259/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:3007748/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:1760756/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:1373623/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:2326417/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:1726148/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:2674086/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:2598218/69‑1 (MQ=255)
                                     gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                 <  1:2417633/69‑1 (MQ=255)
                                             gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTctct          >  1:162879/1‑68 (MQ=255)
                                             gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg         >  1:747118/1‑69 (MQ=255)
                                             gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg         >  1:2140024/1‑69 (MQ=255)
                                             gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg         >  1:2056288/1‑69 (MQ=255)
                                             gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg         >  1:1714688/1‑69 (MQ=255)
                                             gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg         >  1:2642852/1‑69 (MQ=255)
                                             gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg         >  1:1885879/1‑69 (MQ=255)
                                                     gCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCg  >  1:1293087/1‑68 (MQ=255)
                                                     gCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCg  >  1:15786/1‑68 (MQ=255)
                                                     gCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCg  >  1:3199430/1‑68 (MQ=255)
                                                     gCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCg  >  1:3121724/1‑68 (MQ=255)
                                                     gCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCg  >  1:1635813/1‑68 (MQ=255)
                                                      cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                    >  1:1364616/1‑49 (MQ=37)
                                                      cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                    >  1:83642/1‑49 (MQ=37)
                                                      cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                    >  1:2076962/1‑49 (MQ=37)
                                                      cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                    >  1:1859793/1‑49 (MQ=37)
                                                      cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                    >  1:1875369/1‑49 (MQ=37)
                                                                  |                                                      
GTTGGCGCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCG  >  minE/693954‑694074

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