Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,282,469 A→G 24.7% E819G (GAA→GGA)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,282,4690A→G24.7% 5.6 ‑3.5 0.0322456narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (13/10):  total (22/10)
Fisher's exact test strand distribution p-value = 3.17e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCAT  >  REL606/1282436‑1282502
                                 |                                 
aCCCGCTGTCTGCGGCGGTTGCTCCGGCATGGGGAg                                 >  1:6608270/1‑36 (MQ=12)
   cGCTGTCTGCGGCGGGCGATCCGGCCTGGGGAGCgg                              >  1:6817858/1‑35 (MQ=18)
    gCTGTCTGCGGCGGTTGCGGCCGCCTGGGAAGCgga                             >  1:4484601/1‑34 (MQ=9)
    gCTGTCTGCGGCGGTCGATCCGGCCTTGGGAGCgga                             >  1:4060084/1‑34 (MQ=18)
    gCTGTCTGCGGCGGTCCCTCCCGACTGGGGAGTGaa                             >  1:6167495/1‑36 (MQ=9)
     cTGTCTGCGGCGGTTGATCCCGCCTTGGGAGCGaaa                            >  1:3693892/1‑36 (MQ=19)
      tGTCTGCGGCGGTCCATCCGGCCTTGGGAGCGgaaa                           >  1:1529600/1‑36 (MQ=16)
      tGTCTGCGGCGGCCGCTCCGGCCTTGGGAGCCaaaa                           >  1:6748336/1‑36 (MQ=16)
       gTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAg                          <  1:6427203/36‑1 (MQ=255)
        tCTGCGGCGGTCGGTCCCGCCTTGGGAGCGgaaaag                         >  1:5704422/1‑34 (MQ=12)
           gcggcgGTCGATCCGGCCTGGGAAGCGAAAAGCGAc                      <  1:1208333/36‑1 (MQ=255)
            cggcggTCGATCCGGCCTGGGAAGCGAAAAGCGACt                     <  1:2274064/36‑1 (MQ=255)
             ggcggTCGATCCGGCCTGGGTAGCGAAAAGCGATTg                    >  1:5040115/1‑36 (MQ=25)
             ggcggTCGATCCGGCCTGGGAAGCGAAAAGCGACTg                    >  1:1219829/1‑36 (MQ=255)
              gcggTCGATCCGGCCTGGGAAGCGAAAAGCGACTgg                   >  1:5684605/1‑36 (MQ=255)
               cggTCGATCCGGCCTGGGAAGCGAAAAGCGACTggg                  <  1:6817461/36‑1 (MQ=255)
                 gTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaa                >  1:5766660/1‑36 (MQ=255)
                    gATCCGGCATGGGAAGCGAAAAGCGACTGGGAAATc             <  1:2942200/36‑1 (MQ=255)
                     aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt            >  1:1467823/1‑36 (MQ=255)
                      tCCGGCCTGGTAAGCGAAAAGCGACTGGGAAATCTa           <  1:7185987/36‑1 (MQ=255)
                        cGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACa         >  1:2553940/1‑36 (MQ=255)
                         ggCCTGGGAAGCGAAAAGCGACTGGGAAATCTACaa        >  1:2441156/1‑36 (MQ=255)
                         ggCCTGGGAAGCGAAAAGCGACTGGGAAATCTACaa        <  1:119716/36‑1 (MQ=255)
                            cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc     >  1:7208024/1‑36 (MQ=255)
                            cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc     <  1:5505508/36‑1 (MQ=255)
                            cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc     <  1:4089112/36‑1 (MQ=255)
                            cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc     >  1:4024677/1‑36 (MQ=255)
                            cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc     >  1:1456510/1‑36 (MQ=255)
                             tGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGcc    >  1:5705495/1‑36 (MQ=255)
                             tGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGcc    >  1:4967283/1‑36 (MQ=255)
                             tGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGcc    >  1:4649146/1‑36 (MQ=255)
                               ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCAt  <  1:271707/36‑1 (MQ=255)
                                 |                                 
ACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCAT  >  REL606/1282436‑1282502

Base quality scores:  ATCG  < 3 ≤  ATCG  < 4 ≤  ATCG  < 12 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG