Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,170,199 C→T 48.5% R46R (CGC→CGT gdhA → glutamate dehydrogenase, NADP‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,170,1990CT48.5% 3.3 / 58.6 70R46R (CGC→CGTgdhAglutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base C (12/24);  new base T (7/27);  total (19/51)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

GCAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTG  >  minE/1170136‑1170266
                                                               |                                                                   
gcAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGAt                                                                <  1:932175/69‑1 (MQ=255)
gcAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGAt                                                                <  1:2189208/69‑1 (MQ=255)
gcAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGAt                                                                <  1:1054553/69‑1 (MQ=255)
gcAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGAt                                                                <  1:1225099/69‑1 (MQ=255)
gcAAGCCGTTCGTGAAGTATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGAt                                                                <  1:3016157/68‑1 (MQ=255)
  aaGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGt                                                              <  1:2046233/69‑1 (MQ=255)
          cGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTAc                                                        >  1:525626/1‑67 (MQ=255)
          cGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTAc                                                        >  1:2391728/1‑67 (MQ=255)
          cGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTAc                                                        >  1:141625/1‑67 (MQ=255)
          cGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTAc                                                        >  1:1994602/1‑67 (MQ=255)
          cGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTAc                                                        >  1:2648505/1‑67 (MQ=255)
          cGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTAc                                                        >  1:1788722/1‑67 (MQ=255)
          cGTGAAGTAATGACCACACTCTGGCCTTTTCTCGAACAAAATCCCAAATATCGTCAGATGTCATTAc                                                        >  1:602639/1‑67 (MQ=255)
                gTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTACTGGAGCg                                                 <  1:1202320/68‑1 (MQ=255)
                gTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTACTGGAGCg                                                 <  1:2614394/68‑1 (MQ=255)
                gTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTACTGGAGCg                                                 <  1:2828766/68‑1 (MQ=255)
                gTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTACTGGAGCg                                                 <  1:1741466/68‑1 (MQ=255)
                gTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTACTGGAGCg                                                 <  1:2830381/68‑1 (MQ=255)
                 tAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCCAAATATCGTCAGATGTCATTACTGGAGCg                                                 <  1:709445/67‑1 (MQ=255)
                   aTGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTg                                             >  1:484075/1‑69 (MQ=255)
                   aTGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTg                                             >  1:2243489/1‑69 (MQ=255)
                        cacaCTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGa                                        <  1:298435/69‑1 (MQ=255)
                        cacaCTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGa                                        <  1:3127649/69‑1 (MQ=255)
                        cacaCTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGa                                        <  1:276975/69‑1 (MQ=255)
                        cacaCTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGa                                        <  1:3337504/69‑1 (MQ=255)
                        cacaCTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGa                                        <  1:1546615/69‑1 (MQ=255)
                        cacaCTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGa                                        <  1:110627/69‑1 (MQ=255)
                        cacaCTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGa                                        <  1:1335055/69‑1 (MQ=255)
                          cacTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAAc                                      <  1:3249411/69‑1 (MQ=255)
                          cacTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAAc                                      <  1:2438054/69‑1 (MQ=255)
                          cacTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAAc                                      <  1:2129622/69‑1 (MQ=255)
                          cacTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAAc                                      <  1:1619628/69‑1 (MQ=255)
                           acTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAAc                                      <  1:2271957/68‑1 (MQ=255)
                             tctGGCCTTTTCTTTAACAAAATCCCAAATATCGTCAGATGTCATTACTGGAGCg                                                 <  1:1243296/55‑1 (MQ=255)
                                   cTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAgcgc                              >  1:2422019/1‑68 (MQ=255)
                                   cTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAgcgc                              >  1:1862881/1‑68 (MQ=255)
                                                     cAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTgg           <  1:2076667/69‑1 (MQ=255)
                                                     cAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTgg           <  1:2509516/69‑1 (MQ=255)
                                                     cAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTgg           <  1:2730289/69‑1 (MQ=255)
                                                     cAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTgg           <  1:3386739/69‑1 (MQ=255)
                                                     cAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTgg           <  1:2776312/69‑1 (MQ=255)
                                                      aaaaTATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTgg           <  1:468876/68‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:3245315/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:1651554/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:1623525/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:3379030/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:132658/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:1281850/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:764162/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:84152/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:911869/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:1062403/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:321758/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:2163735/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:2054400/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:195119/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:2737896/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:2912476/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:2981991/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:1660934/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:222064/69‑1 (MQ=255)
                                                          tatCGTCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAATTTCGCGTAGTATgg      <  1:3221520/69‑1 (MQ=255)
                                                              gCCAGATGTCCTTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGtt   >  1:931455/1‑68 (MQ=255)
                                                              gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTtg  >  1:1427629/1‑69 (MQ=255)
                                                              gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTtg  >  1:1762327/1‑69 (MQ=255)
                                                              gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTtg  >  1:1934434/1‑69 (MQ=255)
                                                              gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTtg  >  1:63553/1‑69 (MQ=255)
                                                              gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTtg  >  1:653797/1‑69 (MQ=255)
                                                              gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTtg  >  1:1683097/1‑69 (MQ=255)
                                                              gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTtg  >  1:2537460/1‑69 (MQ=255)
                                                               |                                                                   
GCAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTG  >  minE/1170136‑1170266

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: