Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,529,781	G→A	8.3%	S72S (TCG→TCA)	PP_4018 →	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,529,781	0	G	A	8.3%	137.2 / 2.4	48	S72S (TCG→TCA)	PP_4018	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base A (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

AGCCGTGTTCGTCGCCGAGCAGCACATTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTA  >  NC_002947/4529634‑4529928                                                                                                                                                     |                                                                                                                                                    aGCCGTGTTCGTCGCCGAGCAGCACATTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCgg                                                                                                                                                    >  1:211532/1‑149 (MQ=255)   ccGTGATGGGCGCCGAGCAGCACAACCCGCCAGAGCTGGAGATCGATTCGGCAGACCACGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTg                                                                                                                                                  <  7:141374/149‑1 (MQ=255)        ttCGTCGCCGAGTAGCACATTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTc                                                                                                                                               >  7:77131/1‑147 (MQ=255)            tcgCCGAGCAGCACATTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGAc                                                                                                                                         >  3:49780/1‑149 (MQ=255)                   gcagcaCATTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTg                                                                                                                                  <  6:51041/149‑1 (MQ=255)                   gcagcaCATTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTg                                                                                                                                  <  3:96630/149‑1 (MQ=255)                       cacaTTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGAc                                                                                                                                         <  8:310614/138‑1 (MQ=255)                       cacaTTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGAc                                                                                                                                         >  7:310614/1‑138 (MQ=255)                       cacaTTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCt                                                                                                                              <  1:46215/149‑1 (MQ=255)                          aTTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGaa                                                                                                                           <  4:36675/149‑1 (MQ=255)                           ttCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAg                                                                                                                          >  4:310209/1‑149 (MQ=255)                           ttCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCAGTGCTCAAGGACTGGCGTGGCCTGAAg                                                                                                                          >  8:5255/1‑149 (MQ=255)                           ttCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCAGTGCTCAAGGACTGGCGTGGCCTGAAg                                                                                                                          >  8:342/1‑149 (MQ=255)                               gCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTgg                                                                                                                      <  4:49780/149‑1 (MQ=255)                                ccAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTggg                                                                                                                     <  2:211532/149‑1 (MQ=255)                                 cAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTGTGAACGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGACCAATCGGGCGCAAATCGGTGCTCAAGGAATGGCGTGGCCTGAAGGTggg                                                                                                                     >  1:89036/1‑149 (MQ=255)                                 aagagCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGc                                                                                                                                     <  6:137632/131‑1 (MQ=255)                                     gCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATg                                                                                                                <  5:171802/149‑1 (MQ=255)                                         gAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGGAGGTGGGCGATGCGtt                                                                                                            >  1:171708/1‑149 (MQ=255)                                         gAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTACTGGCCGACGGCACCATCGGGCACATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGt                                                                                                             <  2:304668/148‑1 (MQ=255)                                          aGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTg                                                                                                           <  5:26413/149‑1 (MQ=255)                                          aGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTg                                                                                                           <  3:310209/149‑1 (MQ=255)                                          aGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCAGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTg                                                                                                           <  7:342/149‑1 (MQ=255)                                          aGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCAGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTg                                                                                                           <  7:5255/149‑1 (MQ=255)                                               gATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCaa                                                                                                      >  6:32334/1‑149 (MQ=255)                                                        gACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCAt                                                                                             <  8:77131/149‑1 (MQ=255)                                                           cAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCAtcgt                                                                                          >  5:159127/1‑149 (MQ=255)                                                                   ccTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCgg                                                                                                 >  2:42106/1‑134 (MQ=255)                                                                        aCTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAAc                                                                             >  8:83457/1‑149 (MQ=255)                                                                          ttCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACcg                                                                           >  7:17575/1‑149 (MQ=255)                                                                             cTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGa                                                                        >  4:267071/1‑149 (MQ=255)                                                                               ggCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGAc                                                                       <  5:110416/148‑1 (MQ=255)                                                                                  ccTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGATGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTg                                                                    >  7:195730/1‑148 (MQ=255)                                                                                       gAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGaaa                                                                 <  8:17575/147‑1 (MQ=255)                                                                                       gAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAAca                                                               <  7:83457/149‑1 (MQ=255)                                                                                       gAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTAGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAAc                                                                <  1:304432/148‑1 (MQ=255)                                                                                           gCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAAcagcag                                                           <  2:240588/149‑1 (MQ=255)                                                                                                tATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGa                                                          <  8:192181/145‑1 (MQ=255)                                                                                                tATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGa                                                          >  7:192181/1‑145 (MQ=255)                                                                                                          gTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAg                                            >  3:161369/1‑149 (MQ=255)                                                                                                           tACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAgg                                           <  6:219012/149‑1 (MQ=255)                                                                                                                 gcgcCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGcacg                                     >  2:57013/1‑149 (MQ=255)                                                                                                                      tgctgGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGcacgccacg                                >  7:122148/1‑149 (MQ=255)                                                                                                                       gctgGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGcacgccacgc                               >  1:105539/1‑149 (MQ=255)                                                                                                                              gACGGCACCATCGGGCGCATCTCGTTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTa                        <  8:195730/149‑1 (MQ=255)                                                                                                                                  gCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAg                    <  2:105539/149‑1 (MQ=255)                                                                                                                                      cATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGc                >  2:109051/1‑149 (MQ=255)                                                                                                                                      cATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGc                >  4:63810/1‑149 (MQ=255)                                                                                                                                        tCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTg              <  3:63810/149‑1 (MQ=255)                                                                                                                                                     ggTGCTCAAGGACTGGCGTGGCCTGGAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTa  <  2:171708/148‑1 (MQ=255)                                                                                                                                                     ggTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTcgc     <  4:161369/145‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGCCGTGTTCGTCGCCGAGCAGCACATTCCGCCAGAGCTGGAGTTCGATTCGGAAGACCAGGACGCCCTGCACTTCCTGGCCCTG‑GAAGGCGACTATCCGATCGGTACTGCGCGCCTGCTGGCCGACGGCACCATCGGGCGCATCTCGGTGCTCAAGGACTGGCGTGGCCTGAAGGTGGGCGATGCGTTGGTCAATGCGGTCATCGTCGAAGCGCAGAACCGCGACCTGAAACAGCAGATGCTCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTA  >  NC_002947/4529634‑4529928

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