Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,823,331 G→T 19.0% 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.0% 140.9 / 14.0 52intergenic (+68/+60)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base T (4/6);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTT  >  minE/1823261‑1823399
                                                                       |                                                                    
tAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCAtg                                                                      >  1:5614669/1‑71 (MQ=255)
  ggTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCAtgtg                                                                    <  1:2626366/71‑1 (MQ=255)
   gTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATgtgt                                                                   <  1:1008766/71‑1 (MQ=255)
   gTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATgtgt                                                                   <  1:5467348/71‑1 (MQ=255)
      ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAgg                                                                <  1:2687919/71‑1 (MQ=255)
      ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAgg                                                                >  1:3997628/1‑71 (MQ=255)
           cTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGc                                                               >  1:3974633/1‑67 (MQ=255)
           cTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGa                                                           <  1:2086222/71‑1 (MQ=255)
              tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGAt                                                          <  1:60408/69‑1 (MQ=255)
              tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATaa                                                        <  1:2136262/71‑1 (MQ=255)
               gCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAg                                                       >  1:5442851/1‑71 (MQ=255)
                 cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCAtgt                                                                     <  1:4765372/55‑1 (MQ=255)
                 cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGc                                                     <  1:839346/71‑1 (MQ=255)
                 cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGc                                                     <  1:61392/71‑1 (MQ=255)
                   aaaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATgtgt                                                                   <  1:4665323/55‑1 (MQ=255)
                    gggtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGc                                                     <  1:3300383/66‑1 (MQ=18)
                    gggtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCg                                                    <  1:5224143/67‑1 (MQ=18)
                    aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGtt                                                  <  1:2659720/71‑1 (MQ=255)
                    aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGtt                                                  >  1:3545806/1‑71 (MQ=255)
                    aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGt                                                   <  1:5583964/70‑1 (MQ=255)
                     ggtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGtt                                                  <  1:2093581/69‑1 (MQ=21)
                      aTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTa                                                >  1:3678413/1‑71 (MQ=255)
                       tGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAc                                               >  1:4195416/1‑71 (MQ=255)
                         tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcg                                              <  1:226097/70‑1 (MQ=255)
                         tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcg                                              <  1:1206975/70‑1 (MQ=255)
                          cGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcgcc                                            >  1:1009605/1‑71 (MQ=255)
                           ggATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcgcc                                            <  1:1886594/70‑1 (MQ=255)
                            gATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccg                                           <  1:3502571/71‑1 (MQ=18)
                            gATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccg                                           <  1:2220199/71‑1 (MQ=18)
                             aTGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCg                                                             >  1:106473/1‑51 (MQ=255)
                               gCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                        <  1:342817/70‑1 (MQ=255)
                                cGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                       <  1:5260542/70‑1 (MQ=255)
                                cGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                       <  1:2446669/70‑1 (MQ=255)
                                 gACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                      >  1:4586736/1‑71 (MQ=18)
                                            gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTc                          >  1:3734995/1‑71 (MQ=255)
                                            gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTc                          >  1:4468053/1‑71 (MQ=255)
                                               ttATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                >  1:2486958/1‑62 (MQ=255)
                                               ttATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                >  1:3773046/1‑62 (MQ=255)
                                                  tCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg                    >  1:1056081/1‑71 (MQ=255)
                                                   ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                           >  1:381741/1‑63 (MQ=255)
                                                    cGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGtt                                                  <  1:521248/40‑1 (MQ=17)
                                                    cGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc                   <  1:2103605/70‑1 (MQ=255)
                                                    cGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc                   <  1:1676696/70‑1 (MQ=255)
                                                    cGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGCGTCCAAATGc                   <  1:150420/70‑1 (MQ=255)
                                                       ccTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTa                                                >  1:170489/1‑39 (MQ=14)
                                                        cTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAacac              >  1:3485102/1‑71 (MQ=255)
                                                          aCGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgcc                                            >  1:617926/1‑40 (MQ=11)
                                                          aCGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgcc                                            >  1:2492745/1‑40 (MQ=11)
                                                          aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGt            >  1:1425518/1‑71 (MQ=255)
                                                           cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                   >  1:1035101/1‑47 (MQ=255)
                                                           cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                   >  1:4126872/1‑47 (MQ=255)
                                                           cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACg             <  1:4216816/69‑1 (MQ=255)
                                                                   gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc                   <  1:2547528/56‑1 (MQ=255)
                                                                   gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc                   <  1:5193773/56‑1 (MQ=255)
                                                                   gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCAACCGGCAATGGTGTCCAAATGCAACACGTTTTATCCg    >  1:1937331/1‑71 (MQ=255)
                                                                     aTGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATc      >  1:364538/1‑67 (MQ=255)
                                                                     aTGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGtt  >  1:2022711/1‑71 (MQ=255)
                                                                     aTGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGt   >  1:5312209/1‑70 (MQ=255)
                                                                       |                                                                    
TAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTT  >  minE/1823261‑1823399

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: