Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,134,189 G→T 17.8% D414Y (GAT→TAT)  yegT → predicted nucleoside transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,134,1890G→T17.8% 5.3 ‑3.8 0.0531240yegTpredicted nucleoside transporter
Reads supporting (aligned to +/- strand):  new base (11/0):  ref base (26/16):  total (37/16)
Fisher's exact test strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAA  >  REL606/2134157‑2134223
                                |                                  
ccGACAACGAAATTACGGCTAACAAAGGCGATTATc                                 >  1:4346569/1‑36 (MQ=19)
   aCAACGAAATTACGGCTATCAAGGTCGATTATcgcg                              >  1:7411304/1‑36 (MQ=255)
   aCAACGAAATTACGGCTATCAAGGTCGATTAGcgcg                              >  1:3318631/1‑36 (MQ=38)
    cAACGAAATTACGGCTATCAAGGTCGATTATCGCGa                             >  1:4428941/1‑36 (MQ=255)
    cAACGAAATTACGGCTATCAAGGTCGATTATCGCGa                             >  1:2424255/1‑36 (MQ=255)
    cAACGAAATTACGGCTATCAAGGTCGATGATCGCGa                             <  1:3209485/36‑1 (MQ=255)
     aaCGAAATTACGGCTATCAAGGTCGATTATcgcgct                            >  1:6685718/1‑34 (MQ=37)
     aaCGAAATTACGGCTATCAAGGTCGATGATCGCGat                            <  1:367569/36‑1 (MQ=255)
     aaCGAAATTACGGCTATCAAGGTCGATGATCGCGat                            <  1:3492094/36‑1 (MQ=255)
     aaCGAAATTACGGCTATCAAGGGCGATGATCGCGat                            >  1:5831/1‑36 (MQ=255)
      aCGAAATTACGGCTATCAAGGTCGATTATCGCGata                           >  1:1818177/1‑36 (MQ=255)
      aCGAAATTACGGCTATCAAAGTCGATGATCGCGata                           >  1:5704893/1‑36 (MQ=255)
      aCGAAATTACGGCTATAAAGGGCGAGTATCGCGata                           >  1:1897726/1‑36 (MQ=19)
       cGAAATTACGGCTATCAAGGTCGATGATCGCGatat                          >  1:54363/1‑36 (MQ=255)
        gAAATTACGGCTATCAAGGTCGATGATCGCGATAtt                         <  1:4037140/36‑1 (MQ=255)
         aaaTTACGGCTATCAAGGTCGATGATCGCGATATTg                        >  1:7096913/1‑36 (MQ=255)
         aaaTTACGGCTATCAAGGTCGATGATCGCGATATTg                        >  1:5945038/1‑36 (MQ=255)
          aaTTACGGCTATCAAGGTCGATGATCGCGATATTGc                       >  1:7268310/1‑36 (MQ=255)
           aTTACGGCTATCAAGGTCGATTATCGCGATATTGCg                      >  1:7350432/1‑36 (MQ=255)
           aTTACGGCTATCAAGGTCGATGATCGCGATATTGCg                      >  1:449250/1‑36 (MQ=255)
            ttACGGCTATCAAGGTCGATGATCGCGATATTGCGt                     <  1:1504084/36‑1 (MQ=255)
             tACGGCTATCAAGGTCTATGATCGCGATATTGCGtt                    <  1:6441566/36‑1 (MQ=255)
             tACGGCTATCAAGGTCGATGATCGCGATATTGCggt                    >  1:5838571/1‑34 (MQ=255)
             tACGGCTATCAAGGTCGATGATCGCGATATTGCGtt                    >  1:5080691/1‑36 (MQ=255)
              aCGGCTATCAAGGTCGATTATCGCGATATTGCGTTg                   >  1:7266128/1‑36 (MQ=255)
              aCGGCTATCAAGGTCGATTATCGCGATATTGCGTTg                   >  1:291542/1‑36 (MQ=255)
              aCGGCTATCAAGGTCGATGATCGCGATATTGCGTTg                   >  1:6529961/1‑36 (MQ=255)
               cGGCTATCAAGGTCGATGATCGCGATATTTTGTTGa                  >  1:5926278/1‑36 (MQ=37)
               cGGCTATCAAGGTCGATGATCGCGATATTGCGTTGa                  <  1:1178896/36‑1 (MQ=255)
               cGGCTATCAAGGTCGATGATCGCGATATTGCGTTGa                  <  1:4497033/36‑1 (MQ=255)
               cGGCTATCAAGGTCGATGATCGCGATATTGCGTTGa                  <  1:7065053/36‑1 (MQ=255)
                ggCTATCAAGGTCGATGATCGCGATATTGCGTTTac                 >  1:5396669/1‑36 (MQ=255)
                ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                 >  1:4904792/1‑36 (MQ=255)
                ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                 >  1:1859907/1‑36 (MQ=255)
                ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                 <  1:5221902/36‑1 (MQ=255)
                ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                 <  1:7433210/36‑1 (MQ=255)
                 gCTATCAAGGTCGATGATCGCGATATTGCGTTGaca                <  1:1025051/36‑1 (MQ=255)
                  cTATCAAGGTCGATGATCGCGATATTGCGTTGacac               <  1:7033518/36‑1 (MQ=255)
                   tATCAAGGTCGATGATCGCGATATTGCGTTGAcaca              <  1:773962/36‑1 (MQ=255)
                   tATCAAGGTCGATGATCGCGATATTGCGTTGAcaca              >  1:5572730/1‑36 (MQ=255)
                    aTCAAGGTCGATGATCGCGATATTGCGTTGACACaa             >  1:4142078/1‑36 (MQ=255)
                    aTCAAGGTCGATGATCGCGATATTCCGTTGACACaa             >  1:3166910/1‑36 (MQ=255)
                       aaGGTCGATGATCGCGATATTGCGTTGACACAAggg          >  1:2857529/1‑36 (MQ=255)
                        aGGTCGATGATCGCGATATTGCGTTGACACAAgggg         >  1:6533716/1‑36 (MQ=255)
                        aGGTCGATGATCGCGATATTGAGTTGAAACAAgggg         >  1:1350791/1‑36 (MQ=37)
                         ggTCGATGATCGCGATATTGCGTTGACACAAGGGGa        <  1:422076/36‑1 (MQ=255)
                         ggTCGATGATCGCGATATTGCGTTGACACAAGGGGa        >  1:2414725/1‑36 (MQ=255)
                         ggTCGATGATCGCGATATTGCGTTGACACAAGGGGa        <  1:203262/36‑1 (MQ=255)
                          gTCGATGATCGCGATATTGCGTTGACACAAGGGGaa       >  1:1067581/1‑36 (MQ=255)
                          gTCGATGATCGCGATATTGCGTTGACACAAGGGGaa       >  1:6878317/1‑36 (MQ=255)
                          gTCGATGATCGCGATATTGCGTTGACACAAGGGGaa       >  1:4393179/1‑36 (MQ=255)
                           tCGATGATCGCGATATTGCGTTGACACAAGGGGAAg      >  1:5105505/1‑36 (MQ=255)
                            cGATGATCGCGATATTGCGTTGACACAAGGGGAAGt     >  1:4188598/1‑36 (MQ=255)
                             gatgatCGCGATATTGCGTTGACACAAGGGGAAGtt    >  1:4715803/1‑36 (MQ=255)
                             gatgatCGCGATATTGCGTTGACACAAGGGGAAGtt    <  1:5620136/36‑1 (MQ=255)
                              atgatCGCGATATTGCGTTGACACAAGGGGAAGTTa   >  1:6789819/1‑36 (MQ=255)
                              atgatCGCGATATTGCGTTGACACAAGGGGAAGTTa   >  1:748216/1‑36 (MQ=255)
                               tgatCGCGATATTGCGTTGACACAAGGGGAAGTTaa  <  1:1315035/36‑1 (MQ=255)
                               tgatCGCGATATTGCGTTGACACAAGGGGAAGTTaa  <  1:849349/36‑1 (MQ=255)
                                |                                  
CCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAA  >  REL606/2134157‑2134223

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG