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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,134,1890G→T14.0% 3.5 ‑3.7 0.0601240yegTpredicted nucleoside transporter
Reads supporting (aligned to +/- strand):  new base (11/0):  ref base (31/18):  total (42/18)
Fisher's exact test strand distribution p-value = 2.45e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AATCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAA  >  REL606/2134154‑2134224
                                   |                                   
aaTCCGACAACGAAATTACGGCTATCAAGGTCgatg                                     <  1:3687779/36‑1 (MQ=255)
 aTCCGACAACGAAATTACGGCTATCAAGGTCgatga                                    <  1:7630526/36‑1 (MQ=255)
 aTCCGACAACGAAATTACGGCTATCAAGGTCgatga                                    >  1:3888543/1‑36 (MQ=255)
 aTCCGACAACGAAATTACGGCTATCAAGGTCgatga                                    <  1:1460896/36‑1 (MQ=255)
 aTCCGACAACGAAATTACGGCTATCAAGGTCgatga                                    >  1:7750032/1‑36 (MQ=255)
  tCCGACAACGAAATTACGGCTATCAAGGTCgattat                                   >  1:7718618/1‑36 (MQ=255)
  tCCGACAACGAAATTACGGCTATCAAGGTCgatgat                                   <  1:5785329/36‑1 (MQ=255)
   ccGACAACGAAATTACGGCTATCAAGGTCGATGATc                                  <  1:2287067/36‑1 (MQ=255)
    cGACAACGAAATTACGGCTATCAAGGTCGATGATcg                                 <  1:3175022/36‑1 (MQ=255)
    cGACAACGAAATTACGGCTATCAAAGTCGATTATcg                                 >  1:1640089/1‑36 (MQ=38)
     gACAACGAAATTACGGCTATCAAGGTCGATTATcgc                                >  1:5136682/1‑36 (MQ=255)
      aCAACGAAATTACGGCTATCAAGGTCGATGATcgcg                               >  1:2172804/1‑36 (MQ=255)
      aCAACGAAATTACGGCTATCAAGGTCGATGATcgcg                               >  1:5693259/1‑36 (MQ=255)
       cAACGAAATTACGGCTATCAAGGTCGATGATCGCGa                              <  1:8049419/36‑1 (MQ=255)
        aaagAAATTACGGCTATCAAGGTCGATGATCGCGat                             <  1:671538/33‑1 (MQ=255)
        aaCGAAATTACGGCTAACAAGGTCGATTATAGCGat                             >  1:2087551/1‑36 (MQ=25)
         aCGAAATTACGGCTATAAAGGTCGATGATTGCGata                            >  1:1320227/1‑36 (MQ=37)
         aCGAAATTACGGCTAGCAAGGTCGATTATCGCGata                            >  1:2013941/1‑36 (MQ=38)
         aCGAAATTAAGGCTATCAAGGTCGATGATCGCGata                            <  1:7129047/36‑1 (MQ=255)
          cGAAATTACGGCTATCAAGGTCGATTATCGCGatat                           >  1:1390132/1‑36 (MQ=255)
           gAAATTACGGCTATCAAGGTCGATGATCGCGATAtt                          <  1:3288031/36‑1 (MQ=255)
           gAAATTACGGCTATCAAGGTCGAGGAAAGCGatagt                          >  1:1941230/1‑34 (MQ=21)
           gAAATTACGGCTATCAAAGACGATGTNNNCGCTAtt                          >  1:1691984/1‑36 (MQ=255)
            aaaTTACGGCTATCAAGGTCGATGATCGCGATATTg                         <  1:7119049/36‑1 (MQ=255)
             aaTTACGGCTATCATTTTTGATGATCGGGATATTGc                        <  1:6763625/36‑1 (MQ=12)
             aaTTACGGCTATCAAGGTCGATTATTGCGGTAtttc                        >  1:4697152/1‑34 (MQ=21)
             aaTTACGGCTATCAAGGTCGATTATCGCGATAtttc                        >  1:4925467/1‑34 (MQ=38)
              aTTACGGCTATCACGGACTATGTTCGCGATTTTgtg                       >  1:2274281/1‑34 (MQ=12)
              aTTACGGCTATCAAGGTCGATTATCGCGATATTGCg                       >  1:836526/1‑36 (MQ=255)
               ttACGGCTATCAAGGTCGATTATCGCGATATTGCGt                      >  1:3843934/1‑36 (MQ=255)
               ttACGGCTATCAAGGTCGATTATCGCGATATTGCGt                      >  1:1629216/1‑36 (MQ=255)
               ttACGGCTATCAAGGTCGATGATCGCGATATTGCGt                      >  1:2827754/1‑36 (MQ=255)
               ttACGGCTATCAAGGTCGATGATCGCGATATTGCGt                      >  1:8208782/1‑36 (MQ=255)
               atACGGCTATCAAGGTCGATGATCGCGATATTGCGt                      <  1:2171519/35‑1 (MQ=255)
                tACGGCTATCAAGGTCGATGATCGCGATATTGCGtt                     <  1:1488133/36‑1 (MQ=255)
                 aCGGCTATCAAGGTCGATTATCGCGATATTGCGTTg                    >  1:8244419/1‑36 (MQ=255)
                 aCGGCTATCAAGGTCGATGATCGCGATATTGCGTTg                    <  1:3696771/36‑1 (MQ=255)
                 aCGGCTATCAAGGTCGATGATCGCGATATTGCGTTg                    >  1:591503/1‑36 (MQ=255)
                  cGGCTATCAAGGTCGATGATCGAGATATTGCGTTGa                   >  1:719234/1‑36 (MQ=255)
                   ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                  >  1:7190631/1‑36 (MQ=255)
                   ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                  >  1:6155533/1‑36 (MQ=255)
                   ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                  <  1:3562689/36‑1 (MQ=255)
                   ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                  >  1:1194623/1‑36 (MQ=255)
                   ggCTATCAAGGTCGATGATCGCGATATTGCGTTGac                  >  1:1201628/1‑36 (MQ=255)
                    gCTATCAAGGTCGATGATCGCGATATTGCGTTGaca                 >  1:8027078/1‑36 (MQ=255)
                     cTATCAAGGTCGATGATCGCGATATTGCGTTGacac                <  1:6847383/36‑1 (MQ=255)
                      tATCAAGGTCGATGATCGCGATATTGCGTTGAcaca               >  1:1558738/1‑36 (MQ=255)
                      tATCAAGGTCGATGATCGCGATATTGCGTTGAcaca               >  1:3948332/1‑36 (MQ=255)
                      tATCAAGGTCGATGATCGCGATATTGCGTTGAcaca               <  1:131666/36‑1 (MQ=255)
                       aTCAAGGTCGATGATCGCGATATTGCGTTGACACaa              >  1:4835350/1‑36 (MQ=255)
                       aTCAAGGTCGATGATCGAGATATTGCGTTGACACaa              >  1:4841353/1‑36 (MQ=255)
                        tCAAGGTCGATGATCGCGATATTGCGTTGACACAAg             >  1:7128117/1‑36 (MQ=255)
                        tCAAGGTCGATGATCGCGATATTGCGTTGACACAAg             <  1:8389100/36‑1 (MQ=255)
                        tCAAGGTCGATGATCGCGATATTGCGTTGACACAAg             >  1:952285/1‑36 (MQ=255)
                         cAAGGTCGATGATCGCGATATTGCGTTGACACAAgg            >  1:6519712/1‑36 (MQ=255)
                         cAAGGTCGATGATCGCGATATTGCGTTGACACAAgg            >  1:7549794/1‑36 (MQ=255)
                         cAAGGTCGATGANNGNGATANTGCGTTGACACAAgg            >  1:2312983/1‑36 (MQ=39)
                          aaGGTCGATGATCGCGATATTGCGTTGACACAAggg           >  1:3798937/1‑36 (MQ=255)
                          aaGGTCGATGATCGCGATATTGCGTTGACACAAggg           >  1:7645486/1‑36 (MQ=255)
                           aGGTCGATGATCGCGATATTGCGTTGACACAAgggg          <  1:5297737/36‑1 (MQ=255)
                           aGGTCGATGATCGCGATATTGCGTTGACACAAgggg          <  1:549033/36‑1 (MQ=255)
                            ggTCGATGATCGCGATATTGCGTTGACACAAGGGGa         <  1:2848930/36‑1 (MQ=255)
                            ggTCGATGATCGCGATATTGCGTTGACACAAGGGGa         >  1:6917759/1‑36 (MQ=255)
                            ggTCGATGATCGCGATATTGCGTTGACACAAGGGGa         >  1:1365492/1‑36 (MQ=255)
                            ggTCGATGATCGCGATATTGCGTTGACACAAGGGGa         >  1:5750009/1‑36 (MQ=255)
                              tCGATGATCGCGATATTGCGTTGACACAAGGGGAAg       >  1:7793783/1‑36 (MQ=255)
                               cGATGATCGCGATATTGCGTTGACACAAGGGGAAGt      >  1:7300032/1‑36 (MQ=255)
                                gatgatCGCGATATTGCGTTGACACAAGGGGAAGtt     <  1:754376/36‑1 (MQ=255)
                                gatgatCGCGATATTGCGTTGACACAAGGGGAAGtt     >  1:3707275/1‑36 (MQ=255)
                                 atgatCGCGATATTGCGTTGACACAAGGGGAAGTTa    >  1:7903832/1‑36 (MQ=255)
                                 atgatCGCGATATTGCGTTGACACAAGGGGAAGTTa    >  1:6904558/1‑36 (MQ=255)
                                  tgatCGCGATATTGCGTTGACACAAGGGGAAGTTaa   >  1:3300373/1‑36 (MQ=255)
                                   gatCGCGATATTGCGTTGACACAAGGGGAAGTTaaa  <  1:292464/36‑1 (MQ=255)
                                   gatCGCGATATTGCGTTGACACAAGGGGAAGTTaaa  >  1:5202761/1‑36 (MQ=255)
                                   |                                   
AATCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAA  >  REL606/2134154‑2134224

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