Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,259,415 A→C 17.5% W11G (TGG→GGG)  alkB ← oxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,259,4150A→C17.5% 3.9 ‑3.3 ‑0.04231alkBoxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (20/14):  total (28/14)
Fisher's exact test strand distribution p-value = 3.68e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

AAATTACCGCTCCAGCCGCCAGTGGCTCTTGCCACGGTTCACCATCGGCAAACAGATCCAACATTACCT  >  REL606/2259382‑2259450
                                 |                                   
aaaTTACCGCTCCAGCCGCCAGTGGCTCTTGCCACg                                   <  1:1016818/36‑1 (MQ=255)
aaaTTACCGCTCCAGCCGCCAGTGGCTCTTGCCACg                                   >  1:2716187/1‑36 (MQ=255)
 aaTTCCCTCTCAAGCCGCCAGTGGCTCTTGCCACgg                                  <  1:7132790/36‑1 (MQ=25)
 aaTTACCGCTCCAGCCGCCAGTGGCTCTTTCCCCgg                                  >  1:5812974/1‑36 (MQ=37)
   ttcccGCCCCAGCCGCCAGTGGCTCTTGCCACGGtt                                <  1:4685562/33‑1 (MQ=37)
   ttACCGCTCCAGCCGCCAGTGGCTCTTTCCACGGtt                                >  1:7486460/1‑36 (MQ=255)
    tACCGCTCCAGCCGCCAGTGGCTCTTTCCCCGGTTc                               >  1:2535325/1‑36 (MQ=37)
     cccGCTCCAGCCGCCAGTGGCTCTTGCCACGGTTCa                              <  1:5077996/35‑1 (MQ=255)
       cGCTCCAGCCGCCAGTGGCTCTTGCCACGGTTCAcc                            >  1:4011981/1‑36 (MQ=255)
        gCTCCAGCCGCCAGTGGCTCTTTCCCCGGTTTACCa                           >  1:5914554/1‑36 (MQ=21)
        gCTCCAGCCGCCAGTGGCTCTTTCCCCGGGTCACCa                           >  1:6934258/1‑36 (MQ=21)
        gCTCCAGCCGCCAGTGGCTCTTGCCACGGTTCACCa                           >  1:4635231/1‑36 (MQ=255)
        gCTCCAGCCGCCAGTGGCTCCTGCCCCGGTTCCccc                           >  1:3360645/1‑35 (MQ=21)
         ccccAGCCGCCAGTGGCTCTTGCCACGGTTCACCAt                          <  1:3880925/34‑1 (MQ=255)
         cTCCAGCCGCCAGTGGCTCCTTCCACGGTTCACCAt                          >  1:7121401/1‑36 (MQ=37)
          tCCAGCCGCCAGTGGCTCTTGCCCCGGTTCCCCATc                         >  1:2466581/1‑36 (MQ=37)
          cccAGCCGCCAGTGGCTCTTGCCACGGTTCACCATc                         <  1:4954000/35‑1 (MQ=255)
          cccAGCCGCCAGTGGCTCTTGCCACGGTTCACCATc                         <  1:4779583/35‑1 (MQ=255)
           ccAGCCGCCAGTGGCTCTTGCCACGGTTCACCATCg                        <  1:4666278/36‑1 (MQ=255)
           ccAGCCGCCAGTGGCTCTTACCACGGTTCACCATCg                        >  1:2183084/1‑36 (MQ=255)
              gccgccAGTGGCTCTTGCCCCGGTTCACCATCGGCa                     >  1:647863/1‑36 (MQ=255)
               ccgccAGTGGCTCCTGCCACCGTTCACCATCGGCaa                    >  1:1460853/1‑36 (MQ=25)
                cgccAGTGGCTCTTTCCACGGTTCACCATCGGCaaa                   >  1:3394238/1‑36 (MQ=255)
                  ccAGTGGCTCTTGCCCCCGTTCACCATCGGCAAACa                 >  1:2481950/1‑36 (MQ=37)
                  ccAGTGGCTCTTGCCACGGTTCACCATCGGCAAACa                 <  1:7339692/36‑1 (MQ=255)
                  ccAGTGGCTCTTGCCACGGTTCACCATCGGCAAACa                 <  1:1367337/36‑1 (MQ=255)
                   cAGTGGCTCTTGCCACGGTTCACCATCGGCAAACAg                <  1:1295934/36‑1 (MQ=255)
                      tGGCTCTTGCCACGGTTCACCATCGGCAAACAGATc             <  1:5503876/36‑1 (MQ=255)
                        gCTCTTGCCACGGTTCACCATCGGCAAACAGATCCa           >  1:3899254/1‑36 (MQ=255)
                        gCTCTTGCCACGGTTCACCATCGCCAAATAGATCCa           >  1:51155/1‑36 (MQ=25)
                         ctctTGCCACGGTTCACCATCGGCAAACAGATCCaa          >  1:4094248/1‑36 (MQ=255)
                         ctctTGCCACGGTTCACCATCGGCAAACAGATCCaa          >  1:1728344/1‑36 (MQ=255)
                          tctTGCCACGGTTCACCATCGGCAAACAGATCCAAc         >  1:5550403/1‑36 (MQ=255)
                           ctTGCCACGGTTCACCATCGGCAAACAGATCCAACa        >  1:6024878/1‑36 (MQ=255)
                            ttGCCACGGTTCACCATCGGCAAACAGATCCAACAt       >  1:5448529/1‑36 (MQ=255)
                            ttGCCACGGTTCACCATCGGCAAACAGATCCAACAt       <  1:7163880/36‑1 (MQ=255)
                              gCCACGGTTCACCATCGGCAAACAGATCCAACATTa     <  1:5692384/36‑1 (MQ=255)
                              gCCACGGTTCACCATCGGCAAACAGATCCAACATTa     >  1:5362466/1‑36 (MQ=255)
                              gCCACGGTTCACCATCGGCAAACAGATCCAACATTa     >  1:1078245/1‑36 (MQ=255)
                               ccACGGTTCACCATCGGCAAACAGATCCAACATTAc    >  1:6200315/1‑36 (MQ=255)
                                cACGGTTCACCATCGGCAAACAGATCCAACATTAcc   >  1:3852159/1‑36 (MQ=255)
                                cACGGTTCACCATCGGCAAACAGATCCAACATTAcc   >  1:4579317/1‑36 (MQ=255)
                                 aCGGTTCACCATCGGCAAACAGATCCAACATTACct  >  1:1334309/1‑36 (MQ=255)
                                 aCGGTTCACCATCGGCAAACAGATCCAACATTACct  <  1:90172/36‑1 (MQ=255)
                                 |                                   
AAATTACCGCTCCAGCCGCCAGTGGCTCTTGCCACGGTTCACCATCGGCAAACAGATCCAACATTACCT  >  REL606/2259382‑2259450

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