Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,050,073	C→T	8.0%	S1712S (AGC→AGT)	yeeJ →	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,050,073	0	C	T	8.0%	130.8 / 2.4	50	S1712S (AGC→AGT)	yeeJ	putative adhesin
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base T (2/2);  total (24/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 = 2.49e-01

GCAAAACTAGCCTCCTTAACCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGCCAGTGCAGAT  >  NC_000913/2049927‑2050214                                                                                                                                                   |                                                                                                                                               gCAAAACTAGCCTCCTTAACCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAagcag                                                                                                                                              >  1:497592/1‑149 (MQ=255)         aGCCTCCTTAACCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTGGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACAAGCg                                                                                                                                      >  1:615199/1‑149 (MQ=255)             tcctTAACCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGa                                                                                                                                  >  1:734190/1‑149 (MQ=255)             tcctTAACCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGa                                                                                                                                  >  1:711018/1‑149 (MQ=255)             tcctTAACCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGa                                                                                                                                  >  2:608054/1‑149 (MQ=255)                   aCCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGga                                                                                                                            <  1:600329/149‑1 (MQ=255)                         gTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑Cg                                                                                                                     >  2:464515/1‑149 (MQ=255)                               tCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTACCGGCAACCCGGTAGAAGGCCTAACAGTTAATTTCCGCGGAACCTACCTCACGCTAAGCAGAACCAGCGTTTAAACGGGTGATGCGGGGtt                                                                                                                >  2:591344/1‑149 (MQ=255)                               tCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGTAGCACCAGCGTTGAAACGGATGAT‑CGGGGttt                                                                                                               >  2:770586/1‑149 (MQ=255)                               tCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGTAGCACCAGCGTTGAAACGGATGAT‑CGGGGttt                                                                                                               >  2:574800/1‑149 (MQ=255)                                  ttCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGc                                                                                                            <  1:808803/149‑1 (MQ=255)                                       cgtcAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGaaa                                                                                                       >  1:187182/1‑149 (MQ=255)                                         tcAGCACGACCGGGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAAtt                                                                                                     <  2:580159/149‑1 (MQ=255)                                             cACGACAGAGGGCGCAACCATAACTGCAATCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTtg                                                                                                 <  2:428284/149‑1 (MQ=255)                                                    gAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCAGGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATACTTGTGACAAg                                                                                          >  2:730996/1‑149 (MQ=255)                                                    gAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAg                                                                                          >  2:709394/1‑149 (MQ=255)                                                    gAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAg                                                                                          <  1:313136/149‑1 (MQ=255)                                                       ggCGCAACCATGGCGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAGCCTCCGTCACGCTAAGCAGCACCAGCGTGGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCAc                                                                                       <  2:5136/149‑1 (MQ=255)                                                           cAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTCTCGCTGAAATTCTTGTGACAAGCCCCGAg                                                                                   <  2:754007/149‑1 (MQ=255)                                                           cAACCATGACGGCAAGCGTCACTGACGCCAACGGCAATGCGGTAGAAGGCATAAAAGTTAATTTCCGCGGACCCTCCGGCACGCTAAGCAGCACCTGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAg                                                                                   <  2:376347/149‑1 (MQ=255)                                                           cAACCAGGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAg                                                                                   <  1:217435/149‑1 (MQ=255)                                                                 tGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGa                                                                             <  2:311877/149‑1 (MQ=255)                                                                           cGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAg                                                                   <  2:191714/149‑1 (MQ=255)                                                                             tCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCGCGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGAAACGAGGTCGGACTGAAGACGGtt                                                                 >  1:601857/1‑149 (MQ=255)                                                                               aCTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTc                                                               >  2:558901/1‑149 (MQ=255)                                                                                    cGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCct                                                          <  1:216934/149‑1 (MQ=255)                                                                                     gCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGGGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGGAAACAGTTTCAGGctc                                                         >  2:446096/1‑149 (MQ=255)                                                                                      cTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTCCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCctct                                                        <  1:398067/149‑1 (MQ=255)                                                                                        aaCGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCctctcc                                                      <  2:653056/149‑2 (MQ=255)                                                                                           ggCAACCAGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGc                                                   <  2:853365/149‑1 (MQ=255)                                                                                              aaCCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGCTGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGa                                                >  1:238115/1‑149 (MQ=255)                                                                                              aaCCCGGTAGAAGGAATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTAGCAGa                                                >  1:159069/1‑149 (MQ=255)                                                                                               aCCCGGTAGAAGGCATAAAAGTTAATTTGCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGAt                                               <  2:212462/149‑1 (MQ=255)                                                                                                cccGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGTAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCTATGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATa                                              <  1:574800/149‑1 (MQ=255)                                                                                                cccGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGTAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATa                                              <  1:770586/149‑1 (MQ=255)                                                                                                        aaGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGACATTCTTGTGACAAGCACCGAGGTCGGACGGAAAACAGTTTCAGCCTCTCTGGCAGATCAAccccct                                      >  1:46510/1‑145 (MQ=255)                                                                                                             aTAAAAGTTAATTCCCGCGGATCCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGt                                 <  2:615199/149‑1 (MQ=255)                                                                                                              tAAAAGTTAATTTCCGCGGAACCTCAGTCACGCTAAGCAGCACCAGCGTTGAAACGGATTAT‑CAGGGTTTCGCTGAAATTCTTGTGACAAGCCCCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTc                                >  1:481194/1‑149 (MQ=255)                                                                                                               aaaaGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGGTTCGCTGCAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGAGAAACCTACTGAAGTCa                               >  2:593778/1‑149 (MQ=255)                                                                                                                aaaGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCAt                              <  2:497592/149‑1 (MQ=255)                                                                                                                      aaTTTTCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCGGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCAtcctgcc                        >  2:678380/1‑144 (MQ=255)                                                                                                                          tcggcgGAAGCTCCGTCACGCTAAGCAGAACCAGCGTTGAAACGGATGAT‑CGGGGTTACGGTGAAATTCTTGTGACAAGCACCGAGGTCGGGCTGCAAACCGTTTCAGCCTCTCTGACCGATACACCTACGGAGGTCTCCTcgcgggtg                    >  2:483753/4‑145 (MQ=255)                                                                                                                          tCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTa                    <  1:217043/149‑1 (MQ=255)                                                                                                                              cgGAACCTCGGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGa                <  1:524717/149‑1 (MQ=255)                                                                                                                                   ccTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGcc           >  1:823616/1‑149 (MQ=255)                                                                                                                                     tCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGGAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCACTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGACAg         >  2:711412/1‑149 (MQ=255)                                                                                                                                      ccGTCACACTAATCACCACCAGCGTTGAAACGGATGAT‑CTGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGCCAGt        <  1:619490/149‑1 (MQ=255)                                                                                                                                       cGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGCCAGTg       <  2:469182/149‑1 (MQ=255)                                                                                                                                          cACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGCCAGTGCAg    <  1:709394/149‑1 (MQ=255)                                                                                                                                            cGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGCCAGTGCAGaa  >  2:207730/1‑148 (MQ=255)                                                                                                                                            cGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGCCAGTGCAGAt  <  2:328623/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GCAAAACTAGCCTCCTTAACCTCTGTATACTCATTCGTCGTCAGCACGACCGAGGGCGCAACCATGACGGCAAGCGTCACTGACGCTAACGGCAACCCGGTAGAAGGCATAAAAGTTAATTTCCGCGGAACCTCCGTCACGCTAAGCAGCACCAGCGTTGAAACGGATGAT‑CGGGGTTTCGCTGAAATTCTTGTGACAAGCACCGAGGTCGGACTGAAAACAGTTTCAGCCTCTCTGGCAGATAAACCTACTGAAGTCATCTCGCGATTACTGAATGCCAGTGCAGAT  >  NC_000913/2049927‑2050214

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