Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,834,428	(G)6→5	coding (110/2319 nt)	yicI ←	putative alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,834,423	0	G	.	98.6%	316.2 / ‑2.7	70	coding (115/2319 nt)	yicI	putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (38/31);  minor base A (1/0);  total (39/31)
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

GCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGT  >  NC_000913/3834327‑3834518                                                                                                 |                                                                                                gCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACgggg                                                                                              <  1:263952/100‑1 (MQ=255)  cTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTCCGTTCACGCACATCACggggg                                                                                             >  1:747795/1‑100 (MQ=255)  cTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACggggg                                                                                             >  1:995807/1‑100 (MQ=255)  cTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACggggg                                                                                             >  1:305514/1‑100 (MQ=255)  cTCAATCCGCACACCGACAATACCTTACTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACggggg                                                                                             >  1:546584/1‑100 (MQ=255)    cAATCCCCACACCGACAATAACTACCTGTGGGGAGAAAAAGCGCAACGTAAATAACAGCGTATCAATATGCCAGGTACGTTCACACACACCACggaggc                                                                                            >  1:366303/1‑98 (MQ=255)      aTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACgggggcagc                                                                                         <  1:945264/100‑5 (MQ=255)      aTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACgggggcagc                                                                                         <  1:986423/100‑5 (MQ=255)       tCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGgcagca                                                                                       <  2:1487930/100‑1 (MQ=255)        ccGCACACCGACAATACCTTCCTGTGGGGGGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGGCAGGTACGTTCACGCACATCAC‑GGGGgcagca                                                                                       >  1:825716/1‑99 (MQ=255)         cGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATa                                                                                     >  1:180314/1‑100 (MQ=255)         cGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATa                                                                                     >  2:1045240/1‑100 (MQ=255)         cGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑‑GGGGCAGCATAg                                                                                    >  2:1134354/1‑100 (MQ=255)            acacCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGac                                                                                  >  1:819818/1‑100 (MQ=255)             cacCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGacc                                                                                 <  1:256325/100‑1 (MQ=255)               ccGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGaccac                                                                               >  1:381071/1‑100 (MQ=255)                cGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGaccacc                                                                              >  1:1705911/1‑100 (MQ=255)                  aCAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCAt                                                                            <  2:1599187/100‑1 (MQ=255)                  aCAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCAt                                                                            <  2:161047/100‑1 (MQ=255)                    aaTACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCAttt                                                                          <  1:610364/100‑1 (MQ=255)                     aTACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTc                                                                         >  2:1430615/1‑100 (MQ=255)                      tACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCAttt                                                                          >  2:412514/1‑98 (MQ=255)                       cccctacTGTAGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATGAAGCTGCCAGGTACGTTCACGCACATGAC‑GGGGGCAGCAGAGACGACCATTTCAt                                                                       >  1:1517171/7‑100 (MQ=255)                        ccTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCAtt                                                                      >  1:855724/1‑100 (MQ=255)                         cTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTa                                                                     <  2:1455704/100‑1 (MQ=255)                             cTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCct                                                                 <  1:390078/100‑1 (MQ=255)                               gtgGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCctgc                                                               <  2:1799708/100‑1 (MQ=255)                               gtgGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCctgc                                                               <  2:625911/100‑1 (MQ=255)                                 ggggAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCctgctg                                                             <  1:176924/100‑1 (MQ=255)                                  gggAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGt                                                            >  1:1778130/1‑100 (MQ=255)                                     agaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCa                                                         <  2:747907/100‑1 (MQ=255)                                     agaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCa                                                         <  1:92590/100‑1 (MQ=255)                                     agaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGACAACATAAAACACCATTTCATTATCCTGCTGttaa                                                         >  2:306333/1‑98 (MQ=255)                                       aaaaaGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCA‑AGGGGGCAGCATAGACCACCATTTCATTATACTGCTGTTCaa                                                        >  1:227334/1‑99 (MQ=255)                                        aaaaGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAAcc                                                      >  1:820127/1‑100 (MQ=255)                                        aaaaGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAAc                                                       >  2:1537737/1‑99 (MQ=255)                                           aGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCg                                                   <  1:94355/100‑1 (MQ=255)                                            gcgcAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGa                                                  >  2:510844/1‑100 (MQ=255)                                             cgcAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGaa                                                 >  2:271912/1‑100 (MQ=255)                                              gcAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAac                                                <  2:414787/100‑1 (MQ=255)                                               cAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAaca                                               >  1:660643/1‑100 (MQ=255)                                                 aCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACAcc                                             <  2:600572/100‑1 (MQ=255)                                                 aCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACAcc                                             <  1:1593672/100‑1 (MQ=255)                                                   gTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTg                                           <  2:1655797/100‑1 (MQ=255)                                                    taaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGa                                          >  2:286975/1‑100 (MQ=255)                                                    taaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGa                                          <  1:1230318/100‑1 (MQ=255)                                                    taaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGa                                          >  2:1472082/1‑100 (MQ=255)                                                        taaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCg                                      >  2:191273/1‑100 (MQ=255)                                                         aaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCgg                                     <  1:164448/100‑1 (MQ=255)                                                           aGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGt                                   >  2:1470284/1‑100 (MQ=255)                                                           aGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGt                                   <  1:1430604/100‑1 (MQ=255)                                                            ggCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTg                                  <  2:723233/100‑1 (MQ=255)                                                            ggCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTg                                  >  1:29600/1‑100 (MQ=255)                                                             gCGTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACAACTGAAGCAGgaga                                 >  2:769786/1‑100 (MQ=255)                                                               gTATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAAt                               <  2:1669260/100‑1 (MQ=255)                                                                tATCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATc                              >  1:574643/1‑100 (MQ=255)                                                                 aTCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCa                             >  1:1426570/1‑100 (MQ=255)                                                                  tCAAGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCaa                            >  1:933405/1‑100 (MQ=255)                                                                    aaGCTGCCAGGTACTTTCACGTACATCAC‑GGGGGTAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCaaa                           <  2:652083/99‑1 (MQ=255)                                                                    aaGCTGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAAt                          >  2:1489046/1‑100 (MQ=255)                                                                        tGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAg                      <  2:1426581/100‑1 (MQ=255)                                                                        tGCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAg                      <  2:381073/100‑1 (MQ=255)                                                                         gCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAgg                     <  2:1478781/100‑1 (MQ=255)                                                                         gCCAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAgg                     <  1:1622245/100‑1 (MQ=255)                                                                          ccAGGTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGc                    >  1:1343115/1‑100 (MQ=255)                                                                              gTACGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCACGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAgg                >  2:114409/1‑100 (MQ=255)                                                                                 cGTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATGTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTg             >  1:971891/1‑100 (MQ=255)                                                                                  gTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGa            >  2:458931/1‑100 (MQ=255)                                                                                  gTTCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGAGGGTGAATCAAATTGAGGCCAGGTTGa            >  2:405805/1‑100 (MQ=255)                                                                                   ttCACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACTTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGaa           <  1:271911/100‑1 (MQ=255)                                                                                     cACGCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAACCGGGTGAATCAAATTGAGACCAGGTTGa            >  1:552179/1‑97 (MQ=255)                                                                                        gCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAAc      <  2:350325/100‑1 (MQ=255)                                                                                        gCACATCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAAc      <  2:1053925/100‑1 (MQ=255)                                                                                         cacaTCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAAcc     >  1:699011/1‑100 (MQ=255)                                                                                         cacaTCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAAcc     >  1:879070/1‑100 (MQ=255)                                                                                           caTCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAg   >  2:469382/1‑100 (MQ=255)                                                                                           caTCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAg   >  2:298784/1‑100 (MQ=255)                                                                                            aTCAC‑GGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGt  <  1:484471/100‑1 (MQ=255)                                                                                                 |                                                                                                GCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGT  >  NC_000913/3834327‑3834518

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