Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,820,162 C→A 16.8% N188K (AAC→AAA gdhA → glutamate dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,820,1620C→A16.8% 2.5 ‑4.4 ‑0.118564gdhAglutamate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (13/2):  total (13/5)
Fisher's exact test strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.19e-01

GGCTTTATGGCGGGGATGATGAAAAAG.CTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCT  >  REL606/1820127‑1820194
                                    |                                
ggttttATGGCGGGGATGATGAAAAAG.CTCTCcaac                                  <  1:5453775/33‑1 (MQ=255)
ggCTTTATGGCGGGGGTGGTGGAAAAA.ATCTCcaac                                  >  1:2372501/1‑36 (MQ=16)
ggCTTTATGGCGGGGATGATGAAAAAG.CTCTCcaac                                  <  1:898412/36‑1 (MQ=255)
 gCTTTATGGCGGGGGTGATGAAAAAA.CTCTCcaaca                                 >  1:5658638/1‑36 (MQ=37)
  cTTTATGGCGGGTATGATGAAAAAG.CTCTCcaacaa                                <  1:3391671/36‑1 (MQ=255)
   tttATGGCGGGGATGATGAAAAAG.CTCTCCAACAAt                               <  1:5070528/36‑1 (MQ=255)
       tGGCGGGGATGATGAAAAAG.CTCTCCAACAATACCg                           >  1:5360658/1‑36 (MQ=255)
       tGGCGGGGATGATGAAAAAG.CTCTCCAACAATACCg                           >  1:5195018/1‑36 (MQ=255)
       tGGCGGGGATGATGAAAAAG.CTATCCAACAAAAcag                           >  1:4922778/1‑34 (MQ=25)
        ggCGGGGATGATGAAAAAG.CTCTCCAACAATACCGc                          >  1:125712/1‑36 (MQ=255)
        ggCGGGGATGATGAAAAAG.CTCTCCAACAATACCGc                          >  1:173460/1‑36 (MQ=255)
        ggCGGGGATGATGAAAAAG.CTCTCCAACAACACCGc                          >  1:5609976/1‑36 (MQ=255)
                  atgaAAAAG.CTCTCCAACAATACCGCCTGCGTCTTc                >  1:2075403/1‑36 (MQ=255)
                  atgaAAAAG.CTCTCCAACAATACCGCCTGCGTCTTc                >  1:2853203/1‑36 (MQ=255)
                  atgaAAAAG.CTCTCCAACAATACCGCCTGCGTCTTc                >  1:3517876/1‑36 (MQ=255)
                   tgaAAAAG.CTCTCCAACAATACCGCCTGCGTCTTCa               >  1:4197310/1‑36 (MQ=255)
                   gaaaaaaC.CTCCCAAAAAATCCCGCCTGCGTCTTCa               <  1:1047127/34‑1 (MQ=12)
                   agaAAAAGTCTCACCAAAAATCCCGCCTGCGTCTTc                <  1:5759436/35‑1 (MQ=12)
                     aaaaaC.CTCCCAAAAAAAACCGCCTGCGTCTTCAcc             <  1:1938175/36‑1 (MQ=12)
                      aaaaG.CTCTCCAACAATACCGCCTGCGTCTTCACCg            >  1:31295/1‑36 (MQ=255)
                       aaaG.CTCTCAAACAATACCGCCTGCGTCTTCACCgg           <  1:484475/36‑1 (MQ=255)
                        aaG.CTCTCCAACAATACCGCCTGCGTCTTCACCGGt          >  1:4896011/1‑36 (MQ=255)
                             tctcCAACAATACCGCCTGCGTCTTCACCGGTAAgg      >  1:4386413/1‑36 (MQ=255)
                                 caacaaTACCGCCTGCGTCTTCACCGGTAAGGGCCt  >  1:3468068/1‑36 (MQ=255)
                                    |                                
GGCTTTATGGCGGGGATGATGAAAAAG.CTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCT  >  REL606/1820127‑1820194

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG