Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,977,129 C→A 16.2% Q59K (CAA→AAA)  yedX → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,977,1290C→A16.2% 5.6 ‑3.6 0.048175yedXhypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/8):  ref base (21/18):  total (22/26)
Fisher's exact test strand distribution p-value = 2.76e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TCACTCTTGAAAAGAAGGCGGATAACGGCTGGTTACAACTTAATACCGCCAAAACAGATAAGGATGGACGA  >  REL606/1977094‑1977164
                                   |                                   
tCACTCTTGAAAAGAAGGCGGATAACGGCTGGTTAc                                     >  1:214541/1‑36 (MQ=255)
tCACTCTTGAAAAGAAGGCGGATAACGGCTGGTTAc                                     >  1:7702959/1‑36 (MQ=255)
 cACTCTTGAAAAGAAGGCGGATAACGGCTGGTTACa                                    >  1:4444674/1‑36 (MQ=255)
 cACTCTTGAAAAGAAGGCGGATAACGGCTGGTTACa                                    <  1:804951/36‑1 (MQ=255)
  aCTCTTGAAAAGAAGGCGGATGACGTCTGGTTACaa                                   <  1:8211998/36‑1 (MQ=37)
   ctctTGAAAAGAAGGCGGATAACGGCTGGTTACAAc                                  <  1:6918441/36‑1 (MQ=255)
   ctattGAAAAGAAGGCGGATAACGGCTGGTTACAAc                                  <  1:3751528/33‑1 (MQ=255)
    tctTGAAAAGAAGGCGGATAACGGCTGGTTACAACt                                 <  1:1321087/36‑1 (MQ=255)
    tctTGAAAAGAAGGCGGATAACGGCTGGTTACAACt                                 >  1:4858766/1‑36 (MQ=255)
    tctTGAAAAGAAGGCGGATAACGGCTGGTTACAACt                                 <  1:4754388/36‑1 (MQ=255)
    tctTGAAAAGAAGGCGGATAACGGCTGGTTACAACt                                 <  1:2971359/36‑1 (MQ=255)
    tattGAAAAGAAGGCGGATAACGGCTGGTTACAACt                                 <  1:1584559/34‑1 (MQ=255)
     ctTGAAAAGAAGGCGGATAACGGCTGGTTACAACtt                                <  1:499250/36‑1 (MQ=255)
      ttGAAAAGAAGGCGGATAACGGCTGGTTACCACTTa                               >  1:5368949/1‑36 (MQ=255)
      ttGAAAAGAAGGCGGATAACGGCTGGTTACAACTTa                               >  1:3739465/1‑36 (MQ=255)
       tGAAAAGAAGGCGGATAACGGCTGGTTACAACTTaa                              <  1:625016/36‑1 (MQ=255)
        gAAAAGAAGGCGGATAACGGCTGGTTACAACTTAAt                             <  1:3982132/36‑1 (MQ=255)
        gAAAAGAAGGCGGATAACGGCTGGTTACAACTTAAt                             <  1:3465274/36‑1 (MQ=255)
         aaaaGAAGGCGGATAACGGCTGGTTACAACAGAATa                            >  1:1080324/1‑36 (MQ=37)
          aaaGAAGGCGGATAACGGCTGGTTACAACTTAATAc                           <  1:4253911/36‑1 (MQ=255)
          aaaGAAGGCGGATAACGGCTGGTTACAACTTAATAc                           >  1:2526063/1‑36 (MQ=255)
           aagaagGCGGATAACGGCTGGTTACAACTTAATAcc                          >  1:3289682/1‑36 (MQ=255)
           aagaagGCGGATAAAGGCTGGGTAAAACTTAATAcc                          >  1:420876/1‑36 (MQ=25)
            agaagGCGGATAACGGCTGGTTACAACTTAATACCg                         <  1:5079431/36‑1 (MQ=255)
             gaagGCGGATAACGGCTGGTTACAACTTAATACCGc                        <  1:6843279/36‑1 (MQ=255)
              aagGCGGATAACGGCTGGTTACAACTTAATACCGcc                       <  1:7946263/36‑1 (MQ=255)
              aagGCGGATAACGGCTGGTTACAACTTAATACCGcc                       <  1:786602/36‑1 (MQ=255)
              aagGCGGATAACGGCTGGTTACAACTTAATACCGcc                       >  1:1843452/1‑36 (MQ=255)
              aagGCGGATAACGGCTGGTTACAACTTAATACCGcc                       <  1:5909516/36‑1 (MQ=255)
              aagGCGGATAACGGCTGGTTACAACTTAATACCGcc                       >  1:2006900/1‑36 (MQ=255)
                  cGGATAACGGCTGGTTAAAACTTAATACCGCCaaaa                   <  1:1638992/36‑1 (MQ=255)
                   ggATAACGGCTGGTTACAACTTAATACCGCCaaaaa                  >  1:6164512/1‑35 (MQ=255)
                   ggATAACGGCTGGTTACAACTTAATACCGCCAAAAc                  >  1:6243929/1‑36 (MQ=255)
                   ggATAACGGCTGGTTACAACTTAATACCGCCAAAAc                  >  1:6925869/1‑36 (MQ=255)
                     gaaaCGGGTGGTAAACACTTAAAACCGCCAAAACAg                <  1:4373589/34‑1 (MQ=12)
                     aaaaCGGCTGGTTAAAACTAAATACCGCCAAAACAg                <  1:6530883/34‑1 (MQ=37)
                     aTAACGGCTGGTTACAACTTAATACCGCCAAAATAg                >  1:1675806/1‑36 (MQ=255)
                     aTAACGGCTGGTTACAACTTAATACCGCCAAAACAg                >  1:7241349/1‑36 (MQ=255)
                     aTAACGGCTGGTTACAACTTAATACCGCCAAAACAg                >  1:6308042/1‑36 (MQ=255)
                     aTAACGGCTGGGTACAAATTTATACAGCCagaaagg                >  1:6136525/1‑33 (MQ=2)
                      aaaCGGCTGGTTAAAACTAAATACCGCCAAAACAGa               <  1:462126/35‑1 (MQ=37)
                       aaCGGCTGGTTACAACTTAATACCGCCAAAACAGAt              >  1:4062144/1‑36 (MQ=255)
                       aaCGGCTGGTTAAAACTTAATACCGCCAAAACAGAt              <  1:1845473/36‑1 (MQ=255)
                        aCGGGTGGTTAAAACTAAATACCGCCAAAACAGATa             <  1:2564649/36‑1 (MQ=25)
                           gCTGGTTACAACTTAATACCGCCAAAACAGATAAgg          >  1:455844/1‑36 (MQ=255)
                           gCTGGTTAAAATTAAATACCGCCAAAACAGATAAgg          <  1:4446377/36‑1 (MQ=21)
                             tGGTTACAACTTAATACCGCCAAAACAGATAAGGAt        >  1:7535596/1‑36 (MQ=255)
                             tGGTTACAACTTAATACCGCCAAAACAGATAAGGAt        >  1:3422663/1‑36 (MQ=255)
                              ggTTAAAAATAAATACCGCCAAAACAGATAAGGATg       <  1:317393/36‑1 (MQ=21)
                                 tACAACTTAATACCGCCAAAACAGATAAGGATGGAc    >  1:7121130/1‑36 (MQ=255)
                                   cAACTTAATACCGCCAAATCGGATAAGGATGGACGa  >  1:1560354/1‑36 (MQ=37)
                                   |                                   
TCACTCTTGAAAAGAAGGCGGATAACGGCTGGTTACAACTTAATACCGCCAAAACAGATAAGGATGGACGA  >  REL606/1977094‑1977164

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 18 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG