Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 139,942 G→T 18.7% L19F (TTG→TTT cueO → multicopper oxidase (laccase)

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606139,9420G→T18.7% 2.2 ‑4.1 0.04257cueOmulticopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (18/13):  total (28/14)
Fisher's exact test strand distribution p-value = 1.72e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAACCGCGCAGTATTTGCGGCAGAACG  >  REL606/139907‑139977
                                   |                                   
aaaTATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTg                                     <  1:5306207/36‑1 (MQ=255)
  atatTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGcc                                   <  1:987685/36‑1 (MQ=255)
   tatTCCGTCGCGCTGGGTGTGGCTTCGGCTTTTCCg                                  >  1:4653657/1‑36 (MQ=255)
    atTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGc                                 <  1:3020137/36‑1 (MQ=255)
      tCCGTCGCGCTGGGTGTGGCTTCGGCTTTTCCGCtg                               >  1:7209102/1‑36 (MQ=255)
      tCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtg                               <  1:4487120/36‑1 (MQ=255)
      tCCGTCGCGCTGGGTGCGGCTTCGGCTTTGCCCCtg                               >  1:3392797/1‑36 (MQ=37)
        cGTCGCGCTGGGTGTGGCTTCGGCTTTTCCGCtgtg                             >  1:3858459/1‑36 (MQ=255)
        cGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                             <  1:4264724/36‑1 (MQ=255)
         gTCGCGCTGGGTGTGGCTTTGGGTTTGCCGCTGTgg                            >  1:328439/1‑36 (MQ=37)
         gTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg                            <  1:6818598/36‑1 (MQ=255)
         gTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg                            >  1:4721374/1‑36 (MQ=255)
         gTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg                            <  1:1384577/36‑1 (MQ=255)
         gTCGCGCTGGGTGTCGCTTCCGGTTTTCTCCTGTgg                            >  1:592865/1‑36 (MQ=2)
         gTCGCGCTGGGTGAGGCTTCGGCTTTTCCGCTGTgg                            >  1:2381619/1‑36 (MQ=37)
           cgcgCTGGGTGTGGCTTCGGCTTTTCCGCTGTTGaa                          >  1:335681/1‑36 (MQ=37)
           cgcgCTGGGTGTGGCTTCGGCTTTGCCGCTTTGGaa                          <  1:455546/36‑1 (MQ=255)
            gcgcTGGGTGTGGCTTCGGCTTTTCCGCTGTGGAAc                         >  1:1348943/1‑36 (MQ=255)
            gcgcTGGGTGTGGCTTCGGCTTTTCCCCTGTgggaa                         >  1:4642754/1‑33 (MQ=25)
            gcgcTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAAc                         <  1:4541675/36‑1 (MQ=255)
            gcgcTGGGTGTGGCTTCGGCTTTGCACCTGTgggcc                         >  1:4967269/1‑33 (MQ=25)
             cgcTGGGTGTGGCTTCGGCTTTGCCGCTGTTGAAcc                        >  1:1457271/1‑36 (MQ=255)
             cgcTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAAcc                        >  1:4093817/1‑36 (MQ=255)
             cgcTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAAcc                        <  1:1854177/36‑1 (MQ=255)
             cgcTGGGTGTGGCTTCGGCTTTACCGCTGTGGAAcc                        <  1:3214257/36‑1 (MQ=255)
              gcTGGGTGTGGCTTCGGCTTTTCCGCTGTGGAACcg                       >  1:4799730/1‑36 (MQ=255)
              gcTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAAccc                       >  1:3464644/1‑35 (MQ=255)
                 gggTGTGGCTTCGGCTTTGCCGCTGTGGAACCgcgc                    <  1:3915137/36‑1 (MQ=255)
                  ggTGTGGCTTCGGCTTTGCCGCTGTGGAACCGCGCa                   >  1:6074858/1‑36 (MQ=255)
                    tgtgGCTTCGGCTTTTCCGCTTTGGAACCGCGCCGt                 >  1:5376462/1‑36 (MQ=25)
                     gtgGCTTCGGCTTTGCCGATGTGGAACCGCGGAGTa                >  1:6638453/1‑36 (MQ=37)
                       ggCTTCGGCTTTGCCGCTGTGGAACCGCGCAGTAtt              >  1:6866822/1‑36 (MQ=255)
                        gCTTCGGCTTTGCCGCTGTGGAACCGCGCAGTAttt             >  1:4887585/1‑36 (MQ=255)
                        gCTTCGGCTTTGCCGCTGTGAAGCCGCGCAATAttt             <  1:1708787/36‑1 (MQ=21)
                        gCTTCGGCTTTGCAGCTGTGGAACCGCGCAGTAttt             >  1:3231158/1‑36 (MQ=255)
                          ttCGGCTTTGCCGCTGTGGAACCGCGCAGTATTTGc           >  1:5452111/1‑36 (MQ=255)
                          ttCGGCTTTGCCGCTGTGGAACCGCGCAGTATTTGc           >  1:7491439/1‑36 (MQ=255)
                            cGGCTTTGCCGCTGTGGAACCGCGCAGTATTTGCgg         >  1:7052446/1‑36 (MQ=255)
                            cGGCTTTGCCGCTGTGGAACCGCGCAGTATTTGCgg         >  1:7277048/1‑36 (MQ=255)
                            cGGCTTTGCCGCTGTGGAACCGCGCAGTATTTGCgg         >  1:822030/1‑36 (MQ=255)
                              gCTTTGCCGCTGTGGAACCGCGCAGTATTTGCGGCa       <  1:4913357/36‑1 (MQ=255)
                                  tGCCGCTGTGGAACCGCGCAGTATTTGCGGCAGAAc   <  1:2439373/36‑1 (MQ=255)
                                  tGCCGCTGTGGAACCGCGCAGTATTTGCGGCAGAAc   >  1:2744957/1‑36 (MQ=255)
                                   gCCGCTGTGGAACCGCGCAGTATTTGCGGCAGAACg  <  1:3304442/36‑1 (MQ=255)
                                   gCCGCTGTGGAACCGCGCAGTATTTGCGGCAGAACg  >  1:3155776/1‑36 (MQ=255)
                                   |                                   
AAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAACCGCGCAGTATTTGCGGCAGAACG  >  REL606/139907‑139977

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