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

ATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAT  >  minE/1823258‑1823380
                                                                    |                                                       
aTGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGc                                                      <  1:2478489/71‑1 (MQ=255)
 tGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCa                                                     >  1:2126560/1‑71 (MQ=255)
         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:850126/70‑1 (MQ=255)
                 tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATa                                      <  1:1366859/70‑1 (MQ=255)
                 tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATa                                      <  1:2470456/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)
                              ggATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACg                                                       >  1:3024933/1‑40 (MQ=255)
                               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)
                                                                    |                                                       
ATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAT  >  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: