Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 510,400 G→T 12.8% W201L (TGG→TTG)  allP → predicted allantoin transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606510,4000G→T12.8% 3.0 ‑4.6 ‑0.154602allPpredicted allantoin transporter
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (21/26):  total (28/26)
Fisher's exact test strand distribution p-value = 1.04e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.63e-01

GTGCATCTATATCGTTTTCGGCGGTATGGC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCTT  >  REL606/510365‑510433
                                    |                                 
gTGCATCTATATCGTTTTCGGCGGTATGGC.GATTTg                                   <  1:3375409/36‑1 (MQ=255)
  gCATCTATATCGTTTTCGGCGGTATGGC.GATTTggg                                 <  1:2321175/36‑1 (MQ=255)
    aTCTATATCGTTTTCGGCGGTTTGGC.GGTTTtgggg                               >  1:1120220/1‑34 (MQ=21)
    aTCTATATCGTTTTCGGCGGTATTGC.GATTTtgggg                               >  1:4696077/1‑34 (MQ=25)
    aTCTATATCGTTTTCGGCGGTATGGC.GATTTTGGCg                               >  1:2348588/1‑36 (MQ=255)
      cTATATCGTTTTCGGCGGTATGGC.GATTTGGGCGAt                             <  1:3126310/36‑1 (MQ=255)
        atatCGTTTTCGGCTGTATGGC.GATTTGGGCGAttt                           <  1:4442712/36‑1 (MQ=255)
         tatCGTTTTCGGCGGTATGGC.GATTTGGGCGATTTc                          <  1:5830120/36‑1 (MQ=255)
          atCGTTTTCGGCGGTTTGGCGGATTTGGGGGTtttt                          >  1:2789403/1‑35 (MQ=12)
          atCGTTTTCGGCGGTATTGC.GATTTTGGTGGTTTCa                         >  1:2010922/1‑36 (MQ=16)
          atCGTTTTCGGCGGTATGGC.GATTTGGGCGATTTCa                         <  1:1089512/36‑1 (MQ=255)
          atCGTTTTCGGCGGTATCGC.GATTTTGGGGATTTcc                         >  1:6694119/1‑35 (MQ=21)
           tCGTTTTCGGCGGTATGGC.GATTTTGTCGATTTTAc                        >  1:536798/1‑36 (MQ=25)
           tCGTTTTCGGCGGTATGGC.GATTTTGGCGGttttat                        >  1:5692710/1‑33 (MQ=25)
           acGTTTTCGTGTGTATGAC.GATTTGGGCTATTTCAc                        <  1:1934423/35‑1 (MQ=12)
            cGTTTTCGGCGGTATGGC.GATTTGGGCGATTTCACt                       <  1:4940808/36‑1 (MQ=255)
              ttttCGGCGGTATGGC.GATTTGGGCGATTTCACTgg                     <  1:40458/36‑1 (MQ=255)
              tttgCGGCGGTATGGC.GATTTGGGCGATTTCACTgg                     <  1:3608752/36‑1 (MQ=255)
               tttCGGCGGTATTGC.GATTTTGGGGGTTTTACTGGt                    >  1:6895082/1‑36 (MQ=12)
               tttCGGCGGTATGGC.GATTTGGGCGATTTCACTGGt                    <  1:2063237/36‑1 (MQ=255)
                ttCGGCGGTATGGC.GATTTTGGGGATTTTACTTGtt                   >  1:3162512/1‑35 (MQ=16)
                ttCGGCGGTATGGC.GATTTGGGCGATTTCGCTGGTc                   <  1:375217/36‑1 (MQ=255)
                ttCGGCGGTATGGC.GATTTGGGCGATTTCACTGGTc                   <  1:3280160/36‑1 (MQ=255)
                   ggcggTGTGGC.GATTTGGGCGANTTCACTGGTCggg                <  1:4342742/36‑1 (MQ=38)
                   ggcggTATGGC.GATTTGGGCGATTTTACTGGTAggg                >  1:1929524/1‑36 (MQ=37)
                   ggcggTATGGC.GATTTGGGCGATTTCACTGGTCggg                <  1:4739560/36‑1 (MQ=255)
                   ggcggTATGGC.GATTTGGGCGATTTCACTGGTCggg                >  1:763939/1‑36 (MQ=255)
                    gcggTATGGC.GATTTGGGCGATTTCACTGGTCGGGa               <  1:820583/36‑1 (MQ=255)
                     cggTATGGC.GATTTGGGCGATTTCACTGGTCGGGAt              >  1:4827953/1‑36 (MQ=255)
                      ggTATGGC.GATTTGGGCGATTTCACTGGTCGGGATc             <  1:2141325/36‑1 (MQ=255)
                      ggTATGGC.GATTTGGGCGATTTCACTGGTCGGGATc             >  1:6232035/1‑36 (MQ=255)
                       gTATGGC.GTCTTGGGCNATTTCTCTTGTCGGGATCg            <  1:2696483/36‑1 (MQ=255)
                        tATGGC.GATTTGGGCGATTTCACTGGTCAGGATCgg           >  1:1876475/1‑36 (MQ=255)
                         aTGGC.GATTTGGGCGATTTCACTGGTCGGGATCGGt          <  1:6305515/36‑1 (MQ=255)
                          tGGC.GATTTGGGCGATTTCACTGGTCGGGATCGGTc         <  1:1521480/36‑1 (MQ=255)
                          tGGC.GATTTGGGCGATTTCACTGGTCGGGATCGGTc         <  1:4909264/36‑1 (MQ=255)
                           ggC.GATTTGGGCGATTTCACTGGTCGGGATCGGTcc        >  1:2443773/1‑36 (MQ=255)
                           ggC.GATTTGGGCGATTTCACTGGTCGGGATCGGTcc        >  1:6400377/1‑36 (MQ=255)
                            gC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCa       >  1:1073729/1‑36 (MQ=255)
                            gC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCa       >  1:793376/1‑36 (MQ=255)
                            gC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCa       >  1:7322291/1‑36 (MQ=255)
                            gC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCa       >  1:6630171/1‑36 (MQ=255)
                            gC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCa       >  1:6405481/1‑36 (MQ=255)
                            gC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCa       <  1:1669395/36‑1 (MQ=255)
                            gC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCa       <  1:3855205/36‑1 (MQ=255)
                             c.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCaa      <  1:1813815/36‑1 (MQ=255)
                             c.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCaa      <  1:4332481/36‑1 (MQ=255)
                               gATTTGGGCGATTTCACTGGTCGGGATCGGTCCAAt     <  1:5661977/36‑1 (MQ=255)
                               gATTTGGGCGATTTCACTGGTCGGGATCGGTCCAAt     >  1:1554218/1‑36 (MQ=255)
                                 tttGGGCGATTTCACTGGTCGGGATCGGTTCAATCt   >  1:263241/1‑36 (MQ=255)
                                 tttGGGCGATTTCACTGGTCGGGATCGGTCCAATCt   >  1:592454/1‑36 (MQ=255)
                                 tttGGGCGATTTCACTGGTCGGGATCGGTCCAATCt   <  1:5348456/36‑1 (MQ=255)
                                 tttGGGCGATTTCACTGGTCGGGATCGGTCCAATCt   >  1:6488754/1‑36 (MQ=255)
                                 tttGGGCGATTTCACTGGTCGGGATCGGTCCAATCt   >  1:3380246/1‑36 (MQ=255)
                                  ttGGGCGATTTCACTGGTCGTGATCGGTCCAATCtt  >  1:1808841/1‑36 (MQ=255)
                                  ttGGGCGATTTCACTGGTCGGGATCGGTCCAATCtt  >  1:5221788/1‑36 (MQ=255)
                                  ttGGGCGATTTCACTGGTCGGGATCGGTCCAATCtt  <  1:423376/36‑1 (MQ=255)
                                  ttGGGCGATTTCACTGGTCGGGATCGGTCCAATCtt  >  1:1506280/1‑36 (MQ=255)
                                    |                                 
GTGCATCTATATCGTTTTCGGCGGTATGGC.GATTTGGGCGATTTCACTGGTCGGGATCGGTCCAATCTT  >  REL606/510365‑510433

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