Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	307,865:1	+C	100%	intergenic (+880/‑989)	hslO → / → PP_0254	molecular chaperone HslO/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	307,865	1	.	C	98.1%	177.3 / ‑3.2	52	intergenic (+880/‑989)	hslO/PP_0254	molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (25/26);  total (25/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 = 8.04e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAG‑CGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAG  >  NC_002947/307722‑308004                                                                                                                                                   |                                                                                                                                            cATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGGACTACCCTTTCGTCAGCCGAAAG‑GGta                                                                                                                                         <  1:17924/149‑2 (MQ=255)  aTGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAG‑ccgtta                                                                                                                                       >  6:256337/1‑144 (MQ=255)   tGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTAc                                                                                                                                       <  7:305133/149‑1 (MQ=255)       ctcGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGa                                                                                                                                   >  3:90007/1‑149 (MQ=255)          ggTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCg                                                                                                                                >  1:63407/1‑149 (MQ=255)            tGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGt                                                                                                                              <  5:256337/149‑1 (MQ=255)            tGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGt                                                                                                                              <  3:18184/149‑1 (MQ=255)              cAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTga                                                                                                                            <  7:196814/149‑1 (MQ=255)                    ttCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGAcgacga                                                                                                                      <  6:59390/149‑1 (MQ=255)                       ctgctgCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCa                                                                                                                   >  6:58733/1‑149 (MQ=255)                           tgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGt                                                                                                               >  1:352449/1‑149 (MQ=255)                                  tgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGc                                                                                                         >  2:154420/1‑148 (MQ=255)                                    ccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGaa                                                                                                      <  2:63407/149‑1 (MQ=255)                                         gTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGt                                                                                                 >  3:29515/1‑149 (MQ=255)                                         gTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGt                                                                                                 >  3:162245/1‑149 (MQ=255)                                          tACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                <  4:90007/149‑1 (MQ=255)                                              gCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCt                                                                                            >  2:270407/1‑149 (MQ=255)                                               ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGAGCCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCt                                                                                            >  1:282022/1‑149 (MQ=255)                                                cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCtt                                                                                           >  3:236393/1‑148 (MQ=255)                                                cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTc                                                                                          <  3:52718/149‑1 (MQ=255)                                                     cACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACa                                                                                      <  2:282022/148‑1 (MQ=255)                                                     cACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACAt                                                                                     <  5:58733/149‑1 (MQ=255)                                                       cTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATgg                                                                                   <  2:178525/149‑1 (MQ=255)                                                       cTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATgg                                                                                   <  2:51556/149‑1 (MQ=255)                                                        tGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATgg                                                                                   >  6:264742/1‑148 (MQ=255)                                                             ggCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcg                                                                             <  7:311814/149‑1 (MQ=255)                                                              gCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcgg                                                                            >  4:366279/1‑149 (MQ=255)                                                                cAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgct                                                                          <  3:366279/149‑1 (MQ=255)                                                                 aaCATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                >  2:334216/1‑126 (MQ=255)                                                                 aaCATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctg                                                                         >  2:126411/1‑149 (MQ=255)                                                                 aaCATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgctg                                                                         <  2:352449/149‑1 (MQ=255)                                                                 aaCATCGGCGAAGAAGGCGATGTGA‑CCTTGATCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt                                                                                                <  1:334216/126‑1 (MQ=255)                                                                            agaagGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAg                                                               <  1:2609/148‑1 (MQ=255)                                                                                   gATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGa                                                       <  5:42656/149‑1 (MQ=255)                                                                                    aTGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                      <  4:236393/149‑1 (MQ=255)                                                                                      gtgA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCt                                                    >  8:227511/1‑149 (MQ=255)                                                                                        gA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGt                                                  >  4:118544/1‑149 (MQ=255)                                                                                         a‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTc                                                 >  8:213512/1‑149 (MQ=255)                                                                                                  cttcGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                      <  4:69531/136‑1 (MQ=255)                                                                                                  cttcGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                      >  3:69531/1‑136 (MQ=255)                                                                                                  cttcGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGagaag                                          <  4:162245/148‑1 (MQ=255)                                                                                                  cttcGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGagaag                                          <  4:29515/148‑1 (MQ=255)                                                                                                            tCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCgg                                <  3:46953/148‑1 (MQ=255)                                                                                                            tCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGt                               >  2:377316/1‑149 (MQ=255)                                                                                                              cGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCa                             >  7:320008/1‑149 (MQ=255)                                                                                                                    tGGCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTgg                        >  1:327535/1‑148 (MQ=255)                                                                                                                      gCAA‑GACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGaaa                     <  1:270407/149‑1 (MQ=255)                                                                                                                              tACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATc               >  7:340204/1‑148 (MQ=255)                                                                                                                                cccTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCt                          <  3:189453/135‑1 (MQ=255)                                                                                                                                cccTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCt                          >  4:189453/1‑135 (MQ=255)                                                                                                                                cccTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCaa             <  8:320008/148‑1 (MQ=255)                                                                                                                                  cTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCa          >  7:149715/1‑149 (MQ=255)                                                                                                                                          tGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAg  >  2:148754/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                            CATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGA‑CCTTGTTCTTCGGCCTGTCCGGCACTGGCAA‑GACTACCCTTTCGGCTGACGAAAG‑CGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAG  >  NC_002947/307722‑308004

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: ‑