Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,990,234 A→G 70.5% L76L (TTA→TTG mutY → adenine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,990,2340AG70.5% 57.2 / 27.8 44L76L (TTA→TTGmutYadenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base A (8/5);  new base G (13/18);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTA  >  Exported/1990167‑1990303
                                                                   |                                                                     
gCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTgg                                                                    <  1:3250390/71‑1 (MQ=255)
gCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTgg                                                                    <  1:506624/71‑1 (MQ=255)
      tGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCgg                                                               <  1:1121892/70‑1 (MQ=255)
      tGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCgg                                                               <  1:1853017/70‑1 (MQ=255)
            gcTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGc                                                        <  1:2316699/71‑1 (MQ=255)
             cTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCt                                                       <  1:80387/71‑1 (MQ=255)
             cTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCt                                                       <  1:4917634/71‑1 (MQ=255)
               tCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTa                                                      >  1:5110098/1‑70 (MQ=255)
                cccGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTAt                                                     <  1:2926930/70‑1 (MQ=255)
                 ccGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGa                                                                   >  1:1016617/1‑55 (MQ=255)
                   gACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACg                                                 >  1:4917724/1‑71 (MQ=255)
                      ggTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGccc                                              >  1:5493151/1‑71 (MQ=255)
                       gTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTtgtg                                                                     >  1:2470683/1‑47 (MQ=255)
                        tGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCcgc                                            <  1:104304/71‑1 (MQ=255)
                        tGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTAt                                                     <  1:1726608/62‑1 (MQ=255)
                         gACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCcgcg                                           >  1:4425030/1‑71 (MQ=255)
                          aCCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCcgcgc                                          >  1:198280/1‑71 (MQ=255)
                               tctcGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCgcgcgcgc                                      >  1:3293503/1‑70 (MQ=255)
                                ctcGCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCaa                                    >  1:1715202/1‑71 (MQ=255)
                                   gCCAATGCGCCGCTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCa                                     <  1:4852843/67‑1 (MQ=255)
                                   gCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCt                                 <  1:5208161/71‑1 (MQ=255)
                                   gCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCt                                 >  1:414096/1‑71 (MQ=255)
                                       aTGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCa                              >  1:2629534/1‑70 (MQ=255)
                                          cgccgcTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATaa                           <  1:2919316/70‑1 (MQ=255)
                                          cgccgcTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATaa                           <  1:4252385/70‑1 (MQ=255)
                                           gccgcTCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTAt                                                     <  1:941847/43‑1 (MQ=255)
                                             cgcTCGACGAAGTTCTCCACTTGTGGACCGGGCTTg                                                          <  1:363272/36‑1 (MQ=255)
                                             cgcTCGACGAAGTTCTCCACTTGTGGACCGGGCTTg                                                          <  1:3468342/36‑1 (MQ=255)
                                             cgcTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCg                       >  1:2549173/1‑71 (MQ=255)
                                                tCGACGAAGTTCTCCACTTGTGGACCGGGCTTGGCTATTa                                                   >  1:244156/1‑40 (MQ=255)
                                                 cgacgaAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGccc                                              >  1:5244287/1‑44 (MQ=255)
                                                     gaAGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCacaaca               <  1:3661540/71‑1 (MQ=255)
                                                     gaAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGccc                                              >  1:1412877/1‑40 (MQ=255)
                                                     gaAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCcg                                             >  1:3528880/1‑41 (MQ=255)
                                                     gaAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCcg                                             >  1:2628010/1‑41 (MQ=255)
                                                      aaGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCacaac                >  1:4749677/1‑69 (MQ=255)
                                                       aGTTCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAg             <  1:2155732/71‑1 (MQ=255)
                                                        gTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGt            <  1:1103200/71‑1 (MQ=255)
                                                         ttCTCCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTg           <  1:4189981/71‑1 (MQ=255)
                                                           ctcCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGc         >  1:476069/1‑71 (MQ=255)
                                                           ctcCACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGc         >  1:5511793/1‑71 (MQ=255)
                                                           ctcCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCacaac                <  1:5696520/64‑1 (MQ=255)
                                                             ccACTTGTGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGa       >  1:3146200/1‑71 (MQ=255)
                                                              cACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGc         <  1:4405362/68‑1 (MQ=255)
                                                                  tATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTa  >  1:1374745/1‑71 (MQ=255)
                                                                   |                                                                     
GCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCCACTTATGGACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTA  >  Exported/1990167‑1990303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: