Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,609,783 C→T 34.9% Q191Q (CAG→CAA eutC ← ethanolamine ammonia‑lyase, small subunit (light chain)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,609,7830CT34.9% 62.5 / 42.8 80Q191Q (CAG→CAAeutCethanolamine ammonia‑lyase, small subunit (light chain)
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base T (3/25);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 1.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

CGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGATCG  >  minE/1609716‑1609838
                                                                   |                                                       
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:239113/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:820452/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:1544919/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:1273264/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:1718466/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:2964260/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:882118/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:1362597/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:1492850/69‑1 (MQ=255)
cGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:1463393/69‑1 (MQ=255)
 gatctgatctTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:396515/68‑1 (MQ=255)
       atcttcaatcttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAg                                                 >  1:2408053/1‑69 (MQ=255)
       atcttcaatcttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAg                                                 >  1:510371/1‑69 (MQ=255)
       atcttcaatcttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAg                                                 >  1:40782/1‑69 (MQ=255)
         cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                                                      <  1:1705113/62‑1 (MQ=255)
         cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                                                      <  1:2610546/62‑1 (MQ=255)
         cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                                                      <  1:3065587/62‑1 (MQ=255)
         cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                                                      <  1:374880/62‑1 (MQ=255)
         cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                                                      <  1:2532493/62‑1 (MQ=255)
         cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAgg                                                >  1:919867/1‑68 (MQ=255)
         cttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGc                                               >  1:3042205/1‑69 (MQ=255)
          ttcaatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGt                                                      <  1:2937952/61‑1 (MQ=255)
            caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTc                                                   >  1:2121026/1‑62 (MQ=255)
            caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGccc                                             >  1:1259770/1‑68 (MQ=255)
            caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGccc                                             >  1:628993/1‑68 (MQ=255)
            caatcttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGccc                                             >  1:2182979/1‑68 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:2785520/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:3029553/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:366460/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:3122675/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:3244234/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:2274686/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:1001378/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:2094780/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:1297474/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:1559236/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:1175814/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:681650/69‑1 (MQ=255)
                cttcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:657807/69‑1 (MQ=255)
                cttcaCGCGACCATAACATACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:2685841/69‑1 (MQ=255)
                 ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc                                       >  1:680004/1‑69 (MQ=255)
                 ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc                                       >  1:1030650/1‑69 (MQ=255)
                 ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc                                       >  1:1926970/1‑69 (MQ=255)
                 ttcaCGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATc                                       >  1:2290453/1‑69 (MQ=255)
                  tcaCGCGACCATAACGTACAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCAt                                        <  1:3042524/67‑1 (MQ=255)
                           cATAACGAACAAAGAATGGTGTACCACTTTCAGCCCGGCCTGt                                                      <  1:1458480/43‑1 (MQ=38)
                             tAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCt                                                        <  1:1638221/40‑1 (MQ=255)
                               aCGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGt                                 >  1:834728/1‑61 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:3287852/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:3238224/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:1238814/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:3408091/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:1253697/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:927461/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:2856868/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:2747526/60‑1 (MQ=255)
                                cgtaCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:177821/60‑1 (MQ=255)
                                  taCAAAGAATGGCGTACCGACTTTCAGCCCGGCTTGCTTCAGGCCCGCCATCAGCggcggc                              <  1:249569/60‑1 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTTGCAGGATCTCTTCGt                 >  1:120797/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:430095/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:1022606/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:562216/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:1361907/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:2930773/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:2866608/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:2792325/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:1979504/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:1917726/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:1716998/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:155805/1‑60 (MQ=255)
                                                tACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:1528359/1‑60 (MQ=255)
                                                  ccGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGt       <  1:1352729/68‑1 (MQ=255)
                                                  ccGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGt       <  1:2767967/68‑1 (MQ=255)
                                                  ccGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGt       <  1:1685902/68‑1 (MQ=255)
                                                  ccGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTg      >  1:1894031/1‑69 (MQ=255)
                                                    gACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGt                 >  1:1778973/1‑56 (MQ=255)
                                                    gACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGAt    >  1:941801/1‑69 (MQ=255)
                                                    gACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGAt    >  1:3101441/1‑69 (MQ=255)
                                                      cTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGATCg  <  1:3103229/69‑1 (MQ=255)
                                                      cTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGATCg  <  1:3292912/69‑1 (MQ=255)
                                                                   |                                                       
CGATCTGATCTTCAATCTTCACGCGACCATAACGAACAAAGAATGGTGTACCGACTTTCAGCCCGGCCTGTTTCAGGCCCGCCATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGATCG  >  minE/1609716‑1609838

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: