Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,854	G→T	9.5%	G115V (GGC→GTC)	PP_4562 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,854	0	G	T	9.5%	154.3 / 5.0	52	G115V (GGC→GTC)	PP_4562	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base T (2/3);  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 = 7.39e-01

AGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCT  >  NC_002947/5182707‑5183001                                                                                                                                                      |                                                                                                                                                     aGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                                     <  5:9065/149‑1 (MQ=255) aGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                                     <  7:163742/149‑1 (MQ=255)  gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                                     >  4:362869/1‑148 (MQ=255)  gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                                    <  5:112467/149‑1 (MQ=255)  gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                                    <  5:107108/149‑1 (MQ=255)  gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                                    <  7:52364/149‑1 (MQ=255)    tGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTCGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑ag                                                                                                                                                 >  5:121753/1‑149 (MQ=255)     gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCGCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCCGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCg‑cgc                                                                                                                                                <  6:268373/149‑1 (MQ=255)      cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCG‑cgca                                                                                                                                               <  6:38399/149‑1 (MQ=255)                 ggTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCT‑AGCAGCTGGACTTcc                                                                                                                                    <  4:326024/149‑1 (MQ=255)                      ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCg‑cgc                                                                                                                                                <  4:301503/132‑1 (MQ=255)                      ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCg‑cgc                                                                                                                                                >  3:301503/1‑132 (MQ=255)                       caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACt                                                                                                                              >  4:2269/1‑149 (MQ=255)                         ccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGAGCCGACCGGGACTCCCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGCCTTCCGCAACTTc                                                                                                                            >  3:324566/1‑149 (MQ=255)                          caATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCg                                                                                                                           >  7:287107/1‑149 (MQ=255)                                ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAgg                                                                                                                     >  7:97911/1‑149 (MQ=255)                                ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAg                                                                                                                      <  3:50989/148‑1 (MQ=255)                                ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAg                                                                                                                      <  6:217587/148‑1 (MQ=255)                                 cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGc                                                                                                                    <  3:362869/149‑1 (MQ=255)                                       gCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAgg                                                                                                              <  2:125814/149‑1 (MQ=255)                                         tCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGcc                                                                                                            >  4:6588/1‑149 (MQ=255)                                                gAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                               <  2:221934/139‑1 (MQ=255)                                                gAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                               >  1:221934/1‑139 (MQ=255)                                                     aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtct                                                                                                <  3:2269/149‑1 (MQ=255)                                                     aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtct                                                                                                >  2:87915/1‑149 (MQ=255)                                                        ggAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGa                                                                                             >  5:360979/1‑149 (MQ=255)                                                         gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                            <  6:360979/149‑1 (MQ=255)                                                         gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                            <  6:121753/149‑1 (MQ=255)                                                         gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                            <  2:47904/149‑1 (MQ=255)                                                         gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCCGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                            <  1:65188/149‑1 (MQ=255)                                                             tGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGa                                                                                        <  8:209790/149‑1 (MQ=255)                                                                       tGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGa                                                                              >  7:31779/1‑149 (MQ=255)                                                                         gCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGAt                                                                             <  1:87915/148‑1 (MQ=255)                                                                          cAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGa                                                                           <  8:287107/149‑1 (MQ=255)                                                                               gccgcTCCACCATCAGACTGGCCCTGATCGCCGACCGTGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATc                                                                      <  4:324566/149‑1 (MQ=255)                                                                                        ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGc                                                             <  6:211966/149‑1 (MQ=255)                                                                                                       tCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGcc                                               >  5:113925/1‑148 (MQ=255)                                                                                                               gACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCctt                                      >  3:100654/1‑149 (MQ=255)                                                                                                                        gtccagctgctCGCCACCCGGACGGGTG‑‑TCG‑CGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGtt                             >  3:192586/6‑149 (MQ=38)                                                                                                                          gctgctGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                           >  8:62475/1‑149 (MQ=255)                                                                                                                          gctgctGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                           >  6:381283/1‑149 (MQ=255)                                                                                                                          ccagctgctCGCCACCCGGACGGGTGTCGCGCAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGt                              >  3:242308/4‑149 (MQ=38)                                                                                                                             gctgcGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCt                        >  3:137367/1‑149 (MQ=255)                                                                                                                               tgcGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg                      >  3:112162/1‑149 (MQ=255)                                                                                                                               tgcGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTAGAACAACCCGCCAAGGGCCTTCTTGAAGTTGACCtt                       >  7:126322/1‑148 (MQ=255)                                                                                                                                  gcgACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTg                   <  4:137367/149‑1 (MQ=255)                                                                                                                                         ccGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                           >  7:265293/1‑134 (MQ=255)                                                                                                                                         ccGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                           <  8:265293/134‑1 (MQ=255)                                                                                                                                         ccGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGCCTATCTCGATGCGATCGTCATGGATGACGATCAGCCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTg             >  2:145153/1‑148 (MQ=255)                                                                                                                                             ccGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAg        >  5:212389/1‑149 (MQ=255)                                                                                                                                                 gTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg    >  4:295928/1‑149 (MQ=255)                                                                                                                                                 gTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg    <  8:97911/149‑1 (MQ=255)                                                                                                                                                   g‑‑gCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCt  >  7:265500/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCT  >  NC_002947/5182707‑5183001

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