Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,532,120:1	+C	100%	coding (110/762 nt)	yeaE ←	methylglyoxal reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,532,120	1	.	C	88.1%	121.7 / 11.4	42	coding (110/762 nt)	yeaE	methylglyoxal reductase
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base C (21/16);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

AAGTCGATGCAATCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑CCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCT  >  NC_002947/3531973‑3532254                                                                                                                                                      |                                                                                                                                        aaGTCGATGCAATCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑c                                                                                                                                        >  2:302248/1‑149 (MQ=255)   gTCGATGCAATCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑ccc                                                                                                                                      <  4:94605/149‑1 (MQ=255)   gTCGATGCAATCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑ccc                                                                                                                                      >  6:180892/1‑149 (MQ=255)    tCGATGCAATCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑cccc                                                                                                                                     >  2:356183/1‑149 (MQ=255)            aTCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑cccccccccttc                                                                                                                             <  5:103537/149‑6 (MQ=255)             tCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTGGGCCCAACGACCTCTTCTG‑ccccccccttcg                                                                                                                             <  1:204964/149‑6 (MQ=255)                   cccAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACa                                                                                                                      >  2:80736/1‑149 (MQ=255)                     cAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACAtc                                                                                                                    >  7:117430/1‑149 (MQ=255)                             tGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑CCCCCCC‑TTCGGCATACATCTCGGCGgtg                                                                                                           >  3:95546/1‑149 (MQ=255)                               aGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGgtg                                                                                                           <  1:80736/148‑1 (MQ=255)                                 gCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCg                                                                                                        <  5:180892/149‑1 (MQ=255)                                  cTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑CCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGAt                                                                                                      <  1:356183/149‑1 (MQ=255)                                        tCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGt                                                                                                 >  7:232390/1‑149 (MQ=255)                                            aGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAg                                                                                             >  1:249469/1‑149 (MQ=255)                                             ggCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAgg                                                                                            <  2:263616/149‑1 (MQ=255)                                              gCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGc                                                                                           <  2:249469/149‑1 (MQ=255)                                                     ccATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCt                                                                                     <  7:93186/148‑1 (MQ=255)                                                     ccATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCt                                                                                     >  8:93186/1‑148 (MQ=255)                                                     ccATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTc                                                                                    >  2:269595/1‑149 (MQ=255)                                                          ccGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATg                                                                                <  8:179461/148‑1 (MQ=255)                                                                  gCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTg                                                                        >  6:371444/1‑148 (MQ=255)                                                                    tGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑CCCCCCCGTTCGGCATACATCTCGGGGGTGTCGAGAAGTTTCAGGCCCAGCGCGGTGCCCTGTTGCAg                                                                     >  6:321434/1‑149 (MQ=255)                                                                           gtgCGGGTAGACATTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCATGCCCAGCTCGATTCCCTGTTGCAGcgcctc                                                               <  8:232390/148‑3 (MQ=255)                                                                              cGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAgg                                                                                                  >  7:108254/1‑110 (MQ=255)                                                                              cGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAgg                                                                                                  <  8:108254/110‑1 (MQ=255)                                                                               gggTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCa                                                                                              >  7:258626/1‑113 (MQ=255)                                                                               gggTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCa                                                                                              <  8:258626/113‑1 (MQ=255)                                                                               gggTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑CCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTc                                                         >  2:109364/1‑149 (MQ=255)                                                                                    gACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑CCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAg                                                    >  3:129632/1‑149 (MQ=255)                                                                                      ccTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCTGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGc                                                   <  8:117430/149‑1 (MQ=255)                                                                                        ttGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTc                                                         <  8:240894/141‑1 (MQ=255)                                                                                        ttGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTc                                                         >  7:240894/1‑141 (MQ=255)                                                                                            tGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGc                                              >  8:23157/1‑148 (MQ=255)                                                                                            tGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCg                                             >  1:96536/1‑149 (MQ=255)                                                                                                 aGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATg                                                                                >  1:200630/1‑109 (MQ=255)                                                                                                 aGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATg                                                                                <  2:200630/109‑1 (MQ=255)                                                                                                  gAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCGTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTAAGGCCCCGCGCGAGGCCCTGTTGCCGCGCCGCGCCTTCGGCAGCCGTGCGACCTgg                                       >  8:161041/1‑149 (MQ=255)                                                                                                         ccGGGCGCCGCGACCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCg                                <  1:78725/149‑1 (MQ=255)                                                                                                         ccGGGCGCCGCGACCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCg                                >  2:78725/1‑149 (MQ=255)                                                                                                                cggcgCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATg                          >  7:51858/1‑148 (MQ=255)                                                                                                                      ccGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAgg                   >  6:22355/1‑149 (MQ=255)                                                                                                                       cGGCAATGGCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGt                  <  8:51858/149‑1 (MQ=255)                                                                                                                               gCTTGG‑CCAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGc          <  1:269595/149‑1 (MQ=255)                                                                                                                                      ccAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATATCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTAGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATg     >  6:236274/1‑148 (MQ=255)                                                                                                                                       cAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCa                                                                      <  4:189103/82‑1 (MQ=255)                                                                                                                                       cAACGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCa                                                                      >  3:189103/1‑82 (MQ=255)                                                                                                                                        aaCGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCt  >  6:166752/1‑149 (MQ=255)                                                                                                                                        aaCGACCTCTTCTGCCCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCt  >  6:161512/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                        AAGTCGATGCAATCGGTGCCCAGCCGAGTGAGGCTGCGCTCGCAGGCTGCCGCCATGCCGCGCAGGCTGGCATTGTGCGGGTAGACCTTGCTGACCAGAAATACCCGGTCGCGGCGCCCGGCAATGGCTTGG‑CCAACGACCTCTTCTG‑CCCCCCC‑TTCGGCATACATCTCGGCGGTGTCGATCAGGTTCAGGCCCAGCTCGATGCCCTGTTGCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCT  >  NC_002947/3531973‑3532254

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