Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 694,044 T→C 32.5% D199D (GAT→GAC 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,0440TC32.5% 71.6 / 47.0 89D199D (GAT→GAChyaFprotein involved in nickel incorporation into hydrogenase‑1 proteins
Reads supporting (aligned to +/- strand):  ref base T (38/22);  new base C (15/14);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGC  >  minE/693979‑694101
                                                                 |                                                         
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:1910631/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:2362975/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:1785991/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:1860381/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:250439/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:1861363/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:3324717/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:1334089/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:586782/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:719847/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:1211086/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:868518/69‑1 (MQ=255)
gTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:1139349/69‑1 (MQ=255)
 tACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCgg                                                        <  1:2811832/68‑1 (MQ=255)
      aaCCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                                                  >  1:1330084/1‑69 (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:2877106/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:2151829/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:2155610/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:3004169/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:2328235/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:2376249/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:2826724/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:1088841/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:3274832/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:3427225/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:890571/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:1902930/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:1388050/1‑69 (MQ=255)
          ccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                              >  1:130072/1‑69 (MQ=255)
            gACGCGCAGCCGCACATCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:360955/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:2674086/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:1373623/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:3331295/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:3007748/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:848753/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:400596/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:1726148/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:2598218/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:1760756/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:1144259/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:2417633/69‑1 (MQ=255)
            gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                            <  1:2326417/69‑1 (MQ=255)
                    gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTctct                                     >  1:162879/1‑68 (MQ=255)
                    gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg                                    >  1:747118/1‑69 (MQ=255)
                    gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg                                    >  1:2056288/1‑69 (MQ=255)
                    gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg                                    >  1:2642852/1‑69 (MQ=255)
                    gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg                                    >  1:1714688/1‑69 (MQ=255)
                    gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg                                    >  1:1885879/1‑69 (MQ=255)
                    gCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtg                                    >  1:2140024/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:83642/1‑49 (MQ=37)
                             cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                                               >  1:1875369/1‑49 (MQ=37)
                             cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                                               >  1:1364616/1‑49 (MQ=37)
                             cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                                               >  1:2076962/1‑49 (MQ=37)
                             cATTAATCTGACTCAATTACCCATCAGCGAAGCTGACCGGCTTTTtctc                                               >  1:1859793/1‑49 (MQ=37)
                                               aCCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGc                               >  1:466414/1‑47 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  >  1:590710/1‑69 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:738712/69‑1 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:416475/69‑1 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:2882265/69‑1 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:3317556/69‑1 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:320304/69‑1 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:1664980/69‑1 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:2719868/69‑1 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  >  1:1785809/1‑69 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  >  1:1838803/1‑69 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  >  1:1941565/1‑69 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCAATGGc  >  1:1045802/1‑69 (MQ=255)
                                                      aGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTATGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:1299035/69‑1 (MQ=255)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCggg                           >  1:2831162/2‑39 (MQ=37)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:975441/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:921319/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:702531/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:449828/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:3104742/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:1632107/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:2131597/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:2118493/2‑40 (MQ=38)
                                                           agCTGACCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGa                          >  1:1961590/2‑40 (MQ=38)
                                                              tGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGc  <  1:2478130/61‑1 (MQ=255)
                                                                 |                                                         
GTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTACCATTGGC  >  minE/693979‑694101

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