Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,919,248 G→A 7.4% Q205Q (CAG→CAA tas → predicted oxidoreductase, NADP(H)‑dependent aldo‑keto reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,919,2480GA 7.4% 259.1 / 7.9 81Q205Q (CAG→CAAtaspredicted oxidoreductase, NADP(H)‑dependent aldo‑keto reductase
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base A (4/2);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CTGCATTTGGCGTAATGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATGAAGGGGT  >  minE/1919178‑1919316
                                                                      |                                                                    
cTGCATTTGGCGTAATGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAg                                                                      <  1:2803764/71‑1 (MQ=255)
 tGCATTTGGCGTAATGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAg                                                                      >  1:1274206/1‑70 (MQ=255)
       tGGCGTAATGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCt                                                                >  1:1192451/1‑70 (MQ=255)
       tGGCGTAATGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCt                                                                >  1:348601/1‑70 (MQ=255)
          cGTAATGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACa                                                            <  1:1464185/71‑1 (MQ=255)
              aTGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCt                                                        >  1:2189466/1‑71 (MQ=255)
              aTGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCt                                                        >  1:2190759/1‑71 (MQ=255)
              aTGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCt                                                        >  1:1352023/1‑71 (MQ=255)
              aTGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCt                                                        >  1:2945940/1‑71 (MQ=255)
                gcgcTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACa                                                            >  1:1681486/1‑65 (MQ=255)
                gcgcTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACa                                                            >  1:672148/1‑65 (MQ=255)
                 cgcTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGt                                                      >  1:2165623/1‑70 (MQ=255)
                 cgcTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGt                                                      >  1:933319/1‑70 (MQ=255)
                 cgcTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAAAACCCTTACAGTCTGt                                                      >  1:1075667/1‑70 (MQ=255)
                   cTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTaa                                                   <  1:1148983/71‑1 (MQ=255)
                        tGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCt                                                        >  1:2251939/1‑61 (MQ=255)
                        tGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTaaa                                                  <  1:2313461/67‑1 (MQ=255)
                        tGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAAc                                                 >  1:599156/1‑68 (MQ=255)
                         gCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCg                                               >  1:268449/1‑69 (MQ=255)
                         gCATCTGGCGGACAAACACGATCTGACGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCg                                               >  1:2361501/1‑69 (MQ=255)
                            tCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTa                                                    <  1:505989/61‑1 (MQ=255)
                              tGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGtt                                          >  1:2790697/1‑69 (MQ=255)
                               ggCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAAAACCCTTACAGTCTGTTAAACCGCAg                                            <  1:3054044/66‑1 (MQ=255)
                                gCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAAc                                                 >  1:3422590/1‑60 (MQ=255)
                                gCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTg                                       <  1:2515201/70‑1 (MQ=255)
                                gCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTg                                       <  1:205389/70‑1 (MQ=255)
                                 cGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGaa                                     <  1:774165/71‑1 (MQ=255)
                                  ggACAAACACGATCTGCCGCGTATTGTCACCATTCAGAAc                                                                   <  1:960069/40‑1 (MQ=255)
                                   gACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTg                                       <  1:2060279/67‑1 (MQ=255)
                                   gACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGt                                   <  1:4093384/71‑1 (MQ=255)
                                     cAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGttt                                         >  1:4116732/1‑63 (MQ=255)
                                     cAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGt                                           <  1:273396/61‑1 (MQ=255)
                                      aaaCACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAgg                                <  1:3675146/71‑1 (MQ=255)
                                       aaCACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTAc                                                             <  1:1693907/41‑1 (MQ=255)
                                       aaCACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTAc                                                             <  1:856205/41‑1 (MQ=255)
                                        acacGATCTGCCGCGTATTGTCACCATTCAGAACCCtt                                                               <  1:1400396/38‑1 (MQ=255)
                                         cacGATCTGCCGCGTATTGTCACCATTCAGAACCCTTa                                                              >  1:1176094/1‑38 (MQ=255)
                                         cacGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGtt                                                     >  1:2519426/1‑47 (MQ=255)
                                         cacGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTaaa                                                  >  1:720925/1‑50 (MQ=255)
                                         cacGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTaaa                                                  >  1:159509/1‑50 (MQ=255)
                                         cacGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTaaa                                                  >  1:1667785/1‑50 (MQ=255)
                                         cacGATCTGCCGCGTATTGTCACCATTCAAAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCt                             >  1:1674363/1‑71 (MQ=255)
                                             aTCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGtttt                                        <  1:231049/56‑1 (MQ=255)
                                             aTCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGt                                           <  1:941752/53‑1 (MQ=255)
                                               cTGCCGCGTATTGTCACCATTCAGAACCCTTACAg                                                           >  1:3520672/1‑35 (MQ=255)
                                                  ccGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAg                        <  1:1795153/67‑1 (MQ=255)
                                                  ccGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAg                        <  1:2965657/67‑1 (MQ=255)
                                                  ccGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAg                     <  1:3957901/70‑1 (MQ=255)
                                                  ccGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGt                    >  1:2117798/1‑71 (MQ=255)
                                                  ccGCGTATTGTCACCATTCAAAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAg                        <  1:262363/67‑1 (MQ=255)
                                                   cgcgTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAg                                    >  1:3382998/1‑54 (MQ=255)
                                                    gcgTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCa                  <  1:1456395/71‑1 (MQ=255)
                                                    gcgTATTGTCACCATTCAAAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTa                                  >  1:1860371/1‑55 (MQ=255)
                                                    gcgTATTGTCACCATTCAAAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTa                                  >  1:3027881/1‑55 (MQ=255)
                                                      gTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGt                                   <  1:367577/52‑1 (MQ=255)
                                                      gTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGc                <  1:2843355/71‑1 (MQ=255)
                                                       tATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGc                >  1:2056869/1‑70 (MQ=255)
                                                         ttGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGt                    >  1:2656004/1‑64 (MQ=255)
                                                         ttGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCa              >  1:1177571/1‑70 (MQ=255)
                                                         ttGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAg             <  1:2491931/71‑1 (MQ=255)
                                                          tGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAg             >  1:3586519/1‑70 (MQ=255)
                                                           gTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAg                                 >  1:2583454/1‑49 (MQ=255)
                                                           gTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTa           >  1:1603447/1‑71 (MQ=255)
                                                            tCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGaa                                     >  1:2226197/1‑44 (MQ=255)
                                                            tCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGa                                      <  1:651859/43‑1 (MQ=255)
                                                            tCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTAt          >  1:3093933/1‑71 (MQ=255)
                                                              aCCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAg                        <  1:3568470/55‑1 (MQ=255)
                                                              aCCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATg         <  1:4114133/70‑1 (MQ=255)
                                                               ccATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGt                               <  1:3429933/47‑1 (MQ=255)
                                                               ccATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGt                               <  1:545102/47‑1 (MQ=255)
                                                               ccATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATg         >  1:544132/1‑69 (MQ=255)
                                                               ccATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATGaa       <  1:1225654/71‑1 (MQ=255)
                                                               ccATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGACAGCCAGTATGa        <  1:3983659/70‑1 (MQ=255)
                                                               ccATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTAGAAGTAGGTCTGGCAGAAGTCAGCCAGTATg         >  1:589816/1‑69 (MQ=255)
                                                                cATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCt                             <  1:1562552/48‑1 (MQ=255)
                                                                 aTTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATGAAgg     >  1:1984431/1‑71 (MQ=255)
                                                                  ttCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTg                                       <  1:868389/36‑1 (MQ=255)
                                                                  ttCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAg                                    >  1:756432/1‑39 (MQ=255)
                                                                  ttCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAg                        <  1:1031413/51‑1 (MQ=255)
                                                                  ttCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATGAAggg    <  1:3315114/71‑1 (MQ=255)
                                                                   tCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGc                >  1:3911131/1‑58 (MQ=255)
                                                                    cAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATg         >  1:3666834/1‑64 (MQ=255)
                                                                    cAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATGAAGGGGt  >  1:2369806/1‑71 (MQ=255)
                                                                      |                                                                    
CTGCATTTGGCGTAATGCGCTACCTGCATCTGGCGGACAAACACGATCTGCCGCGTATTGTCACCATTCAGAACCCTTACAGTCTGTTAAACCGCAGTTTTGAAGTAGGTCTGGCAGAAGTCAGCCAGTATGAAGGGGT  >  minE/1919178‑1919316

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: