Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 481,510 +G 4.8% coding (2038/2505 nt) copA ← copper transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606481,5101.→G 4.8% 3.5 ‑2.0 ‑0.442coding (2038/2505 nt)copAcopper transporter
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (14/26):  total (16/26)
Fisher's exact test strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.23e-01

CTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGAC  >  REL606/481475‑481545
                                    |                                     
cTCATCAATCCCTGCTTCTTTGGCGATCGCATTTGc                                        >  1:641853/1‑36 (MQ=255)
 tcatcaATCCCTGCTTCTTTGGCGATCGCATTGGC.g                                      <  1:1305143/36‑1 (MQ=255)
  catcaATCCCTGCTTCTTTGGCGATCGCATTGGC.gg                                     >  1:3183630/1‑36 (MQ=255)
   atcaATCCCTGCTTCTTTGGCGATCGCATTGGC.ggt                                    <  1:6478393/36‑1 (MQ=255)
    tcaATCCCTGCTTCTTTGGCGATCGCATTGGC.ggt.g                                  >  1:5921143/1‑36 (MQ=255)
    tcaATCCCTGCTGCTTTGGCGATCGCATTGGC.ggt.g                                  <  1:4813010/36‑1 (MQ=255)
     caATCCCTGCTTCTTTGGCGATCGCATTGGC.ggt.gg                                 >  1:1814074/1‑36 (MQ=255)
      aaTCCCTGCTTCTTTGGCGATCGCATTGGC.ggt.gg.t                               >  1:7824557/1‑36 (MQ=255)
       aTCCCTGCTTCTTTGGCGATCGCATTGGC.ggg.gggt                               >  1:8076324/1‑34 (MQ=37)
       aTCCCTGCTTCTTTGGCGATCGCATTGGC.ggg.gg.tg                              >  1:5351349/1‑35 (MQ=255)
       aTCCCTGCTTCTTTGGCGATCGCATTGGC.GGT.GG.tt                              <  1:1414288/36‑1 (MQ=255)
       aTCCCTGCTTCTTTGGCGATCGCATTGGC.GGG.GG.tt                              >  1:1434791/1‑36 (MQ=255)
         cccTGCTTCTTTGGCTATCGCATTGGC.GGT.GG.TTgg                            <  1:8145352/36‑1 (MQ=255)
         cccTGCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTgg                            <  1:2063527/36‑1 (MQ=255)
         cccTGCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTgg                            <  1:2680595/36‑1 (MQ=255)
         cccTGCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTgg                            <  1:1484836/36‑1 (MQ=255)
          ccTGCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTtgg                           >  1:7044436/1‑36 (MQ=255)
             gCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTGGGTTa                        <  1:1473252/36‑1 (MQ=255)
             gCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTGGGTTa                        <  1:5655238/36‑1 (MQ=255)
               ttcttTGGCGATCGCATTGGC.GGT.GG.TTGGGTTATc                      <  1:1623776/36‑1 (MQ=255)
                tcttTGGCGATCGCATTGGC.GGT.GG.TTGGGTTATcc                     <  1:518359/36‑1 (MQ=255)
                  tttGGCGATCGCATTGGC.GGT.GG.TTGGGTTATcccc                   <  1:3833883/36‑1 (MQ=255)
                  tttGGCGATCGCATTGGC.GGT.GG.TTGGGTTATcccc                   <  1:7613152/36‑1 (MQ=255)
                   ttGGCGATCGCATTGGC.GGT.GG.TTGGGTTATCCCCg                  <  1:4469164/36‑1 (MQ=255)
                   ttGGCGATCGCATTGGC.GGT.GG.TTGGGTTATCCCCg                  <  1:6268853/36‑1 (MQ=255)
                   tgggCGATCGCATTTGC.GGT.GT.TTGGGTTATCCCCg                  <  1:1988431/34‑1 (MQ=37)
                      gCGATCGCATTGGCGGGG.GG.TTGTGTTTTCACCGGt                >  1:5469209/1‑36 (MQ=9)
                      gCGATCGCATTGGCGGGG.GG.TTGGTTTTTTCCCggc                >  1:5223669/1‑35 (MQ=9)
                      gCGATCGCATTGGC.GGT.GG.TTGGGTTATCCCCGGTc               <  1:4657352/36‑1 (MQ=255)
                      gCGATCGCATTGGC.GGG.GG.GTTGGTTTTTCCCGGtg               >  1:6842539/1‑35 (MQ=12)
                      gCGATCGCATTGGC.GGG.GG.GTGGGTTTTGCCCgggg               >  1:6143800/1‑34 (MQ=16)
                      gCGATCGAATTGGC.GGT.GG.TTGGGTTATCCCCGGTc               <  1:2472862/36‑1 (MQ=255)
                       cGATCGCATTGGC.GGT.GG.TTGGGTTATCCCCGGTCa              <  1:1717705/36‑1 (MQ=255)
                        gATCGCATTGGC.GGTGGG.TGGGGTTTTCCCCGGCCa              >  1:7670646/1‑36 (MQ=16)
                         aTCGGATTGGC.GGT.GG.TTGGGTTATCCCCGGTCAAc            <  1:2256248/36‑1 (MQ=255)
                         aTCGCATTGGC.GGC.GGTTTGGCTTTTCCCCGGTCac             >  1:6218754/1‑35 (MQ=12)
                           cGCATTGGC.GGT.TG.GTGTGGTATTCCCGGGCCACac          >  1:476174/1‑35 (MQ=2)
                            gCATTGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATc         <  1:6487283/36‑1 (MQ=255)
                             ccttGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCa        <  1:1218006/34‑1 (MQ=255)
                             ccttGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCa        <  1:1118706/34‑1 (MQ=255)
                             cATTGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCa        <  1:7809314/36‑1 (MQ=255)
                              aTTGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCAc       <  1:7687438/36‑1 (MQ=255)
                               ttGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCAcc      <  1:6071870/36‑1 (MQ=255)
                               ttGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCAcc      <  1:8423661/36‑1 (MQ=255)
                               ttGGC.GGT.GG.TTGGGTTATCCCCGGTCAACACCAcc      >  1:68420/1‑36 (MQ=255)
                                 ggC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAg    <  1:5684759/36‑1 (MQ=255)
                                 ggC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAg    >  1:5381840/1‑36 (MQ=255)
                                  gC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGa   >  1:4349109/1‑36 (MQ=255)
                                  gC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGa   >  1:4313158/1‑36 (MQ=255)
                                  gC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGa   >  1:782128/1‑36 (MQ=255)
                                  gC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGa   >  1:2468523/1‑36 (MQ=255)
                                  gC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGa   >  1:8278047/1‑36 (MQ=255)
                                   c.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGAc  <  1:3882465/36‑1 (MQ=255)
                                    |                                     
CTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGC.GGT.GG.TTGGGTTATCCCCGGTCAACATCACCAGAC  >  REL606/481475‑481545

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