Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,531,764 A→G 8.9% C156R (TGC→CGC)  yeaE ← methylglyoxal reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,531,7640AG 8.9% 120.5 / 2.8 45C156R (TGC→CGC) yeaEmethylglyoxal reductase
Reads supporting (aligned to +/- strand):  ref base A (23/18);  new base G (2/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

GTCACCCAGGCCAGGCAGACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGCCCCAACGCTTGATCTTGCCT  >  NC_002947/3531628‑3531900
                                                                                                                                         |                                                                                                                                        
gTCACCCAGGCCAGGCAGACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGc                                                                                                                                >  2:12468/1‑147 (MQ=255)
 tCACCCAGGCCAGGCAGACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAAAAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAg                                                                                                                             >  6:200420/1‑149 (MQ=255)
    cccAGGCCAGGCAGACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTc                                                                                                                           >  5:341634/1‑148 (MQ=255)
    cccAGGCCAGGCAGACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTc                                                                                                                           >  3:274111/1‑148 (MQ=255)
     ccaggccaggCAGACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCg                                                                                                                          >  3:163328/1‑148 (MQ=255)
                 gACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCcgc                                                                                                             <  2:111965/149‑1 (MQ=255)
                    cTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGtt                                                                                                          >  3:280130/1‑149 (MQ=255)
                             tgtgGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACgg                                                                                                                                    >  2:88520/1‑114 (MQ=255)
                                    ccATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTacaa                                                                                          >  4:328405/1‑149 (MQ=255)
                                    ccATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTacaa                                                                                          <  4:202170/149‑1 (MQ=255)
                                          cgcTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtg                                                                                    <  3:327372/149‑1 (MQ=255)
                                             tCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtggtt                                                                                 >  2:212272/1‑149 (MQ=255)
                                              cGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtggtt                                                                                 >  1:87311/1‑149 (MQ=255)
                                               ggCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtggttgg                                                                               <  6:341634/149‑1 (MQ=255)
                                                gCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCGCCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtggttgg                                                                               >  1:174921/1‑148 (MQ=255)
                                                gCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCGCCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtggttgg                                                                               <  2:180539/148‑1 (MQ=255)
                                                gCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCGCCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtggttgg                                                                               <  2:174921/148‑1 (MQ=255)
                                                gCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCGCCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCtggttgg                                                                               >  1:180539/1‑148 (MQ=255)
                                                 cGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTtggtg                                                                             <  3:154937/149‑1 (MQ=255)
                                                    tttCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGca                                                                          <  3:138035/149‑1 (MQ=255)
                                                                  gggTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCTAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTtgtg                                                            <  5:307282/149‑1 (MQ=255)
                                                                  gggTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTtgtg                                                            <  1:12468/149‑1 (MQ=255)
                                                                          caacaaGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTACAGTTcgc                                                    >  8:138140/1‑149 (MQ=255)
                                                                              aaGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTtgtg                                                            >  2:36858/1‑137 (MQ=255)
                                                                              aaGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTtgtg                                                            <  1:36858/137‑1 (MQ=255)
                                                                               aGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAgg                                               >  7:287060/1‑149 (MQ=255)
                                                                                gcgcCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGgt                                              >  6:142236/1‑149 (MQ=255)
                                                                                   cccGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAAGTCGATGCCGCGTTGGGCCGGGTTGTAGAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGgtcgt                                           >  4:17086/1‑149 (MQ=255)
                                                                                    ccGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGgtcgtc                                          <  3:10754/149‑1 (MQ=255)
                                                                                         cTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCg                                                                                                                    >  4:197017/1‑70 (MQ=255)
                                                                                             gcCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGaa                                 <  4:163328/149‑1 (MQ=255)
                                                                                                 aGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGtt                              <  4:280130/149‑1 (MQ=255)
                                                                                                 aGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGtt                              <  2:87311/149‑1 (MQ=255)
                                                                                                   cGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGATCGCGCAGGTCGTCGACGTCGAAGTTgg                            <  1:212272/149‑1 (MQ=255)
                                                                                                    gggCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGa                           >  5:169533/1‑149 (MQ=255)
                                                                                                     ggCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGac                          <  5:94654/149‑1 (MQ=255)
                                                                                                      gCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTCGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGaca                         >  5:218355/1‑149 (MQ=255)
                                                                                                       cAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGacac                        <  4:274111/149‑1 (MQ=255)
                                                                                                        aGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACg                       >  8:198094/1‑149 (MQ=255)
                                                                                                          tAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGcc                     <  5:142236/149‑1 (MQ=255)
                                                                                                              ccATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGCCCCaa                 >  2:141757/1‑149 (MQ=255)
                                                                                                                       gCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGCCCCAACGCTTGAt         <  8:287060/148‑1 (MQ=255)
                                                                                                                             cACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGCCCCAAGGCTTGATcttgcct  >  4:269226/1‑149 (MQ=255)
                                                                                                                             cACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGCCCCAACGCTTGATcttgcct  >  2:102989/1‑149 (MQ=255)
                                                                                                                             cACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATAAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGCCCCAAAGCTTGATcttgcct  >  7:289797/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GTCACCCAGGCCAGGCAGACCTGGGCCGGTGTGGCGCCATGGCGCTCGGCGATTTCCGCCAGTACCGGGTGCTGCAACAAGCGCCCGGCCTGCGCCAGCGGGCAGTAGGCCATGGTTGGCAAGCCACGTTTTCCGCACCACGGTAGCAAGTCGAATTCGATGCCGCGTTGGGCCGGGTTGTACAACACCTGGTTGGTGGCACAATCAGGGTTGTGCAGTTCGCGCAGGTCGTCGACGTCGAAGTTGGACACGCCCCAACGCTTGATCTTGCCT  >  NC_002947/3531628‑3531900

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