Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,294,000 G→T 25.8% intergenic (+150/‑92) nrdA → / → nrdB ribonucleotide‑diphosphate reductase alpha subunit/ribonucleotide‑diphosphate reductase beta subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,294,0000G→T25.8% 6.4 ‑3.8 ‑0.036intergenic (+150/‑92)nrdA/nrdBribonucleotide‑diphosphate reductase alpha subunit/ribonucleotide‑diphosphate reductase beta subunit
Reads supporting (aligned to +/- strand):  new base (9/1):  ref base (11/15):  total (20/16)
Fisher's exact test strand distribution p-value = 2.19e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

GCCTTATCCGGCCTACGGCTCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATC  >  REL606/2293965‑2294031
                                   |                               
gCCTTATCCGGCCTACGGCTCTACGGCTCGGATTTg                                 <  1:2592420/36‑1 (MQ=255)
 ccTTATCCGGCCTACGGCTCTACGGCTCGGATTTGt                                <  1:2726074/36‑1 (MQ=255)
  cTTATCCGGCCTACGGCTCTACGGCTCGGGTTTTTa                               >  1:3183833/1‑36 (MQ=21)
    tATCCGGCCTACGGCTCTACGGCTCGGATTTGTAgg                             >  1:296342/1‑36 (MQ=255)
     aTCCGGCCTACGGCTCTACGGCTCGGATTTGTAGGc                            <  1:8268237/36‑1 (MQ=255)
        cGGCCTACGGCTCTACGGCTCGGATTTGTAGGCCTg                         >  1:6876630/1‑36 (MQ=255)
         ggCCTACGGCTCTACGGCTCGGATTTTTAGGCCTgt                        >  1:5290290/1‑35 (MQ=255)
           ccTACGGCTCTACGTCTCGGATTTGTAGGCCTGATa                      <  1:3393321/36‑1 (MQ=255)
           ccTACGGCTCTACGGCTCGGATTTTTAGGCCTGATa                      >  1:7643406/1‑36 (MQ=255)
           ccTACGGCTCTACGGCTCGGATTTGTAGGCCTGATa                      <  1:8217633/36‑1 (MQ=255)
            ctacggctctacggctCGGGTTTTTAGGCCTGATaa                     >  1:7803609/1‑36 (MQ=21)
            ctacggctctacggctCGGGTTTTTAGGCCTGATaa                     >  1:4421614/1‑36 (MQ=21)
            ctacggctctacggctCGGGTTTTTAGGCCTGANna                     >  1:2339724/1‑36 (MQ=21)
            ctacggctctacggctCCGATTTTTAGGCCTTATaa                     >  1:5168467/1‑36 (MQ=37)
             tacggctctacggctcGGGTTTTTAGGCCTGATAAg                    >  1:5288293/1‑36 (MQ=21)
             tacggctctacggctcGGATTTGTAGGCCTGATAAg                    <  1:7311681/36‑1 (MQ=255)
                ggctctacggctcGGATTTGTAGGCCTGATAAGAcg                 <  1:5588670/36‑1 (MQ=255)
                 gctctacggctcGGATTTTTAGGCCTGATAAGAcgc                <  1:374261/36‑1 (MQ=255)
                 gctctacggctcGGATTTGTAGGCCTGATAAGAcgc                >  1:3110828/1‑36 (MQ=255)
                 gctctacggctcGGATTTGTAGGCCTGATAAGAcgc                <  1:5887232/36‑1 (MQ=255)
                  ctctactgctcGGATTTGTAGGCCTGATAAGAcgcg               <  1:6453372/36‑1 (MQ=255)
                  ctctacggctcGGCTTTGTAGGCCTGATAAGAcgcg               >  1:1981639/1‑36 (MQ=255)
                  ctctacggctcGGATTTGTAGGCCTTATAAGAcgcg               >  1:4092758/1‑36 (MQ=255)
                  ctctacggctcGGATTTGTAGGCCTGATAAGAcgcg               <  1:5065802/36‑1 (MQ=255)
                   tctacggctcGGATTTTTATGCCTGATAAGACgcgc              >  1:4183940/1‑36 (MQ=18)
                   tctacggctcGGATTTGTAGGCCTTATAAGACgcgc              >  1:4355195/1‑36 (MQ=255)
                    ctacggctcGGATTTGTAGGCCTGATAAGACGCGcc             <  1:8110509/36‑1 (MQ=255)
                       cggctcGGATTTGTAGGCCTGATAAGACGCGCCAGc          <  1:5391808/36‑1 (MQ=255)
                       cggctcGGATTTGTAGGCCTGATAAGACGCGCCAGc          <  1:7095006/36‑1 (MQ=255)
                         gctcGGATTTGTAGGCCTGATAAGACGCGCCAGCGt        <  1:5452769/36‑1 (MQ=255)
                          ctcGGATTTGTAGGCCTGATAAGACGCGCCAGCGTc       <  1:1531100/36‑1 (MQ=255)
                            cGGATTTGTAGGCNNGATNNGACGCGCCCGCGTCGc     >  1:3442269/1‑36 (MQ=14)
                             ggATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa    >  1:7424902/1‑36 (MQ=38)
                              gATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt   >  1:3397023/1‑36 (MQ=35)
                              gATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt   >  1:2640834/1‑36 (MQ=35)
                              gATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt   >  1:2256022/1‑36 (MQ=35)
                               aTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATc  <  1:7941624/36‑1 (MQ=32)
                                   |                               
GCCTTATCCGGCCTACGGCTCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATC  >  REL606/2293965‑2294031

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG