Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,823,325:1 +G 17.5% 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.G17.5% 175.0 / 45.8 57intergenic (+62/+66)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base . (20/27);  new base G (4/6);  total (24/33)
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.35e-01

GATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGT  >  minE/1823257‑1823389
                                                                     |                                                                
gATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACg                                                                <  1:2565379/71‑1 (MQ=255)
gATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACg                                                                <  1:5339693/71‑1 (MQ=255)
 aTGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGc                                                               <  1:3050926/71‑1 (MQ=255)
 aTGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGc                                                               <  1:5402383/71‑1 (MQ=255)
   gTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCAt                                                             >  1:555788/1‑71 (MQ=255)
    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:3997628/1‑71 (MQ=255)
          ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAgg                                                      <  1:2687919/71‑1 (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:61392/71‑1 (MQ=255)
                     cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGc                                           <  1:839346/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:3773046/1‑62 (MQ=255)
                                                   ttATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                      >  1:2486958/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:2492745/1‑40 (MQ=11)
                                                              aCGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgcc                                  >  1:617926/1‑40 (MQ=11)
                                                              aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGt  >  1:1425518/1‑71 (MQ=255)
                                                               cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                         >  1:4126872/1‑47 (MQ=255)
                                                               cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                         >  1:1035101/1‑47 (MQ=255)
                                                               cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACg   <  1:4216816/69‑1 (MQ=255)
                                                                     |                                                                
GATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGT  >  minE/1823257‑1823389

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: