Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,889,980 G→A 24.4% V83V (GTC→GTT gudX ← predicted glucarate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,889,9800GA24.4% 80.8 / 19.7 65V83V (GTC→GTTgudXpredicted glucarate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (30/19);  new base A (11/5);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAACCATCGGAATAGCATCGACAAGC  >  minE/1889913‑1890040
                                                                   |                                                            
cAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTgg                                                              >  1:3194024/1‑68 (MQ=255)
    ggTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCAc                                                         >  1:2693080/1‑69 (MQ=255)
    ggTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCAc                                                         >  1:1845274/1‑69 (MQ=255)
    ggTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCAc                                                         >  1:1829129/1‑69 (MQ=255)
    ggTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCAc                                                         >  1:1705266/1‑69 (MQ=255)
            cACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCa                                                 >  1:3048498/1‑69 (MQ=255)
              ccTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAg                                                >  1:1298162/1‑68 (MQ=255)
              ccTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGt                                               >  1:121495/1‑69 (MQ=255)
              ccTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGt                                               >  1:2620724/1‑69 (MQ=255)
              ccTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGt                                               >  1:517931/1‑69 (MQ=255)
              ccTTTACCGAAGGGATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGt                                               >  1:2954271/1‑69 (MQ=255)
               cTTTGCCGAAGGTATCAAAATCGGCTGCCTGATTACCTTTATGCACCTGCTGAACCACTTTATTCAggc                                              <  1:2894705/69‑3 (MQ=255)
               cTTTGCCGAAGGTATCAAAATCGGCTGCCTGATTACCTTTATGCACCTGCTGAACCACTTTATTCAggc                                              <  1:2109756/69‑3 (MQ=255)
               cTTTGCCGAAGGTATCAAAATCGGCTGCCTGATTACCTTTATGCACCTGCTGAACCACTTTATTCAggc                                              <  1:2047569/69‑3 (MQ=255)
               cTTTGCCGAAGGTATCAAAATCGGCTGCCTGATTACCTTTATGCACCTGCTGAACCACTTTATTCAggc                                              <  1:1720740/69‑3 (MQ=255)
               cTTTGCCGAAGGTATCAAAATCGGCTGCCTGATTACCTTTATGCACCTGCTGAACCACTTTATTCAggc                                              <  1:865937/69‑3 (MQ=255)
                tttACCGAAGGTATCAAAATTGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:2102713/1‑69 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTgg                                                              >  1:3019990/1‑52 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGAcc                                                           >  1:3023961/1‑55 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:2766431/1‑69 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:329358/1‑69 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:2663852/1‑69 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:2975036/1‑69 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:2433386/1‑69 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:18946/1‑69 (MQ=255)
                tttACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTcg                                             >  1:3288732/1‑69 (MQ=255)
                             tCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:327141/1‑69 (MQ=255)
                             tCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:2142947/1‑69 (MQ=255)
                             tCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:328387/1‑69 (MQ=255)
                             tCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:1352474/1‑69 (MQ=255)
                             tCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:1680222/1‑69 (MQ=255)
                               aaaaTCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:2351237/1‑67 (MQ=255)
                               aaaaTCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:1322797/1‑67 (MQ=255)
                               aaaaTCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:2573658/1‑67 (MQ=255)
                               aaaaTCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:122658/1‑67 (MQ=255)
                               aaaaTCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:1537825/1‑67 (MQ=255)
                               aaaaTCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTgg                                >  1:2762513/1‑67 (MQ=255)
                                   tCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:321450/69‑1 (MQ=255)
                                   tCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:47070/69‑1 (MQ=255)
                                   tCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:374034/69‑1 (MQ=255)
                                   tCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:659062/69‑1 (MQ=255)
                                   tCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:1876174/69‑1 (MQ=255)
                                   tCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:783945/69‑1 (MQ=255)
                                   tCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:1178957/69‑1 (MQ=255)
                                    cggcggCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGa                          <  1:2572450/68‑1 (MQ=255)
                                       cggcCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAAcc                       <  1:2166006/68‑1 (MQ=255)
                                       cggcCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAAcc                       <  1:3032271/68‑1 (MQ=255)
                                       cggcCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAAcc                       <  1:2129564/68‑1 (MQ=255)
                                       cggcCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAAcc                       <  1:1874592/68‑1 (MQ=255)
                                       cggcCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAAcc                       <  1:331309/68‑1 (MQ=255)
                                        ggcCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAAcc                       <  1:3358900/67‑1 (MQ=255)
                                           cTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGCCCCAGAACCATCgg                  <  1:1305043/69‑1 (MQ=255)
                                                    tttATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAACCATCGGAATAGCATc         <  1:279827/69‑1 (MQ=255)
                                                         gCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAAcc                       <  1:684906/50‑1 (MQ=255)
                                                          cACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAACCATCGGAATAGCATCGACaa    <  1:1229372/68‑1 (MQ=255)
                                                          cACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAACCATCGGAATAGCATCGACaa    <  1:1107622/68‑1 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTg                                 >  1:1755263/1‑38 (MQ=38)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:377226/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:902915/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:356606/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:3346508/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:3248914/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:2810242/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:2133437/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:17784/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:1253292/1‑68 (MQ=255)
                                                           aCCTGCTGAACCACTTTATTCAGGCGCGCAACTTCCTGGCCCAGCACCATCGGAATAGCCTTGACCAgt  >  1:113033/1‑68 (MQ=255)
                                                                   |                                                            
CAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCTGGACCACTTTATTCAGTCGCGCAACTTCCTGGCCCAGAACCATCGGAATAGCATCGACAAGC  >  minE/1889913‑1890040

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: