Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,054,020 C→A 14.5% S107S (TCC→TCA yihX → phosphatase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,054,0200C→A14.5% 4.9 ‑3.9 0.056321yihXphosphatase
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (21/26):  total (21/35)
Fisher's exact test strand distribution p-value = 1.94e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CTGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCAATACCAACCGCCTGCATACCACCTTCTGGCCGG  >  REL606/4053985‑4054054
                                   |                                  
cttcGTGAGCAGGGGCATCGCGTGGTGGTGCTTTcc                                    <  1:4939081/33‑1 (MQ=255)
cTGCGTGAGCAGGGGCATCTCGTGGTGGTGCTTTcc                                    <  1:474377/36‑1 (MQ=255)
cTGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTcc                                    <  1:3932155/36‑1 (MQ=255)
cTGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTcc                                    <  1:4677378/36‑1 (MQ=255)
cTGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTcc                                    <  1:6325836/36‑1 (MQ=255)
 tGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCa                                   <  1:1480557/36‑1 (MQ=255)
 tGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCa                                   <  1:1959931/36‑1 (MQ=255)
 tGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCa                                   <  1:2201057/36‑1 (MQ=255)
  gCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCaa                                  <  1:4784203/36‑1 (MQ=255)
  gCGTGAGCAGGGGCATCGCGTGGTGGTGCTTACCaa                                  <  1:2809999/36‑1 (MQ=255)
   cGTGAGCCGGGGCATCGCGTGGTGGTGCTTTCCAAt                                 <  1:5326901/36‑1 (MQ=255)
   cGTGAGCATGGGCATCGCGTGGTGGTGCTTTCCAAt                                 <  1:6390611/36‑1 (MQ=255)
   cGTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCAAt                                 <  1:5983058/36‑1 (MQ=255)
   cGTGAGCAGGGGCATAGCGTGGTGGTGCTTTCCAAt                                 >  1:7317753/1‑36 (MQ=255)
    gTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCAATa                                <  1:3438983/36‑1 (MQ=255)
    gTGAGAAGGGGCATCGAGTGGTGGTGATTTCCAAtc                                >  1:6805374/1‑35 (MQ=21)
     tGAGCATGGGCATCGCGTGGTGGTGCTTTCCAATAc                               <  1:6645947/36‑1 (MQ=255)
     tGAGCAGGGGCATCGCGTGGTGGTGCTTTCCAATAc                               <  1:5828871/36‑1 (MQ=255)
      gAGCAGGGGCATCGCGTGGTGGTGCTTTCCAATAcc                              <  1:4425177/36‑1 (MQ=255)
      gAGCAGGGGCATCGCGTGGTGGTGCTTTCCAATAcc                              <  1:2557813/36‑1 (MQ=255)
      gAGCAGGGGCATCGCGTGGTGGTGCTTTCCAATAcc                              >  1:7812057/1‑36 (MQ=255)
       aGCAGGGGCATCGCGTGGTGGTGCTTTCCAATACCa                             <  1:47279/36‑1 (MQ=255)
       aGCAGGGGCATCGCGTGGTGGTGCTTTCCAATACCa                             <  1:7217743/36‑1 (MQ=255)
         cAGGGGCATCGCGTGGTGGTGCTTTCCAATACCAAc                           <  1:6890435/36‑1 (MQ=255)
         cAGGGGCATCGCGTGGTGGTGCTTTCCAATACCAAc                           >  1:4223448/1‑36 (MQ=255)
          aGGGGCATCGCGTGGTGGTGCTTTCCAATACCAAcc                          >  1:5411245/1‑36 (MQ=255)
          aGGGGCATCGCGTGGTGGTGCTTTCCAATACCAAcc                          >  1:1047432/1‑36 (MQ=255)
          aGGGGCATCGCGTGGTGGTGCTTTCCAATACCAAcc                          <  1:3277511/36‑1 (MQ=255)
           ggggCATCGCGTGGTGGTGCTTTCCAATACCAACCg                         <  1:427956/36‑1 (MQ=255)
           ggggCATCGCGTGGTGGTGCTTTCCAATACCAACCg                         <  1:3958740/36‑1 (MQ=255)
            gggCATCGCGTGGTGGTGCTTTCCAATACCAACCGc                        >  1:772028/1‑36 (MQ=255)
             ggCATCGCGTGGTGGTGCTTTCCAATACCAACCGcc                       <  1:918956/36‑1 (MQ=255)
             ggCATCGCGTGGTGGTGCTTTCCAATACCAACCGcc                       <  1:2743091/36‑1 (MQ=255)
             ggCATCGCGTGGTGGTGCTTTCCAATACCAACCGcc                       >  1:5834285/1‑36 (MQ=255)
              gCATCGCGTGGTGGTGCTTTCCAATACCAACCGCCt                      <  1:1211824/36‑1 (MQ=255)
              gCATCGCGTGGTGGTGCTTTCCAATACCAACCGCCt                      <  1:7105172/36‑1 (MQ=255)
               cATCGCGTGGTGGTGCTTTCCAATACCAACCGCCTg                     <  1:276020/36‑1 (MQ=255)
                aTCGCGTGGTGGTGCTTTCCAATACCAACCGCCTGc                    >  1:3983153/1‑36 (MQ=255)
                aTCGCGTGGTGGTGCTTTCCAATACCAACCGCCTGc                    >  1:5802055/1‑36 (MQ=255)
                  cgcgTGGTGGTCCTTCCAAAAACCAACCGCCTGCAt                  <  1:1980706/36‑1 (MQ=16)
                   gcgTGGTGGTGCTTTCCAATACCAACCGCCCGCATa                 >  1:2176858/1‑36 (MQ=255)
                   gcgGGGTGGTGTTTTCCAATCCCAACCGCCTGCATa                 <  1:6824787/36‑1 (MQ=25)
                    cgTGGTGGTGTTTTCAAAAACAAACCGCCTGCATac                <  1:5928138/36‑1 (MQ=16)
                    cgTGGTGGTGCTTTCCAATACCAACCGCCTGCATac                >  1:6923079/1‑36 (MQ=255)
                     gtggtgttgCTTTCAAATACAACCCGCCTGCATacc               <  1:6733104/36‑1 (MQ=19)
                     gtggtggtgTTTCCAAATACCAACCGCCTGCATacc               <  1:1697183/36‑1 (MQ=25)
                     gtggtggtgCTTTCCAATACCAACCGCCTGCATacc               >  1:4258969/1‑36 (MQ=255)
                     gtggtggtgCTTTCCAATACCAACCGCCTGCATacc               >  1:2576250/1‑36 (MQ=255)
                     gtggtggtgCTTTCAAATACCAACCGCCTGCATacc               <  1:5936567/36‑1 (MQ=255)
                       ggtggtgCTTTCCAATACCAACCGCCTGCATaccac             >  1:5503987/1‑36 (MQ=255)
                       gggggtgCTTTCAAAAACAAACCGCCTGCATaccac             <  1:2482717/33‑1 (MQ=16)
                        gtggtgCTTTCAAATACCAACCGCCTGCATaccacc            <  1:109528/36‑1 (MQ=255)
                        gtggtgCTTTCAAAACCCACCCGCCTGCATaccacc            <  1:3846941/36‑1 (MQ=19)
                        gtggtaCTTACCAACACCAACCGCCTGCATaccacc            <  1:1374055/36‑1 (MQ=25)
                        ggggtgCTTCCAAATACCAACCGCCTGCATaccacc            <  1:6881944/34‑1 (MQ=37)
                           gtgCTTTCCAATACCAACCGCCTGCATACCACCTTc         >  1:7037029/1‑36 (MQ=255)
                             gCTTTCCAATACCAACCGCCTGCATACCACCTTCTg       >  1:4105215/1‑36 (MQ=255)
                              cTTTCCAATACCAACCGCCTGCATACCACCTTCTgg      >  1:2596758/1‑36 (MQ=255)
                               tttCCAATACCAACCGCCTGCATACCACCTTCTGGc     >  1:6420598/1‑36 (MQ=255)
                               tttCCAATACCAACCGCCTGCATACCACCTTCTGGc     >  1:2985558/1‑36 (MQ=255)
                                 tCCAATACCAACCGCCTGCATACCACCTTCTGGCCg   >  1:2115793/1‑36 (MQ=255)
                                  ccAATACCAACCGCCTGCATACCACCTTCTGGCCgg  >  1:2133835/1‑36 (MQ=255)
                                   |                                  
CTGCGTGAGCAGGGGCATCGCGTGGTGGTGCTTTCCAATACCAACCGCCTGCATACCACCTTCTGGCCGG  >  REL606/4053985‑4054054

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