Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,823,331 G→T 19.3% intergenic (+68/+60) norW → / ← hypF NADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,823,3310GT19.3% 99.2 / 7.8 36intergenic (+68/+60)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base G (7/22);  new base T (2/5);  total (9/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 = 1.00e+00

GGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTC  >  minE/1823267‑1823400
                                                                 |                                                                     
ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAgg                                                                 <  1:3860173/71‑1 (MQ=255)
ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAgg                                                                 <  1:1731577/71‑1 (MQ=255)
   tACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCg                                                              >  1:1062286/1‑71 (MQ=255)
     cTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGa                                                            <  1:2140307/71‑1 (MQ=255)
        tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATa                                                          <  1:2470456/70‑1 (MQ=255)
        tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATa                                                          <  1:850126/70‑1 (MQ=255)
        tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATa                                                          <  1:1366859/70‑1 (MQ=255)
          ccTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAgg                                                       <  1:1193108/71‑1 (MQ=255)
              aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGtt                                                   >  1:523859/1‑71 (MQ=255)
                gtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGttt                                                  <  1:1518391/70‑1 (MQ=21)
                 tGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAc                                                <  1:1564648/71‑1 (MQ=255)
                   tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCAtg                                                                       <  1:77009/46‑1 (MQ=255)
                   tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcg                                               <  1:3682640/70‑1 (MQ=255)
                   tCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGAt                                                           <  1:3729480/59‑1 (MQ=18)
                   tCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcg                                               <  1:3630201/71‑1 (MQ=21)
                      gATTCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCAtttg                                                                     <  1:2979983/46‑1 (MQ=11)
                      gATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCAtttg                                                                     <  1:1020986/46‑1 (MQ=11)
                        tGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                         <  1:2192944/71‑1 (MQ=255)
                        tGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                         >  1:2065098/1‑71 (MQ=255)
                             cGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAgg                                                                 >  1:3869039/1‑42 (MQ=255)
                             cGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                    <  1:2743751/71‑1 (MQ=255)
                             cGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                        >  1:3271078/1‑68 (MQ=18)
                              gCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                   <  1:2999668/71‑1 (MQ=255)
                                 gacgcgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCac                                  >  1:2978610/3‑69 (MQ=18)
                                  acgcgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCacc                                 <  1:3187666/69‑3 (MQ=18)
                                     cgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATa                                                          >  1:3005508/1‑41 (MQ=255)
                                       tCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                  <  1:1946645/63‑1 (MQ=255)
                                           aTCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgt                              <  1:3298407/63‑1 (MQ=255)
                                           aTCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCaa                        <  1:2020963/69‑1 (MQ=255)
                                             ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                            >  1:1546459/1‑63 (MQ=255)
                                                aCCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt                      <  1:3934155/66‑1 (MQ=255)
                                                aCCTACGGGGACGCATGTGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt                      <  1:863585/66‑1 (MQ=255)
                                                   tACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAc                                                <  1:3781593/38‑1 (MQ=14)
                                                    aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCa                                          <  1:3070990/42‑1 (MQ=255)
                                                    aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCaaa                       >  1:3220937/1‑61 (MQ=255)
                                                            cGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATcc      <  1:3725954/71‑1 (MQ=255)
                                                            cGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATcc      <  1:2998487/71‑1 (MQ=255)
                                                             gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCg     <  1:1856761/71‑1 (MQ=255)
                                                               aTGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGtt   <  1:3069867/71‑1 (MQ=255)
                                                                tgtgTAGTCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTc  <  1:174558/71‑1 (MQ=255)
                                                                tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTc  <  1:2321273/71‑1 (MQ=255)
                                                                 |                                                                     
GGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTC  >  minE/1823267‑1823400

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: