Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 348,859 C→G 17.4% R349G (CGT→GGT)  mhpT → predicted 3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606348,8590C→G17.4% 4.8 ‑4.6 0.0471045mhpTpredicted 3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  new base (8/1):  ref base (15/23):  total (23/24)
Fisher's exact test strand distribution p-value = 1.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCGCAACAGGTGTGGGAACAGCCGTGGCGGTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTGGCCGGG  >  REL606/348824‑348894
                                   |                                   
gcgcACGAGGTGTGGGAACAGCCGTGACGGTAGGGc                                     <  1:328148/36‑1 (MQ=25)
  gcAACAGGTGTGGGAACAGCCGTGGCGGGAGGGCGt                                   >  1:615496/1‑36 (MQ=255)
     agaGGTTTGGGAACCGCCGTGGCGGTAGGGCGTCTg                                <  1:497994/34‑1 (MQ=37)
     aCAGGTGTGGGAACAGCCGTGGCGGGAGGGGGGCTg                                >  1:4133887/1‑36 (MQ=255)
      cAGGTGTGGGAACAGCCGTGGCGGTGGGGGGTCTgg                               >  1:3061065/1‑36 (MQ=37)
       aGGTGTGGGAACAGCCGTGGCGGTAGGGGGGCTggg                              >  1:3096630/1‑36 (MQ=37)
        ggggtgGGAACAGCCGTGGCGGTAGGGGGTTTtggg                             >  1:6228377/4‑36 (MQ=19)
         gtgtGGGAACAGCCGTGGCGGTGGGGGGTCTggggg                            >  1:1832128/1‑36 (MQ=37)
         gtggGGGAACAGCCGTGGCGGTAGGGCGTCTggggg                            <  1:651377/36‑1 (MQ=255)
         gtggGGGAACAGCCGTGGCGGTAGGGCGTCTggggg                            <  1:3376022/36‑1 (MQ=255)
         gggtgGGAACAGCCGTGGCGGTAGGGCGTCTggggg                            <  1:6020538/34‑1 (MQ=255)
          tgtgGGAACAGCCGTGGCGGGGGGGGGTCTGGTGGc                           >  1:2742952/1‑36 (MQ=255)
           gtgGGAACAGCCGTGGCGGTAGGGCGGCTggggggg                          >  1:5316449/1‑34 (MQ=37)
           gtgGGAAAAGCCGTGGCGGTAGGGGGGCTggggggc                          >  1:5475007/1‑34 (MQ=21)
            tgGGAACAGCCGTGGCGGTAGGGCGTCTGGGGGCTa                         <  1:3593788/36‑1 (MQ=255)
            ggggAACAGCCGTGGCGGTAGGGCGTCTGGGGGCTa                         <  1:6420786/35‑1 (MQ=255)
             gggAACAGCCGTGGCGGTAGGGCGTCTGGGGGCTAt                        <  1:5629331/36‑1 (MQ=255)
             gggAACAGCCGTGGCGGTAGGGCGTCTGGGGGCTAt                        <  1:3170094/36‑1 (MQ=255)
              ggAACAGCCGTGGCGGTAGGGCGTCTGGGGGCTATg                       <  1:1495498/36‑1 (MQ=255)
               gaccAGCCGTGGCGGTAGGGCGTCTGGGGGCTATGa                      <  1:5597440/33‑1 (MQ=255)
                aaCAGCCGTGGCGGTAGGGCGCCTGGGGGGTATGAg                     >  1:6771641/1‑36 (MQ=38)
                aaCAGCCGTGGCGGTAGGGCGCCTGGGGGCTATGAg                     >  1:7573681/1‑36 (MQ=255)
                aaCAGCCGTGGCGGGAGGGGGTGTGGGGGGGATGAg                     >  1:2546488/1‑36 (MQ=255)
                  cAGCCGTGGCGGTAGGGCGTCTGGGGGCTATGAGCg                   <  1:4017799/36‑1 (MQ=255)
                   aGCCGTGGCGGTTGGGCGTCTGGGGGGTTTGAGCgg                  >  1:6449319/1‑36 (MQ=25)
                   aGCCGTGGCGGTAGGGCGTCTGGGGGCTATGAGCgg                  <  1:3933180/36‑1 (MQ=255)
                     ccGTGGCGGTAGGGCGTCTGGGGGCTATGAGCGGTc                <  1:3898569/36‑1 (MQ=255)
                     ccGTGGCGGTAGGGCGTCTGGGGGCTATGAGCGGTc                <  1:6311648/36‑1 (MQ=255)
                     ccGTGGCGGTAGGGCGTCTGGGGGCTATGAGCGGTc                <  1:3262931/36‑1 (MQ=255)
                      cGTGGCGGTAGGGCGTCTGGGGGCTATGAGCGGTcc               <  1:5657761/36‑1 (MQ=255)
                        tGGCGGTAGGGCGTCTGGGGGCTATGAGCGGTCCGt             <  1:6712472/36‑1 (MQ=255)
                        tGGCGGTAATGGCTCTGGGGGCTTTGAGCGGTCCGt             <  1:1669254/36‑1 (MQ=16)
                         ggCGGTAGGGCGTCTGGGGGCTATGAGCGGTCCGtt            <  1:3825982/36‑1 (MQ=255)
                         ggCGGTAGGGCGTCTGGGGGCTATGAGCGGTCCGtt            >  1:254467/1‑36 (MQ=255)
                          gCGGTGGGGCGTCTGGGGGCTATGAGCGGTCCGTTa           <  1:2671572/36‑1 (MQ=255)
                           cGGTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTAc          >  1:5853526/1‑36 (MQ=255)
                            ggTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTACt         >  1:2400279/1‑36 (MQ=255)
                            ggTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTACt         <  1:5871348/36‑1 (MQ=255)
                             gTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTg        >  1:3033362/1‑36 (MQ=255)
                             gTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTg        >  1:1467148/1‑36 (MQ=255)
                             gTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTg        <  1:1104078/36‑1 (MQ=255)
                               aGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTGGc      <  1:4911231/36‑1 (MQ=255)
                               aGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTGGc      >  1:4049062/1‑36 (MQ=255)
                               aGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTGGc      >  1:5035068/1‑36 (MQ=255)
                                gggCGTCTGGGGGCTATGAGCGGTCCGTTACTGGcc     <  1:4658126/36‑1 (MQ=255)
                                gggCGTCTGGGGGCTATGAGCGGTCCGTTACTGGcc     >  1:1708449/1‑36 (MQ=255)
                                gggCGTCTGGGGGCTATGAGCGGTCCGTTACTGGcc     >  1:4312157/1‑36 (MQ=255)
                                 ggCGTCTGGGGGCTATGAGCGGTCCGTTACTGGCCg    >  1:7043342/1‑36 (MQ=255)
                                   cGTCTGGGGGCTATGAGCGGTCTGTCAATGGCCggg  >  1:4520351/1‑36 (MQ=25)
                                   cGTCTGGGGGCTATGAGCGGTCCGTTACTGGCCggg  >  1:1415989/1‑36 (MQ=255)
                                   |                                   
GCGCAACAGGTGTGGGAACAGCCGTGGCGGTAGGGCGTCTGGGGGCTATGAGCGGTCCGTTACTGGCCGGG  >  REL606/348824‑348894

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