Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 694,025 T→A 7.7% L193* (TTA→TAA)  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,0250TA 7.7% 181.2 / 7.0 65L193* (TTA→TAA) hyaFprotein involved in nickel incorporation into hydrogenase‑1 proteins
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base A (3/2);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

GCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTAC  >  minE/693960‑694094
                                                                 |                                                                     
gaaTGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTAc                                                                     >  1:1018324/3‑68 (MQ=255)
   tGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTa                                                                      >  1:3603487/1‑64 (MQ=255)
   tGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCAt                                                                 >  1:993049/1‑69 (MQ=255)
   tGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCAt                                                                 >  1:4209001/1‑69 (MQ=255)
        tACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGc                                                             >  1:2406717/1‑68 (MQ=255)
        tACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGc                                                             >  1:3431197/1‑68 (MQ=255)
        tACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGc                                                             >  1:1819717/1‑68 (MQ=255)
        tACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGc                                                             >  1:182797/1‑68 (MQ=255)
         aCTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAgg                                                         >  1:1266944/1‑71 (MQ=255)
         aCTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAg                                                          <  1:2929394/70‑1 (MQ=255)
         aCTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAg                                                          <  1:485524/70‑1 (MQ=255)
               acaTGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATc                                                   >  1:758742/1‑71 (MQ=255)
               acaTGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATc                                                   >  1:3841775/1‑71 (MQ=255)
               acaTGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATc                                                   >  1:3778880/1‑71 (MQ=255)
                  tGTACGTAACCCCGACGCGCATCCGCACAGCATTAATCTGACGCAATTa                                                                      <  1:1404315/49‑1 (MQ=255)
                  tGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTg                                                      <  1:2312783/65‑1 (MQ=255)
                    tACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGAt                                                    >  1:3507128/1‑65 (MQ=255)
                    tACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGAt                                                    >  1:432234/1‑65 (MQ=255)
                    tACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGAt                                                    >  1:406104/1‑65 (MQ=255)
                    tACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCtt                                              >  1:3010933/1‑71 (MQ=255)
                    tACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCtt                                              >  1:1175800/1‑71 (MQ=255)
                      cGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTg                                                      <  1:287444/61‑1 (MQ=255)
                      cGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCt                                               <  1:565606/68‑1 (MQ=255)
                      cGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCt                                               <  1:2542835/68‑1 (MQ=255)
                       gTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGa                                                           <  1:4182692/55‑1 (MQ=255)
                       gTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                                           <  1:2078415/71‑1 (MQ=255)
                       gTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                                           <  1:114993/71‑1 (MQ=255)
                       gTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                                           <  1:2611395/71‑1 (MQ=255)
                       gTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCt                                               <  1:970226/67‑1 (MQ=255)
                       gTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGCCttttt                                           <  1:163212/71‑1 (MQ=255)
                        tAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCttt                                             >  1:2272770/1‑68 (MQ=255)
                        tAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtc                                          <  1:3993553/71‑1 (MQ=255)
                        tAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtc                                          <  1:2039706/71‑1 (MQ=255)
                        tAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtc                                          <  1:1221620/71‑1 (MQ=255)
                        tAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATAACCCATCAGCGAGGCTGATCGGCTTTTtc                                          <  1:3948967/71‑1 (MQ=255)
                        tAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATAACCCATCAGCGAGGCTGATCGGCTTTTtc                                          <  1:3821746/71‑1 (MQ=255)
                           ccccGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                       <  1:3672447/71‑1 (MQ=255)
                               gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                     >  1:1562969/1‑69 (MQ=255)
                               gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                     >  1:3204616/1‑69 (MQ=255)
                               gACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCa                                     >  1:1546028/1‑69 (MQ=255)
                                   cgcAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGtctc                               <  1:2256891/71‑1 (MQ=255)
                                             caGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGtctc                               <  1:1970266/61‑1 (MQ=255)
                                                 aTTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTtctct                                       <  1:2481856/49‑1 (MQ=255)
                                                 aTTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtgt                            <  1:1127775/60‑1 (MQ=255)
                                                 aTTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtgt                            <  1:3552341/60‑1 (MQ=255)
                                                  ttAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCtgt                            <  1:3675541/59‑1 (MQ=255)
                                                    aaTCTGACGCAATTACCCATCAGCGAGGCTGATCGGCtttt                                            >  1:1153475/1‑41 (MQ=255)
                                                    aaTCTGACGCAATTACCCATCAGCGAGGCTGATCGGCtttt                                            >  1:1755178/1‑41 (MQ=255)
                                                    aaTCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAttat              >  1:2626963/1‑71 (MQ=255)
                                                    aaTCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAatat              >  1:3846800/1‑71 (MQ=255)
                                                    aaTCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAatat              >  1:940315/1‑71 (MQ=255)
                                                        tGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                                           >  1:3043602/1‑38 (MQ=255)
                                                        tGACGCAATTACCCATCAGCGAGGCTGATCGGCttttt                                           >  1:1532564/1‑38 (MQ=255)
                                                        tGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATAtt             >  1:516745/1‑68 (MQ=255)
                                                        tGACGCAATAACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATAtt             >  1:1516013/1‑68 (MQ=255)
                                                        tGACGCAATAACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATAtt             >  1:2712543/1‑68 (MQ=255)
                                                        tGACGCAATAACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATAtt             >  1:2288556/1‑68 (MQ=255)
                                                         gACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATAtt             >  1:1622142/1‑67 (MQ=255)
                                                            gCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTc      >  1:3155778/1‑71 (MQ=255)
                                                            gCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTc      >  1:726336/1‑71 (MQ=255)
                                                            gCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTc      >  1:2938781/1‑71 (MQ=255)
                                                            gCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTc      >  1:608579/1‑71 (MQ=255)
                                                            gCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTc      >  1:1135076/1‑71 (MQ=255)
                                                             cAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCg     <  1:1017662/71‑1 (MQ=255)
                                                             cAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCg     <  1:292288/71‑1 (MQ=255)
                                                                ttACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTAc  <  1:3725214/71‑1 (MQ=255)
                                                                ttACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTAc  >  1:2636081/1‑71 (MQ=255)
                                                                 tACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTa   >  1:2845051/1‑69 (MQ=255)
                                                                 tACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTAc  >  1:538600/1‑70 (MQ=255)
                                                                 tACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTAc  >  1:1845386/1‑70 (MQ=255)
                                                                 tACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTAc  >  1:1419695/1‑70 (MQ=255)
                                                                 |                                                                     
GCATGAGTTACTGGCACATGTACGTAACCCCGACGCGCAGCCGCACAGCATTAATCTGACGCAATTACCCATCAGCGAGGCTGATCGGCTTTTTCTCTCACGTCTCTGTGGGCCGGGAAATATTCAGATTCGTAC  >  minE/693960‑694094

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: