Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 348,392 +G 6.9% coding (578/1212 nt) mhpT → predicted 3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606348,3921.→G 6.9% 4.3 ‑3.0 ‑0.129coding (578/1212 nt)mhpTpredicted 3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (5/22):  total (7/22)
Fisher's exact test strand distribution p-value = 5.17e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.79e-01

GGTGCCGCTATTAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCTGGCGAAAAACAGTCTGCGCCACCAC  >  REL606/348357‑348427
                                    |                                   
ggTGCCGCTATTAATGCGCTGGCTGCCGGAGTcggc                                      <  1:5109890/36‑1 (MQ=255)
ggTGCCGCTATTAATGCGCTGGCTGCCGGAGTcggc                                      >  1:417076/1‑36 (MQ=255)
  tGCCGCTATTAATGCGCTGGCTGCCGGAGTcggc.gg                                   <  1:3003899/36‑1 (MQ=255)
   gCCGCTATTAATGCGCTGGCTGCCGGAGTcggc.ggg                                  >  1:6662604/1‑35 (MQ=255)
   gCCGCTATTAATGCGCTGGCTGCCGGAGTCGGC.GGt                                  <  1:2551425/36‑1 (MQ=255)
     cacTATTAATGCGCTGGCTGCCGGAGTCGGC.GTttt                                <  1:3014289/34‑1 (MQ=37)
     cGCTATTAATGCGCTGGCTGCCGGAGTCGGC.GGttt                                <  1:6704677/36‑1 (MQ=255)
       cTATTAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTc                              <  1:6254352/36‑1 (MQ=255)
        tATTAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTCg                             <  1:2445522/36‑1 (MQ=255)
         aTTAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGc                            <  1:4852417/36‑1 (MQ=255)
          ttAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCt                           <  1:4745224/36‑1 (MQ=255)
          ttAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCt                           <  1:3036774/36‑1 (MQ=255)
           tAATGCGCTGGCTGCCGGCGAAGGC.GGGGTTAGCTg                          >  1:799742/1‑36 (MQ=9)
           tAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCTg                          <  1:4898130/36‑1 (MQ=255)
            aaagcgcTGGCTGCCGGATTCGGC.GGTTTTCGCTgg                         <  1:7780900/33‑1 (MQ=25)
            aaTGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCTgg                         <  1:6533203/36‑1 (MQ=255)
            aaTGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCTgg                         <  1:5553458/36‑1 (MQ=255)
             aTGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCTGGc                        <  1:923949/36‑1 (MQ=255)
             aTGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCTGGc                        <  1:1001137/36‑1 (MQ=255)
               tcgcTGGCTGCCGGAGTCGGC.GGTTTTCGCTGGCGa                      <  1:6526515/35‑1 (MQ=255)
               gcgcTGGCTGCCGGAGTCGGC.GGTTTTCGCTGGCGa                      <  1:6109705/36‑1 (MQ=255)
                 gctggctgCCGGAGTCGGC.GGTTTTCGCTGGCGaaa                    <  1:1718392/36‑1 (MQ=255)
                  ctggctgCCGGAGTCGGC.GGTTTTTGCTGGGGgaaa                   >  1:6044077/1‑36 (MQ=25)
                  ctggctgCCGGAGTCGGC.GGGTTTTTCTTGCGgaaa                   >  1:3728469/1‑36 (MQ=12)
                   tggctgCCGGAGTCGGCGGGTTTTTGTTGGGGaaaa                   >  1:4388842/1‑36 (MQ=12)
                   tggctgCCGGAGTCGGC.GGTTTTCGCTGGCGaaaaa                  <  1:2590298/36‑1 (MQ=255)
                    ggctgCCGGAGTCGGCGGGTTTTGGTTGGGGaaaac                  >  1:6737800/1‑35 (MQ=12)
                     tctgCCGGAGTCGGC.GGTTTTCGCTGGCGAAAACCa                <  1:4975707/35‑1 (MQ=255)
                     gctgCCGGAGTCGGC.GGTTTTTTCTGGGGGAAAACa                >  1:6418592/1‑36 (MQ=16)
                     gctgCCGGAGTCGGC.GGATTTCGGTGGCGGAAAACa                >  1:7591422/1‑36 (MQ=21)
                       tgCCGGAGTCGGC.GGTTTTCGCTGGCGAAAAACAGt              <  1:4912282/36‑1 (MQ=255)
                        gCCGGAGTCGGC.GGTTTTCGCTGGCGAAAAACAGTc             <  1:2579337/36‑1 (MQ=255)
                           ggAGTCGGC.GGTTTTCGCTGGCGAAAAACAGTCTgc          <  1:5265941/36‑1 (MQ=255)
                                 ggc.ggTTTTCGCTGGCGAAAAACAGTCTGCGccacc    <  1:328293/36‑1 (MQ=255)
                                   c.ggTTTTCGCTGGCGAAAAACAGTCTGCGCcaccac  >  1:1415381/1‑36 (MQ=255)
                                    |                                   
GGTGCCGCTATTAATGCGCTGGCTGCCGGAGTCGGC.GGTTTTCGCTGGCGAAAAACAGTCTGCGCCACCAC  >  REL606/348357‑348427

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