Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 801,414 A→G 53.1% S186S (AGT→AGC dhaL ← dihydroxyacetone kinase, C‑terminal domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported801,4140AG53.1% 7.6 / 65.0 49S186S (AGT→AGCdhaLdihydroxyacetone kinase, C‑terminal domain
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (12/14);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGCTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGCG  >  Exported/801348‑801480
                                                                  |                                                                  
tGCGGCTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACt                                                                  <  1:640557/69‑1 (MQ=255)
  cGGCTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGt                                                              <  1:3227647/71‑1 (MQ=255)
  cGGCTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGt                                                              <  1:1038425/71‑1 (MQ=255)
  cGGCTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGt                                                              <  1:1426204/71‑1 (MQ=255)
     cTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCa                                                           >  1:502354/1‑71 (MQ=255)
         cGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCga                                                       <  1:3312800/71‑1 (MQ=255)
         cGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCga                                                       <  1:1541288/71‑1 (MQ=255)
         cGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCga                                                       <  1:308969/71‑1 (MQ=255)
          gCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCAcc                                                         >  1:1413652/1‑68 (MQ=255)
            cAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCgaga                                                     >  1:2173700/1‑70 (MQ=255)
              aCATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGCCCAATGCTGCGTTCACCGAGATaa                                                  >  1:2285922/1‑71 (MQ=255)
              aCATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGt                                                              >  1:130966/1‑59 (MQ=255)
              aCATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATaa                                                  <  1:1197867/71‑1 (MQ=255)
               cATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGtt                                                             <  1:2730717/59‑1 (MQ=255)
                 tttGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTg                                               >  1:282632/1‑71 (MQ=255)
                  ttGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATa                                                                    >  1:3119515/1‑49 (MQ=255)
                  ttGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTg                                               >  1:2542898/1‑70 (MQ=255)
                        catAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGc                                        >  1:2657089/1‑71 (MQ=255)
                        catAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGc                                        >  1:1877207/1‑71 (MQ=255)
                           aaaCATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCtt                                     >  1:2590730/1‑71 (MQ=255)
                            aaCATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                    <  1:2054963/71‑1 (MQ=255)
                            aaCATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                    <  1:1372414/71‑1 (MQ=255)
                            aaCATCACCGAGGTCGCACCGGGATCCTGGTAACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                    <  1:1387929/71‑1 (MQ=255)
                                 cACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGc                                                                <  1:3665787/38‑1 (MQ=255)
                                 cACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGc                                                                <  1:110968/38‑1 (MQ=255)
                                  aCCGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGc                              <  1:1403273/71‑1 (MQ=255)
                                   ccGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCgag                                                      <  1:1144559/46‑1 (MQ=255)
                                   ccGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCgag                                                      <  1:3598511/46‑1 (MQ=255)
                                   ccGAGGTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCgag                                                      <  1:3329800/46‑1 (MQ=255)
                                    cGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCtt                            <  1:37911/71‑1 (MQ=255)
                                        gTCGCGCCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCat                        >  1:2793871/1‑71 (MQ=255)
                                           gcgcCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAAc                                                 >  1:262879/1‑43 (MQ=255)
                                           gcgcCGGGATCCTGGTGACCAATGCTGCGTTCACCGAGATAAc                                                 >  1:2375547/1‑43 (MQ=255)
                                               cGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgtaatc                 >  1:2097357/1‑71 (MQ=255)
                                                 ggATCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcgtaatcg                >  1:278947/1‑70 (MQ=255)
                                                    tCCTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACt            <  1:485381/71‑1 (MQ=255)
                                                      cTGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACttt          <  1:536259/71‑1 (MQ=255)
                                                       tGGTGACCAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTtg         >  1:2701589/1‑71 (MQ=255)
                                                        ggTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatcg                      >  1:897953/1‑57 (MQ=255)
                                                        ggTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCatc                       >  1:1994328/1‑56 (MQ=255)
                                                        ggTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACt            <  1:3546951/67‑1 (MQ=255)
                                                           gACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTtgtg       >  1:3163649/1‑69 (MQ=255)
                                                           gACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTtgtg       >  1:3032907/1‑69 (MQ=255)
                                                             ccAATGCTGCGTTCACCGAGATAACTGGCGCGGCCttt                                    <  1:680294/38‑1 (MQ=255)
                                                              cAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTg                           >  1:2807847/1‑46 (MQ=255)
                                                              cAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTg                           >  1:384798/1‑46 (MQ=255)
                                                              cAATGCTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGCg  >  1:3605620/1‑71 (MQ=255)
                                                              cAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGCg  <  1:3060173/71‑1 (MQ=255)
                                                              cAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGCg  <  1:2371419/71‑1 (MQ=255)
                                                               aaTACTGCGTTCACCGAGATAACTGGCGCGGCCtt                                     <  1:438752/35‑1 (MQ=255)
                                                                  |                                                                  
TGCGGCTAACGCCAACATTTGCATCATAAACATCACCGAGGTCGCGCCGGGATCCTGGTGACCAATACTGCGTTCACCGAGATAACTGGCGCGGCCTTTGCGGGCTTGCATCGTAATCGTACTTTGTGCAGCG  >  Exported/801348‑801480

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: