Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,409,316 T→G S41S (TCT→TCG yheS → fused predicted transporter subunits of ABC superfamily: ATP‑binding components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,409,3160TG98.7% 292.2 / ‑4.1 78S41S (TCT→TCGyheSfused predicted transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (28/49);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CCCCGGGCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCTACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGACGGCGGCAG  >  REL606/3409268‑3409363
                                                |                                               
ccccGGGCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGAc                                               <  1:595826/51‑1 (MQ=255)
  ccGGGCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGAccc                                             >  1:6228298/1‑51 (MQ=255)
   cGGGCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCct                                            >  1:4920024/1‑51 (MQ=255)
   cGGGCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCct                                            <  1:7121797/51‑1 (MQ=255)
   cGGGCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCct                                            >  1:2195178/1‑51 (MQ=255)
   cGGGCAGAAAGTCGGCCTGGGGGGTAAAAACGGCTGTGGTAAATCGAccc                                             >  1:4597331/1‑50 (MQ=39)
    gggCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctg                                           <  1:3926404/51‑1 (MQ=255)
    gggCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctg                                           >  1:2038028/1‑51 (MQ=255)
    gggCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctg                                           >  1:2341384/1‑51 (MQ=255)
     ggCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgc                                          <  1:5615653/51‑1 (MQ=255)
      gCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgct                                         <  1:6908931/51‑1 (MQ=255)
      gCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgct                                         <  1:5493483/51‑1 (MQ=255)
      gCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgct                                         >  1:6701251/1‑51 (MQ=255)
      gCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgct                                         <  1:4555637/51‑1 (MQ=255)
      gCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgct                                         <  1:3152075/51‑1 (MQ=255)
      gCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgct                                         <  1:6517557/51‑1 (MQ=255)
       cAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCctgctg                                        <  1:4671442/51‑1 (MQ=255)
        ggAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTgg                                       <  1:1501041/50‑1 (MQ=255)
           aaGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCAt                                    <  1:4508624/51‑1 (MQ=255)
           aaGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCAt                                    <  1:255296/51‑1 (MQ=255)
            aGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCAtt                                   <  1:4486582/51‑1 (MQ=255)
            aGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCAtt                                   >  1:4916078/1‑51 (MQ=255)
            aGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCAtt                                   <  1:2942861/51‑1 (MQ=255)
             gTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTg                                  <  1:2778469/51‑1 (MQ=255)
             gTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTg                                  <  1:4845318/51‑1 (MQ=255)
              tCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGc                                 <  1:5089127/51‑1 (MQ=255)
                ggCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTg                               >  1:6439318/1‑51 (MQ=255)
                 gCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGa                              <  1:6354502/51‑1 (MQ=255)
                 gCCTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGa                              >  1:3775654/1‑51 (MQ=255)
                  ccTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGaa                             <  1:3705676/51‑1 (MQ=255)
                  ccTGGTGGGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGaa                             <  1:4784976/51‑1 (MQ=255)
                    tggtggGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGaaaa                           >  1:4614945/1‑51 (MQ=255)
                     ggtggGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGaaaaa                          <  1:4835409/51‑1 (MQ=255)
                     ggtggGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGaaaaa                          <  1:758907/51‑1 (MQ=255)
                      gtggGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAAt                         >  1:5718224/1‑51 (MQ=255)
                       tggGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATg                        <  1:1138442/51‑1 (MQ=255)
                       tggGTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATg                        >  1:5602824/1‑51 (MQ=255)
                         ggTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGaa                      <  1:4213635/51‑1 (MQ=255)
                         ggTAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGaa                      <  1:1895688/51‑1 (MQ=255)
                           tAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAAt                    >  1:5314204/1‑51 (MQ=255)
                           tAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAAt                    <  1:1268114/51‑1 (MQ=255)
                           tAAAAACGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAAt                    <  1:5717454/51‑1 (MQ=255)
                            aaaaaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATc                   >  1:4153723/1‑51 (MQ=255)
                            aaaaaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATc                   <  1:2218681/51‑1 (MQ=255)
                             aaaaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCa                  >  1:2635445/1‑51 (MQ=255)
                             aaaaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCa                  <  1:4711421/51‑1 (MQ=255)
                             aaaaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCa                  <  1:3982049/51‑1 (MQ=255)
                              aaaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAg                 <  1:5126852/51‑1 (MQ=255)
                               aaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgc                >  1:2232222/1‑51 (MQ=255)
                               aaCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgc                >  1:2212508/1‑51 (MQ=255)
                                aCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgcg               <  1:3119328/51‑1 (MQ=255)
                                aCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgcg               <  1:2535837/51‑1 (MQ=255)
                                aCGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgcg               <  1:3114149/51‑1 (MQ=255)
                                 cGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgcgc              >  1:974054/1‑51 (MQ=255)
                                 cGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgcgc              >  1:1315411/1‑51 (MQ=255)
                                 cGGCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgcgc              <  1:239929/51‑1 (MQ=255)
                                  ggCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGcc             <  1:2663185/51‑1 (MQ=255)
                                  ggCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGcc             <  1:2587435/51‑1 (MQ=255)
                                   gCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAgc                >  1:3650061/1‑47 (MQ=255)
                                   gCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCg            <  1:5712976/51‑1 (MQ=255)
                                   gCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCg            <  1:3378528/51‑1 (MQ=255)
                                   gCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCg            <  1:5995554/51‑1 (MQ=255)
                                   gCTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCg            >  1:1773987/1‑51 (MQ=255)
                                    cTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGa           >  1:1600889/1‑51 (MQ=255)
                                    cTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGa           <  1:4304101/51‑1 (MQ=255)
                                    cTGTGGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGa           <  1:3981403/51‑1 (MQ=255)
                                     tgtgGTAAATCTACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAc          >  1:666067/1‑51 (MQ=255)
                                     tgtgGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAc          >  1:2918993/1‑51 (MQ=255)
                                      gtgGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAcg         >  1:5967297/1‑51 (MQ=255)
                                      gtgGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAcg         <  1:72098/51‑1 (MQ=255)
                                       tgGTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAcgg        <  1:2456047/51‑1 (MQ=255)
                                         gTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAcggcg      <  1:4585280/51‑1 (MQ=255)
                                         gTAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAcggcg      <  1:6392709/51‑1 (MQ=255)
                                          tAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAcggcgg     >  1:6318965/1‑51 (MQ=255)
                                          tAAATCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGAcggcgg     <  1:3572918/51‑1 (MQ=255)
                                            aaTCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGACGGCGGCa   >  1:5020505/1‑51 (MQ=255)
                                            aaTCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGACGGCGGCa   <  1:4162884/51‑1 (MQ=255)
                                             aTCGACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGACGGCGGCAg  >  1:3138625/1‑51 (MQ=255)
                                                |                                               
CCCCGGGCAGAAAGTCGGCCTGGTGGGTAAAAACGGCTGTGGTAAATCTACCCTGCTGGCATTGCTGAAAAATGAAATCAGCGCCGACGGCGGCAG  >  REL606/3409268‑3409363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: