Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	618,312	C→T	10.1%	G86E (GGG→GAG)	PP_0532 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	618,312	0	C	T	10.1%	134.5 / 3.3	49	G86E (GGG→GAG)	PP_0532	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/24);  new base T (2/3);  total (22/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.00e+00

CGCTGCTGCACCTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTG  >  NC_002947/618165‑618460                                                                                                                                                      |                                                                                                                                                     cGCTGCTGCACCTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTccc                                                                                                                                                     >  1:126533/1‑149 (MQ=255)      ctgcACCTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAgcg                                                                                                                                                >  7:56226/1‑149 (MQ=255)        gcACCTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCAGTCCAGTTGCTGCTGATCCAGCGTCCCCAgcg                                                                                                                                                >  3:345262/1‑147 (MQ=255)        gcACCTGCTCCGGCAGCGCCCGGGTGTAGTCCAGCGTGCGCCCGGCCCTGTCGATGCGCAGTTCGGACCGGTAACCAACC‑GCCACCAGCACGG‑GCCCCTGGCCCACCGCGAGCCAGTCCAGTAGCTGCTGATCCAGCGTCCCCAccg                                                                                                                                                >  3:339782/1‑147 (MQ=255)         cACCTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCgcag                                                                                                                                             >  3:50639/1‑149 (MQ=255)         cACCTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCgca                                                                                                                                              <  8:272597/148‑1 (MQ=255)           ccTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCgcagca                                                                                                                                           >  7:208361/1‑149 (MQ=255)               cTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCTCCAGCGCAGCATTcc                                                                                                                                       >  2:12094/1‑149 (MQ=255)               cTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCTCCAGCGCAGCATTcc                                                                                                                                       >  2:17705/1‑149 (MQ=255)                 ccGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTcccc                                                                                                                                                    <  6:139099/134‑1 (MQ=255)                 ccGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTcccc                                                                                                                                                    >  5:139099/1‑134 (MQ=255)                   ggCAGCGCCGGGGTGTAGTCCAGCGTGCGGCGGGCCGTGTCGAGGCGCAGGTCGCTCCGGTAACCCACC‑GCGACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATc                                                                                                                                   <  4:339782/149‑1 (MQ=255)                   ggCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCacc                                                                                                                                   <  4:158506/149‑3 (MQ=255)                   ggCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATc                                                                                                                                   <  8:56226/149‑1 (MQ=255)                   ggCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATc                                                                                                                                   <  4:345262/149‑1 (MQ=255)                            ggggTGTAGTCCAGCGTGCGGCCGGCCCGGTCGATGCGCAGGTCGCTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCAGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCca                                                                                                                          <  8:252176/149‑1 (MQ=255)                              ggTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCcaca                                                                                                                        <  7:266409/149‑1 (MQ=255)                                     tCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGGTGCGCCTGGCCCACGGCGAGCCTGTCAAGTTGCTGCTGATCAAGCGTCCCCAGCGCAGCATTCCCAACGCGCATCCACAGCACCg                                                                                                                  >  1:286465/1‑149 (MQ=255)                                       cAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCggg                                                                                                                >  2:56140/1‑148 (MQ=255)                                              cggccggcCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGa                                                                                                        <  1:279536/149‑1 (MQ=255)                                               ggccggccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGAt                                                                                                       >  2:109998/1‑149 (MQ=255)                                                 ccggccCTGTCGATGCGCAGTTCGGTCCGGTATCCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGc                                                                                                     >  1:243759/1‑149 (MQ=255)                                                      ccTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCTCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGgcg                                                                                                <  8:170463/149‑1 (MQ=255)                                                        tGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCt                                                                                              >  6:164868/1‑149 (MQ=255)                                                          gcctagcgcAGTTCGGTCCCGTAACCCCCCGGCCACC‑CCCCGG‑GCGGCGGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTgg                                                                                            <  8:140480/144‑1 (MQ=255)                                                              tGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTccc                                                                                                                                                     <  5:83734/88‑1 (MQ=255)                                                              tGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTccc                                                                                                                                                     >  6:83734/1‑88 (MQ=255)                                                                       cGGTCCGGTATCCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCtt                                                                                <  2:243759/148‑1 (MQ=255)                                                                        ggTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTg                                                                               <  2:126533/148‑1 (MQ=255)                                                                        ggTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTg                                                                               <  8:207534/148‑1 (MQ=255)                                                                        ggTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTACAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGc                                                                              <  8:330456/149‑1 (MQ=255)                                                                           ccGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCaa                                                                           <  6:69220/149‑1 (MQ=255)                                                                              gTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAAcc                                                                         >  7:172892/1‑148 (MQ=255)                                                                                         gCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCgg                                                                >  2:55218/1‑147 (MQ=255)                                                                                         gCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTa                                                              <  1:55218/149‑1 (MQ=255)                                                                                         gCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTa                                                              <  2:286465/149‑1 (MQ=255)                                                                                             ccaGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCTCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGa                                                          <  1:17705/149‑1 (MQ=255)                                                                                             ccaGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCTCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGa                                                          <  1:12094/149‑1 (MQ=255)                                                                                               agcagGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACg                                                        <  4:50639/149‑1 (MQ=255)                                                                                                gcagGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGc                                                       >  5:268957/1‑149 (MQ=255)                                                                                                 cagGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATGGAGCATCCACATCACGGGGTACCCGGATGCTGGCGCTGGACGATTCGGCTTGCCAAccc                                                                        <  3:163994/131‑1 (MQ=255)                                                                                                  agGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAg                                                     >  7:168554/1‑149 (MQ=255)                                                                                                    gg‑gCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGACTTGCCAACCCACAGGCGGTAGTGACGCAGGc                                                   >  6:43146/1‑149 (MQ=255)                                                                                                          ccTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCg                                              <  5:164868/149‑1 (MQ=255)                                                                                                           cTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCgg                                             >  5:281159/1‑149 (MQ=255)                                                                                                              gCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGcc                                           <  1:56140/148‑1 (MQ=255)                                                                                                                   cGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCTAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCAGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTcc                                     <  2:195194/149‑1 (MQ=255)                                                                                                                      cGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGa                                  <  6:268957/149‑1 (MQ=255)                                                                                                                               ccAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTgg                         <  1:319897/149‑1 (MQ=255)                                                                                                                                        gctgATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGa                >  7:7343/1‑149 (MQ=255)                                                                                                                                        gctgATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGa                >  7:30923/1‑149 (MQ=255)                                                                                                                                              ccAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATa              <  8:118801/145‑1 (MQ=255)                                                                                                                                              ccAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATa              >  7:118801/1‑145 (MQ=255)                                                                                                                                               cAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGc           >  1:337511/1‑147 (MQ=255)                                                                                                                                                      cccAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTg  >  5:134742/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CGCTGCTGCACCTGCTCCGGCAGCGCCGGGGTGTAGTCCAGCGTGCGGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACC‑GCCACCAGCAGGG‑GCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTG  >  NC_002947/618165‑618460

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