Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,823,325:1 +G 23.7% intergenic (+62/+66) 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,3251.G23.7% 95.6 / 42.9 38intergenic (+62/+66)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base . (9/20);  new base G (2/7);  total (11/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 = 9.37e-01

ATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAT  >  minE/1823258‑1823380
                                                                    |                                                       
aTGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGc                                                      <  1:2478489/71‑1 (MQ=255)
 tGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCa                                                     >  1:2126560/1‑71 (MQ=255)
         ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAgg                                             <  1:3860173/71‑1 (MQ=255)
         ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAgg                                             <  1:1731577/71‑1 (MQ=255)
            tACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCg                                          >  1:1062286/1‑71 (MQ=255)
              cTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGa                                        <  1:2140307/71‑1 (MQ=255)
                 tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATa                                      <  1:850126/70‑1 (MQ=255)
                 tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATa                                      <  1:1366859/70‑1 (MQ=255)
                 tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATa                                      <  1:2470456/70‑1 (MQ=255)
                   ccTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAgg                                   <  1:1193108/71‑1 (MQ=255)
                       aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGtt                               >  1:523859/1‑71 (MQ=255)
                         gtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGttt                              <  1:1518391/70‑1 (MQ=21)
                          tGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTAc                            <  1:1564648/71‑1 (MQ=255)
                            tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCAtg                                                   <  1:77009/46‑1 (MQ=255)
                            tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTAcg                           <  1:3682640/70‑1 (MQ=255)
                            tCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGAt                                       <  1:3729480/59‑1 (MQ=18)
                            tCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcg                           <  1:3630201/71‑1 (MQ=21)
                              ggATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑Cg                                                       >  1:3024933/1‑40 (MQ=255)
                               gATTCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCAtttg                                                 <  1:2979983/46‑1 (MQ=11)
                               gATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCAtttg                                                 <  1:1020986/46‑1 (MQ=11)
                                 tGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCAt                     <  1:2192944/71‑1 (MQ=255)
                                 tGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCAt                     >  1:2065098/1‑71 (MQ=255)
                                      cGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAgg                                             >  1:3869039/1‑42 (MQ=255)
                                      cGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                <  1:2743751/71‑1 (MQ=255)
                                      cGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                    >  1:3271078/1‑68 (MQ=18)
                                       gCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa               <  1:2999668/71‑1 (MQ=255)
                                          gacgcgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCac              >  1:2978610/3‑69 (MQ=18)
                                           acgcgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCacc             <  1:3187666/69‑3 (MQ=18)
                                              cgTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATa                                      >  1:3005508/1‑41 (MQ=255)
                                                tCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa              <  1:1946645/63‑1 (MQ=255)
                                                    aTCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgt          <  1:3298407/63‑1 (MQ=255)
                                                    aTCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCaa    <  1:2020963/69‑1 (MQ=255)
                                                      ccGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt        >  1:1546459/1‑63 (MQ=255)
                                                         aCCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt  <  1:3934155/66‑1 (MQ=255)
                                                         aCCTACGGGGA‑CGCATGTGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt  <  1:863585/66‑1 (MQ=255)
                                                            tACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAc                            <  1:3781593/38‑1 (MQ=14)
                                                             aCGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCa                      <  1:3070990/42‑1 (MQ=255)
                                                             aCGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCaaa   >  1:3220937/1‑61 (MQ=255)
                                                                    |                                                       
ATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAT  >  minE/1823258‑1823380

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: ‑