Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,050,074	A→G	19.4%	S1713G (AGC→GGC)	yeeJ →	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,050,074	0	A	G	19.4%	140.8 / 29.8	72	S1713G (AGC→GGC)	yeeJ	putative adhesin
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base G (8/6);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

ACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAA  >  NC_000913/2049975‑2050172                                                                                                    |                                                                                                   aCCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAagca                                                                                                    <  2:397625/100‑1 (MQ=255)    gAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAgcagca                                                                                                 >  1:173011/1‑100 (MQ=255)        gcgcAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAg                                                                                             >  1:981160/1‑100 (MQ=255)         cgcAACCATGACGTCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGTAACCTCCGTCACGCTAAGCAGCAct                                                                                               <  2:1970366/97‑2 (MQ=255)         cgcAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGc                                                                                            >  2:2336046/1‑100 (MQ=255)           cAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGt                                                                                          >  1:1965458/1‑100 (MQ=255)            attcATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCg                                                                                           <  1:1467766/95‑1 (MQ=255)             aCCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTg                                                                                        >  1:2311548/1‑100 (MQ=255)             aCCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTg                                                                                        <  2:1789150/100‑1 (MQ=255)               cATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCTTCACGCTAAGCAGCACCAGCGTTGaa                                                                                      >  1:1933382/1‑100 (MQ=255)                  ccgatctaGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACg                                                                                   <  2:1164410/93‑1 (MQ=255)                   aCGGCAAGCGTCACTGGCGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACgg                                                                                  <  2:654895/100‑1 (MQ=255)                     ggCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCATCAGCGTTGAAACGgat                                                                                <  2:2361395/100‑1 (MQ=255)                     ggCAAGCGTCACTGACGCTAAAGGAAACACGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCAGTAACGCTAAGAAGAACCAGAGTTGAAACGgat                                                                                >  2:642333/1‑100 (MQ=255)                      gCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGgat                                                                                <  1:2263695/99‑1 (MQ=255)                         aGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACgagatcgg                                                                            >  1:1164421/1‑93 (MQ=255)                            gTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCggg                                                                         <  1:1039638/100‑1 (MQ=255)                             tCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCgggg                                                                        <  1:126411/100‑1 (MQ=255)                             tCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCgggg                                                                        <  1:472982/100‑1 (MQ=255)                              cACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGt                                                                       <  1:2079176/100‑1 (MQ=255)                                  gACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCg                                                                   >  2:2082388/1‑100 (MQ=255)                                   aCGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTAAAACGGATGATCGGGGTTTCGc                                                                  >  1:1021559/1‑100 (MQ=255)                                    cGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCt                                                                 <  1:2313883/100‑1 (MQ=255)                                    cGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCt                                                                 >  2:1982631/1‑100 (MQ=255)                                     gCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAGCCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCt                                                                 <  1:1049226/99‑1 (MQ=255)                                      cTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGa                                                               >  2:128223/1‑100 (MQ=255)                                       tAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGaa                                                              <  1:2276043/100‑1 (MQ=255)                                       tAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGaa                                                              <  2:535920/100‑1 (MQ=255)                                          cGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAAtt                                                           >  1:2030798/1‑100 (MQ=255)                                          cGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAAtt                                                           >  1:2087191/1‑100 (MQ=255)                                          cGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAAtt                                                           >  1:688744/1‑100 (MQ=255)                                           ggCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTc                                                          <  1:1861058/100‑1 (MQ=255)                                            gcatccatgTAGAAGGCATAAAAGTCAAGTTCTGCTGCACCTCAGTCACGCTAAGCAGCACCAGCGTTGAAACGCATCATCGGGGTTTCGCCGAAATCCt                                                         <  2:1183136/92‑1 (MQ=255)                                            gCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCt                                                         >  1:1367975/1‑100 (MQ=255)                                               aCCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTtgt                                                      >  1:1187267/1‑100 (MQ=255)                                               aCCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTtgt                                                      >  1:2084103/1‑100 (MQ=255)                                                cccGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTtgtg                                                     <  2:591110/100‑1 (MQ=255)                                                cccGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTtgtg                                                     <  1:1199001/100‑1 (MQ=255)                                                 ccGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTtgtg                                                     <  2:1789669/99‑1 (MQ=255)                                                 ccGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTtgtg                                                     <  1:901812/99‑1 (MQ=255)                                                  cGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGAc                                                   >  2:1363702/1‑100 (MQ=255)                                                  cGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGATTCGCTGAAATTCTTGTGAc                                                   >  2:76279/1‑100 (MQ=255)                                                  cGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGAc                                                   >  1:2317876/1‑100 (MQ=255)                                                   ggTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACa                                                  >  1:1157530/1‑100 (MQ=255)                                                    gTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTtgt                                                      <  1:33709/95‑1 (MQ=255)                                                     tAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAg                                                <  1:2114692/100‑1 (MQ=255)                                                      aGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGc                                               >  2:2016569/1‑100 (MQ=255)                                                          ggCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCg                                           >  2:864890/1‑100 (MQ=255)                                                           gtaTAAAAGTTTATTTTCGTGGAGCCTCCGTCACGCTAAGCAGTACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGc                                               <  1:705368/93‑1 (MQ=255)                                                           gCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGa                                          >  1:552146/1‑100 (MQ=255)                                                            cATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAg                                         >  1:1658347/1‑100 (MQ=255)                                                            cATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAg                                         <  1:1488805/100‑1 (MQ=255)                                                              tAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAgg                                        >  2:184647/1‑99 (MQ=255)                                                                 aaGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCgg                                    <  1:1518377/100‑1 (MQ=255)                                                                 aaGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCgg                                    >  1:1153137/1‑100 (MQ=255)                                                                 aaGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCgg                                    <  1:1669859/100‑1 (MQ=255)                                                                   gTTAATTTCCGCGGAACCTCCGTCACGCTAAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGAc                                  <  2:981149/100‑1 (MQ=255)                                                                   gTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGAc                                  >  1:540976/1‑100 (MQ=255)                                                                   gTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGAc                                  <  1:628522/100‑1 (MQ=255)                                                                             gcgGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGt                        >  2:433689/1‑100 (MQ=255)                                                                                 aaCCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCa                    <  1:864899/100‑1 (MQ=255)                                                                                  aCCTCCGTCACGCTAAGCAGCAACAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAg                   <  2:2119748/100‑1 (MQ=255)                                                                                   catCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAg                   <  2:1919920/97‑1 (MQ=255)                                                                                    cTCCGTCACGCTGAGCGGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCa                    <  2:140120/97‑1 (MQ=255)                                                                                    cTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGcc                 <  2:1153126/100‑1 (MQ=255)                                                                                    cTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGcc                 <  1:2082410/100‑1 (MQ=255)                                                                                    cTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGcc                 <  1:601064/100‑1 (MQ=255)                                                                                      ccGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCctc               <  2:1105548/100‑1 (MQ=255)                                                                                         tCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCctctct            >  2:1562478/1‑100 (MQ=255)                                                                                          cACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCctctccg           <  2:12933/100‑3 (MQ=255)                                                                                          cACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTg           >  2:2402363/1‑100 (MQ=255)                                                                                           aCGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTgg          >  1:741094/1‑100 (MQ=255)                                                                                           aCGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTgg          >  1:1102346/1‑100 (MQ=255)                                                                                              cTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAg       >  1:481037/1‑100 (MQ=255)                                                                                                 agcagcACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATa    >  1:737740/1‑100 (MQ=255)                                                                                                 agcagcACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATa    >  1:738185/1‑100 (MQ=255)                                                                                                   cagcaCCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATaaa  <  1:2337422/100‑1 (MQ=255)                                                                                                    |                                                                                                   ACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGATCGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAA  >  NC_000913/2049975‑2050172

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