Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,542,103 +G 5.1% coding (319/1218 nt) yhhS ← predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,542,1031.→G 5.1% 3.4 ‑3.0 ‑0.039coding (319/1218 nt)yhhSpredicted transporter
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (7/30):  total (9/30)
Fisher's exact test strand distribution p-value = 4.86e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.05e-01

GCAAAGTAATAACAGGCTGATGACAGGCAGACTGGC.GGTTAATCCTGCCGTCAGATACCCCAGACCGCTCA  >  REL606/3542068‑3542138
                                    |                                   
gCAAAGTAATAACAGGCTGATGACAGGCAGACTGGc                                      <  1:195228/36‑1 (MQ=255)
gCAAAGTAATAACAGGCTGATGACAGGCAGACTGGc                                      >  1:3859040/1‑36 (MQ=255)
gCAAAGTAATAACAGGCTGATGACAGGCAGACTGGc                                      <  1:2876991/36‑1 (MQ=255)
gCAAAGTAATAACAGGCTGATGACAGGCAGACTGGc                                      <  1:3420116/36‑1 (MQ=255)
gCAAAGTAATAACAGGCTGATGACAGGCAGACTGGc                                      >  1:2952396/1‑36 (MQ=255)
gCAAAGTAATAACAGGCTGATGACAGGCAGACTGGc                                      <  1:3192149/36‑1 (MQ=255)
gCAAAGTAATAACAGGATGATGACAGGCAGACTGGc                                      <  1:1312668/36‑1 (MQ=255)
 cAAAGTAATAACAGGCTGATGTCAGGCAGACTGGC.g                                    >  1:5602756/1‑36 (MQ=255)
 cAAAGTAATAACAGGCTGATGACAGGCAGACTGGC.g                                    <  1:2039994/36‑1 (MQ=255)
  aaaGTAATAACAGGCTGATGACAGGCAGACTGGC.gg                                   >  1:5689/1‑36 (MQ=255)
   aaGTAATAACAGGCTGATGACAGGCAGACTGTC.gga                                  >  1:7413409/1‑35 (MQ=255)
   aaGTAATAACAGGCTGATGACAGGCAGACTGGC.ggg                                  >  1:897248/1‑35 (MQ=255)
   aaGTAATAACAGGCTGATGACAGGCAGACTGGC.ggg                                  >  1:4951411/1‑35 (MQ=255)
   aaGTAATAACAGGCTGATGACAGGCAGACTGGC.ggg                                  >  1:391090/1‑35 (MQ=255)
   aaGTAATAACAGGCTGATGACAGGCAGACTGGC.ggg                                  >  1:1511374/1‑35 (MQ=255)
   aaGTAATAACAGGCTGATGACAGGCAGACTGGC.GGt                                  <  1:2536120/36‑1 (MQ=255)
    aGTAATAACAGGCTGATGACAGGCAGACTGGC.gggt                                 >  1:1757129/1‑34 (MQ=255)
    aGTAATAACAGGCTGATGACAGGCAGACTGGC.gggg                                 >  1:595909/1‑34 (MQ=255)
    aGTAATAACAGGCTGATGACAGGAAGACGGGC.gggg                                 >  1:4290350/1‑34 (MQ=37)
     gTAATAACAGGCTGATGACAGGCAGACTTGC.gggtg                                >  1:34840/1‑33 (MQ=37)
     gTAATAACAGGCTGATGACAGGCAGACTGGC.GGttt                                >  1:7284312/1‑35 (MQ=255)
     gTAATAACAGGCTGATGACAGGCAGACTGGC.GGTTa                                <  1:6991262/36‑1 (MQ=255)
      taataaCAGGCTGATGACAGGCAGACTGGC.GGtttt                               >  1:7610297/1‑34 (MQ=255)
      taataaCAGGCTGAGGACAGGCAGACTGGC.GGGTag                               >  1:626344/1‑35 (MQ=37)
       aataaCAGGCTGATGACAGGCAGACTGGC.GGTTAAt                              <  1:5863784/36‑1 (MQ=255)
       aataaCAGGCTGATGACAGGCAGACTGGC.GGTTAAt                              <  1:5135715/36‑1 (MQ=255)
       aataaCAGGCTGATGACAGGCAGACTGGC.GGTTAAt                              <  1:4624734/36‑1 (MQ=255)
        ataaCAGGCTGATGACAGGCAGATTGGCGGGTTAAt                              >  1:2148957/1‑36 (MQ=25)
        ataaCAGGCTGATGACAGGCAGACTGGC.GGTTAATc                             <  1:2959974/36‑1 (MQ=255)
        ataaCAGGCTGATGACAGGCAGACTGGC.GGTTAATc                             <  1:7127986/36‑1 (MQ=255)
         taaCAGGCTGATGACAGGCAGACTGGC.GGTTAATcc                            <  1:37089/36‑1 (MQ=255)
         taaCAGGCTGATGACAGGCAGACTGGC.GGTTAATcc                            <  1:6400774/36‑1 (MQ=255)
         taaCAGGCTGATGACAGGCAGACTGGC.GGTTAATcc                            <  1:3422041/36‑1 (MQ=255)
           aCAGGCTGATGACAGGCAGACTGGC.GGTTAATCCTg                          <  1:1886783/36‑1 (MQ=255)
            cAGGCTGATGACAGGCAGACTGGC.GGTTAATCCTGc                         <  1:42240/36‑1 (MQ=255)
            cAGGCTGATGACAGGCAGACTGGC.GGTTAATCCTGc                         <  1:3562059/36‑1 (MQ=255)
             atgCTGATGACAGGCAGACTGGC.GGTTAATCCTGcc                        <  1:7760479/34‑1 (MQ=255)
              ggCTGATGACAGGCAGACTGGC.GGTTAATCCTGCCg                       <  1:4191116/36‑1 (MQ=255)
              ggCTGATGACAGGCAGACTGGC.GGTTAATCCTGCCg                       <  1:3662710/36‑1 (MQ=255)
               gCTGATGACAGGCAGACTGGC.GGTTAATCCTGCCGt                      <  1:6941268/36‑1 (MQ=255)
                cTGATGACAGGCATACTGGC.GGTTAATCCTGCCGTc                     <  1:6239478/36‑1 (MQ=255)
                   atgaCAGGCAGACTGGCGGGTTATTCCT.CCGTctgc                  >  1:5044465/1‑33 (MQ=12)
                   atgaCAGGCAGACTGGC.GGTTAATCCTGCCGTCAGa                  <  1:552217/36‑1 (MQ=255)
                   atgaCAGGCAGACTGGC.GGGTTATCATGCCGTccgg                  >  1:1685935/1‑33 (MQ=16)
                     gaCAGGCAGACTGGC.GGGTTATTCTGCCGTCGGATa                >  1:5887302/1‑36 (MQ=19)
                        aGGCAGACTGGC.GGTTAATCCTGCCGTCAGATAccc             <  1:3564677/36‑1 (MQ=255)
                         ggCAGACTGGC.GGTTAATCCTGCCGTCGGATAcccc            <  1:17503/36‑1 (MQ=255)
                         ggCAGACTGGC.GGTTAATCCTGCCGTCAGATAcccc            <  1:7746180/36‑1 (MQ=255)
                          gCAGACTGGC.GTTTAATCCTGCCGTCAGATACCCCa           <  1:2886015/36‑1 (MQ=255)
                          gCAGACTGGC.GGTTAATCCTGCCGTCAGATACCCCa           <  1:4838560/36‑1 (MQ=255)
                          gCAGACTGGC.GGGTAATTCTGATGTGTGATAcctcc           >  1:4656019/1‑33 (MQ=2)
                            aGACTGGC.GGTTAATCCTGCCGTCAGATACCCCAGa         <  1:4162963/36‑1 (MQ=255)
                             gACTGTC.GGTTAAACCTGCCGTCAGATACCCCAGAc        <  1:5741373/36‑1 (MQ=37)
                              aCTGGC.GGTTAATCCTGCCGTCAGATACCCCAGAcc       <  1:3812203/36‑1 (MQ=255)
                               cTGGC.GGTTAATCCTGCCGTCAGATACCCCAGACCg      <  1:3223556/36‑1 (MQ=255)
                                 ggC.GGTTAATCCTGCCGTCAGATACCCCAGACCGCt    <  1:3241317/36‑1 (MQ=255)
                                  gC.GGTTAATCCTGCCGTCAGATACCCCAGACCGCTc   >  1:5232106/1‑36 (MQ=255)
                                  gC.GGTTAATCCTGCCGTCAGATACCCCAGACCGCTc   <  1:7443403/36‑1 (MQ=255)
                                   c.GGTTAATCCTGCCGTCAGATACCCCAGACCGCTCa  <  1:4850866/36‑1 (MQ=255)
                                   c.GGTTAATCCTGCCGTCAGATACCCCAGACCGCTCa  <  1:7230364/36‑1 (MQ=255)
                                    |                                   
GCAAAGTAATAACAGGCTGATGACAGGCAGACTGGC.GGTTAATCCTGCCGTCAGATACCCCAGACCGCTCA  >  REL606/3542068‑3542138

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG