Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,484,118 G→T 12.3% D60Y (GAC→TAC)  ydcU → predicted spermidine/putrescine transporter subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,484,1180G→T12.3% 3.2 ‑3.7 ‑0.046178ydcUpredicted spermidine/putrescine transporter subunit
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (19/21):  total (25/21)
Fisher's exact test strand distribution p-value = 2.47e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.43e-01

GCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGAC  >  REL606/1484084‑1484152
                                  |                                  
gctgACACTGTTATGGCAGGGATTTTATACTTTTGa                                   <  1:110013/36‑1 (MQ=255)
 ctgACACTGTTATGGCAGGGGTTTTTTACTTTTTAc                                  >  1:199899/1‑36 (MQ=37)
 ctgACACTGTTATGGCAGGGATTTTATACTTTTGAc                                  <  1:6627710/36‑1 (MQ=255)
   gACACTGTTATGGCAGGGATTTTATACTTTTGACGa                                <  1:1086991/36‑1 (MQ=255)
    acacTGTTATGGCAGGGGTTTTATACTTTTTACGAt                               >  1:2348318/1‑36 (MQ=37)
    acacTGTTATGGCAGGGATTTTATACTTTTGACGAt                               <  1:3491427/36‑1 (MQ=255)
     cacTGTTATGGCAGGGGTTTTATACTTTTTACGAtt                              >  1:3351047/1‑36 (MQ=37)
     cacTGTTATGGCAGGGATTTTTTACTTTTTTCGAtt                              >  1:1956057/1‑36 (MQ=25)
      acTGTTATGGCAGGGGTTTTATACTTTTTTCGAttt                             >  1:1360836/1‑36 (MQ=25)
      acTGTTATGGCAGGGATTTTATACTTTTGACGAttt                             <  1:2978917/36‑1 (MQ=255)
      acTGTTATGGCAGGGATTTTATAATTTTGACGAttt                             <  1:1617913/36‑1 (MQ=255)
       cTGTTATGGCAGGGATTTTATACTTTTGACGATTTc                            <  1:6393629/36‑1 (MQ=255)
          ttATGGCAGGGATTTTATACTTTTTACGATTTCAcc                         >  1:1947021/1‑36 (MQ=255)
          ttATGGCAGGGATTTTATACTTTTGACGATTTCAcc                         <  1:2486479/36‑1 (MQ=255)
           tATGGCAGGGATTTTATACTTTTGACGATTTCACCa                        >  1:6250256/1‑36 (MQ=255)
            aTGGCAGGGATTTTATACTATTGACGATTTCACCAt                       >  1:2647986/1‑36 (MQ=255)
             tGGCAGGGATTTTATACTTTTGACGACTTCACCATg                      >  1:4579534/1‑36 (MQ=255)
              ggCAGGGATTTTATACTTTTGACGATTTCACCATGt                     <  1:5641920/36‑1 (MQ=255)
              ggCAGGGATTTTATACTTTTGACGATTTCACCATGt                     <  1:551381/36‑1 (MQ=255)
              ggCAGGGATTTTATACTTTTGACGATTTCACCATGt                     <  1:86317/36‑1 (MQ=255)
               gCAGGGATTTTATACTTTTGACGATTTCACCATGTc                    <  1:1356869/36‑1 (MQ=255)
               gCAGGGATTTTATACTTTTGACGATTTCACCATGTc                    <  1:5036765/36‑1 (MQ=255)
                 aGGGATTTTATACTTTTGACGATTTCACCATGTCgg                  >  1:3655271/1‑36 (MQ=255)
                 aGGGATTTTATACTTTTGACGATTTCACCATGTCgg                  >  1:5408856/1‑36 (MQ=255)
                  gggATTTTATACTTTTGACGATTTCACCATGTCGGt                 >  1:4621789/1‑36 (MQ=255)
                   ggATTTTATACTTTTGACGATTTCACCATGTCGGTa                <  1:1582724/36‑1 (MQ=255)
                     aTTTTATACTTTTGACGATTTCACCATGTCGGTAAc              >  1:2431977/1‑36 (MQ=255)
                      ttttATACTTTTGACGATTTCACCATGTCGGTAAcc             >  1:2787799/1‑35 (MQ=255)
                      ttttATACTTTTGACGATTTCACCATGTCGGTAACg             >  1:4256532/1‑36 (MQ=255)
                        ttATACTTTTGACGATTTCACCATGTCGGTAACGcc           <  1:1791180/36‑1 (MQ=255)
                        ttATACTTTTGACGATTTCACCATGTCGGTAACGcc           >  1:4083610/1‑36 (MQ=255)
                        ttATACTTTTGACGATTTCACCATGTCGGTAACGcc           <  1:234167/36‑1 (MQ=255)
                         tataCTTTTGACGATTTCACCATGTCGGTAACGCCg          <  1:2902509/36‑1 (MQ=255)
                          ataCTTTTGACGATTTCACCATGCCGGTTACGCCgg         >  1:5734115/1‑36 (MQ=37)
                           taCTTTTGACGATTTCACCATGTCGGTAACGCCGGa        >  1:2142675/1‑36 (MQ=255)
                            aCTTTTGACGATTTCACCATGTCGGTAACGCCGGaa       >  1:6046468/1‑36 (MQ=255)
                            aCTTTTGACGATTTCACCATGTCGGTAACGCCGGaa       >  1:5343671/1‑36 (MQ=255)
                            aCTTTTGACGATTTCACCATGTCGGTAACGCCGGaa       >  1:6449970/1‑36 (MQ=255)
                              ttttGACGATTTCACCATGTCGGTAACGCCGGAACt     >  1:5374213/1‑36 (MQ=255)
                              ttttGACGATTTCACCATGTCGGTAACGCCGGAACt     >  1:5779643/1‑36 (MQ=255)
                              ttttGACGATTTCACCATGTCGGTAACGCCGGAACt     <  1:6538408/36‑1 (MQ=255)
                               tttGACGATTTCACCATGTCGGTAACGCCGGAACTg    >  1:7417340/1‑36 (MQ=255)
                               tttGACGATTTCACCATGTCGGTAACGCCGGAACTg    <  1:3376851/36‑1 (MQ=255)
                                ttGACGATTTCACCATGTCGGTAACGCCGGAACTGa   >  1:4417251/1‑36 (MQ=255)
                                 tGACGATTTCACCATGTCGGTAACGCCGGAACTGAc  <  1:3812971/36‑1 (MQ=255)
                                 tGACGATTTCACCATGTCGGTAACGCCGGAACTGAc  <  1:5671248/36‑1 (MQ=255)
                                  |                                  
GCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGAC  >  REL606/1484084‑1484152

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