Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	587,076	C→G	10.0%	P43R (CCG→CGG)	PP_0497 →	Cro/CI family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	587,076	0	C	G	10.0%	140.7 / 6.3	50	P43R (CCG→CGG)	PP_0497	Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base G (3/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

GTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGT  >  NC_002947/586931‑587219                                                                                                                                                   |                                                                                                                                                gTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGAc                                                                                                                                              <  5:278649/149‑1 (MQ=255)      cTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCtt                                                                                                                                         >  5:173157/1‑149 (MQ=255)         gATTCCGCATATTAGAGATCACTTCGCTACAAACCTCAACCTTTCCTGCAGTCACTCCCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGAGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCttttc                                                                                                                                      <  6:249261/149‑3 (MQ=255)           ttCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                    <  6:50417/149‑1 (MQ=255)               gCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTg                                                                                                                                 >  3:303210/1‑148 (MQ=255)                 atatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAg                                                                                                                              >  6:263157/1‑149 (MQ=255)                   atGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCtttgccc                                                                                                                                    >  6:32595/1‑137 (MQ=255)                   atGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCtttgccc                                                                                                                                    <  5:32595/141‑5 (MQ=255)                    tGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGc                                                                                                                             <  1:242716/148‑1 (MQ=255)                               ttCGCTACAAGCCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGActt                                                                                                                >  6:172238/1‑149 (MQ=255)                               ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTcc                                                                                                                                                 <  6:262142/116‑1 (MQ=255)                               ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTcc                                                                                                                                                 >  5:262142/1‑116 (MQ=255)                               ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGa                                                                                                                   <  3:65093/146‑1 (MQ=255)                               ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGa                                                                                                                   >  4:65093/1‑146 (MQ=255)                               ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGAct                                                                                                                 >  7:24025/1‑148 (MQ=255)                                   cTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCTGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACttcttc                                                                                                            <  5:20615/149‑1 (MQ=255)                                    tACAAACCTCAAACTTGCCTGCAGTCACTACCTCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGc                                                                                                                             <  6:203621/131‑1 (MQ=255)                                    tACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGc                                                                                                                             >  5:203621/1‑131 (MQ=255)                                       aaaCCTCAAACTTGCCTGCAGTCACTACCGGTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGgtg                                                                                                        >  1:26643/1‑149 (MQ=255)                                            tCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑CACAAGTACCTGCGCGGGCAAAGCCGTCCGCCGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                   >  2:253756/1‑149 (MQ=255)                                              aaaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGa                                                                                                 <  8:24025/149‑1 (MQ=255)                                              aaaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGa                                                                                                 <  5:263157/149‑1 (MQ=255)                                               aaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGAt                                                                                                >  8:142176/1‑149 (MQ=255)                                                    gCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTAc                                                                                             <  2:26643/147‑1 (MQ=255)                                                     ccTGAAGTCACTACCGCTCTATTTCCTAAGTTTGCCGCCAGCTGTCTATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCTAGGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGaa                                                                                          <  1:248004/149‑1 (MQ=255)                                                      cTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAAc                                                                                         <  4:303210/149‑1 (MQ=255)                                                        gCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGaa                                                                                          >  4:153839/1‑146 (MQ=255)                                                        gCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGaa                                                                                          <  3:153839/146‑1 (MQ=255)                                                            tCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAAcc                                                                                   <  6:173157/149‑1 (MQ=255)                                                                     ctctATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCccgcc                                                                           >  4:275255/1‑148 (MQ=255)                                                                      tctATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGGGTCGAGGATTACGAACTGAACCTGCccgcc                                                                           >  6:82667/1‑147 (MQ=255)                                                                              cGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCg                                                                 >  4:276591/1‑149 (MQ=255)                                                                               gAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGc                                                                <  5:172238/149‑1 (MQ=255)                                                                                          ccAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGAtt                                                                                               >  6:172828/1‑107 (MQ=255)                                                                                          ccAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGAtt                                                                                               <  5:172828/107‑1 (MQ=255)                                                                                          ccAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCg                                                     <  7:142176/149‑1 (MQ=255)                                                                                                gTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGcc                                                               >  5:206204/1‑133 (MQ=255)                                                                                                gTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGcc                                                               <  6:206204/133‑1 (MQ=255)                                                                                                   gATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCggg                                            >  2:184423/1‑149 (MQ=255)                                                                                                     tCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTa                                          <  7:106342/149‑1 (MQ=255)                                                                                                                             tACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAg                   >  2:265060/1‑149 (MQ=255)                                                                                                                              aCCTGAGCGGGCAAAGCCGTCCGAAGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGACCGACGATCAGTTCGAAAGTCTGATCGGTTCGCGGGTATCGCAACAGGACCATCAAACAAGt                  >  4:82213/1‑149 (MQ=255)                                                                                                                               ccTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGt                  <  3:276591/148‑1 (MQ=255)                                                                                                                                cTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCTGGCCCAGCAACCCAGTGa                >  1:138949/1‑149 (MQ=255)                                                                                                                                cTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGa                >  1:256169/1‑149 (MQ=255)                                                                                                                                cTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGa                >  4:177184/1‑149 (MQ=255)                                                                                                                                  gAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc              >  1:58918/1‑149 (MQ=255)                                                                                                                                   aGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTAACGCCCCAGGCCCAGCAACCCAGTGAc               <  3:275255/147‑1 (MQ=255)                                                                                                                                      gttgaaaGCCGTCGGACCGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGAt          >  5:323743/5‑149 (MQ=255)                                                                                                                                             gCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAg   <  2:58918/149‑1 (MQ=255)                                                                                                                                              ccGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAgt  >  1:269035/1‑149 (MQ=255)                                                                                                                                              ccGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAgt  <  1:253756/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTC‑AACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGT  >  NC_002947/586931‑587219

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