Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	664,902	T→C	7.5%	S589P (TCG→CCG)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	664,902	0	T	C	7.5%	145.4 / 2.7	53	S589P (TCG→CCG)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (2/2);  total (26/27)
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.46e-01

GAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGACCGCGACGGT  >  NC_002947/664764‑665048                                                                                                                                            |                                                                                                                                                   gAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccgg                                                                                                                                         >  2:381849/1‑150 (MQ=255)   aaCCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccggcc                                                                                                                                       <  1:1260291/150‑1 (MQ=255)            gAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCg                                                                                                                              <  2:1026912/150‑1 (MQ=255)            gAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCg                                                                                                                              >  1:1283847/1‑150 (MQ=255)               ccAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCg                                                                                                                                                  >  1:1145460/1‑127 (MQ=255)               ccAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCg                                                                                                                                                  <  2:1145460/127‑1 (MQ=255)                cAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGACCGCCGTTAAGCCGACGC‑CCGCTATCCGCTGTAACTGACCTTGCAGGCCGTCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGgcgc                                                                                                                          >  2:144985/1‑150 (MQ=255)                  ggACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCt                                                                                                                        <  1:1173163/150‑1 (MQ=255)                        ggAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagca                                                                                                                  <  2:1051742/150‑1 (MQ=255)                          aaCTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCACACCATCGAAGT‑CGAAGGAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagcagc                                                                                                                <  2:1008757/150‑1 (MQ=255)                               gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGCGCCGGGTGATCGGCCGGCCACCCgactcct                                                                                                                        >  1:1008135/1‑131 (MQ=37)                               gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGCGCCGGGTGATCGGCCGGCCACCCgactcct                                                                                                                        <  2:1008135/137‑7 (MQ=37)                                   aaGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGgcgc                                                                                                                          <  1:1271642/131‑1 (MQ=255)                                   aaGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGgcgc                                                                                                                          >  2:1271642/1‑131 (MQ=255)                                      tCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCgg                                                                                                    >  1:891258/1‑150 (MQ=255)                                      tCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCgg                                                                                                    <  1:381849/150‑1 (MQ=255)                                        aCCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCt                                                                                                                        <  2:573587/128‑1 (MQ=255)                                        aCCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCt                                                                                                                        >  1:573587/1‑128 (MQ=255)                                        aCCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTg                                                                                                                       >  2:205427/1‑129 (MQ=255)                                        aCCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTg                                                                                                                       <  1:205427/129‑1 (MQ=255)                                        aCCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCt                                                                                                  >  1:361331/1‑150 (MQ=255)                                          ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGa                                                                                                <  2:891258/150‑1 (MQ=255)                                            ttGACCAGCCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAgc                                                                                              <  1:948628/150‑1 (MQ=255)                                               aCCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccggc                                                                                                                                        <  1:827700/105‑1 (MQ=255)                                               aCCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccggc                                                                                                                                        >  2:827700/1‑105 (MQ=255)                                                 cAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTg                                                                                         >  1:360988/1‑150 (MQ=255)                                                             gCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACGGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGa                                                                             <  1:694676/150‑1 (MQ=255)                                                                gtggtCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                          <  2:392764/150‑1 (MQ=255)                                                                   gtCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCGGCGGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTg                                                                       >  1:338036/1‑150 (MQ=255)                                                                                   gTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGgaga                                                       >  2:548967/1‑150 (MQ=255)                                                                                   gTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGCCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGgaga                                                       >  2:228996/1‑150 (MQ=255)                                                                                         ccGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATa                                                 <  2:360988/150‑1 (MQ=255)                                                                                         ccGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGACTGCCGCGCCTGCGAGGACGGTTCGAAGGGGAGGAGGCCggtg                                                 >  2:247744/1‑147 (MQ=255)                                                                                             aGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGCCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGa                                             <  1:228996/150‑1 (MQ=255)                                                                                                  cGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGACTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACggg                                         >  1:210577/1‑150 (MQ=255)                                                                                                       tcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                                      >  1:1090391/1‑132 (MQ=255)                                                                                                       tcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                                      <  2:1090391/132‑1 (MQ=255)                                                                                                          gCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGa                                 >  1:849859/1‑150 (MQ=255)                                                                                                              cAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTc                             <  2:849859/150‑1 (MQ=255)                                                                                                                  tGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGt                         <  1:780652/150‑1 (MQ=255)                                                                                                                    ctcctGCAGGCTGGCCTGCACGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTc                       >  2:725856/5‑150 (MQ=255)                                                                                                                     ccGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                          <  2:928548/98‑1 (MQ=255)                                                                                                                     ccGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                                          >  1:928548/1‑98 (MQ=255)                                                                                                                      cGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCt                     >  1:882282/1‑150 (MQ=255)                                                                                                                           aGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtggtgg                >  2:173488/1‑150 (MQ=255)                                                                                                                                  ccTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGACcg         >  1:120465/1‑150 (MQ=255)                                                                                                                                   cTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                                      <  2:636372/104‑1 (MQ=255)                                                                                                                                   cTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                                      >  1:636372/1‑104 (MQ=255)                                                                                                                                   cTGCAGGAGTCGGGTGGCCGACCGCACCCCCGGCACACAGACCAGTCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGACcgc        <  2:210577/150‑1 (MQ=255)                                                                                                                                     gCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGACCGCGa      <  1:173488/150‑1 (MQ=255)                                                                                                                                         gAGTCGAGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGACCGCGACGGt  <  1:125265/150‑1 (MQ=255)                                                                                                                                          aGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGa            <  2:195875/139‑1 (MQ=255)                                                                                                                                          aGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGa            >  1:195875/1‑139 (MQ=255)                                                                                                                                            |                                                                                                                                                   GAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGC‑CCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGGCCGACCGCGACGGT  >  NC_002947/664764‑665048

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