Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,154,282 C→T 39.1% Q111Q (CAG→CAA astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,154,2820CT39.1% 23.5 / 31.4 51Q111Q (CAG→CAAastEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (9/22);  new base T (8/12);  total (17/34)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

ACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTC  >  minE/1154233‑1154333
                                                 |                                                   
aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACa                                  <  1:680316/69‑1 (MQ=255)
aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACa                                  <  1:1897368/69‑1 (MQ=255)
aCTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACa                                  <  1:3233989/69‑1 (MQ=255)
        aGTTCGCGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:1820818/1‑64 (MQ=255)
        aGTTCGCGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:2462028/1‑64 (MQ=255)
        aGTTCGCGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:2389062/1‑64 (MQ=255)
        aGTTCGCGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:443421/1‑64 (MQ=255)
        aGTTCGCGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:1030149/1‑64 (MQ=255)
        aGTTCGCGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:3281070/1‑64 (MQ=255)
        aGTTCGCGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:118269/1‑64 (MQ=255)
        aGTTCACGCGCCCGACAGGTTTCCCCGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTc                               >  1:546486/1‑64 (MQ=255)
          ttCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCat                        <  1:2649894/69‑1 (MQ=255)
               gcgcCCGACAGTTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGc                   <  1:1629378/69‑1 (MQ=255)
               gcgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGc                   <  1:1158832/69‑1 (MQ=255)
                cgcCCGAGAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:2845555/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:449982/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:108139/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:620180/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:2578673/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:3293734/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:819578/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:2373936/69‑1 (MQ=255)
                cgcCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:1712802/69‑1 (MQ=255)
                  cccGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTAt                >  1:1395677/1‑69 (MQ=255)
                  cccGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTAt                >  1:2188277/1‑69 (MQ=255)
                  cccGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTAt                >  1:3270051/1‑69 (MQ=255)
                  cccGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTAt                >  1:3196990/1‑69 (MQ=255)
                  cccGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTAt                >  1:1743735/1‑69 (MQ=255)
                  cccGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTAt                >  1:2971242/1‑69 (MQ=255)
                  cccGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTAt                >  1:3094922/1‑69 (MQ=255)
                      aCAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACa                                  <  1:253967/47‑1 (MQ=255)
                      aCAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCt                  <  1:2016913/63‑1 (MQ=255)
                          gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATaa        <  1:381129/69‑1 (MQ=255)
                          gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATaa        <  1:2519827/69‑1 (MQ=255)
                          gTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATCGCAATaa        <  1:1668385/69‑1 (MQ=255)
                               tcCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGtt    <  1:3034708/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:826779/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:996418/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:3421224/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:3349418/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:3336030/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:3091488/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:2961620/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:2900568/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:2818303/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:2461468/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTc  <  1:156092/68‑1 (MQ=255)
                                 cGCTTTCAGCAAACAGTTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACATTTc  <  1:2168035/68‑1 (MQ=255)
                                   cTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGc                   <  1:2580323/49‑1 (MQ=255)
                                       cAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCat                        >  1:617174/1‑40 (MQ=255)
                                       cAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCa                         >  1:2473660/1‑39 (MQ=255)
                                                 |                                                   
ACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTC  >  minE/1154233‑1154333

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: