Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	797,013	T→G	8.7%	M9R (ATG→AGG)	PP_0682 →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	797,013	0	T	G	8.7%	136.9 / 3.6	46	M9R (ATG→AGG)	PP_0682	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base G (2/2);  total (24/22)
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.84e-01

GCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATC  >  NC_002947/796866‑797158                                                                                                                                                    |                                                                                                                                                  gCTGTTGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGTCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGAt                                                                                                                                                   <  5:97848/148‑1 (MQ=255)  cTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATg                                                                                                                                                  <  2:188796/148‑1 (MQ=255)                     tgtCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGgctgc                                                                                                                               <  3:84007/148‑1 (MQ=255)                           gATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTGTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGAGGCATACATcgcgttttc                                                                                                                                  >  4:258483/1‑134 (MQ=255)                           gATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTCTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGAGGCATACATcgcgttttc                                                                                                                                  <  3:258483/139‑6 (MQ=255)                             tCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGa                                                                                                                      >  2:62067/1‑149 (MQ=255)                             tCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTAGCACCGa                                                                                                                      >  3:69576/1‑149 (MQ=255)                              cAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGAt                                                                                                                     >  7:203120/1‑149 (MQ=255)                                   tGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACt                                                                                                                >  6:162013/1‑149 (MQ=255)                                     cGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTc                                                                                                              <  4:69576/149‑1 (MQ=255)                                                 cGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTc                                                                                                  >  7:291004/1‑149 (MQ=255)                                                    gCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGc                                                                                               <  1:68875/149‑1 (MQ=255)                                                                aCAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTg                                                                                   <  5:30744/149‑1 (MQ=255)                                                                                aCTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTc                                                                    >  2:253631/1‑148 (MQ=255)                                                                                  tCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAg                                                                 >  1:198383/1‑149 (MQ=255)                                                                                  tCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAg                                                                 >  8:232396/1‑149 (MQ=255)                                                                                     ccGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGaa                                                              <  8:203120/149‑1 (MQ=255)                                                                                        tCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTg                                                           <  8:291004/149‑1 (MQ=255)                                                                                        tCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGAGATGTATGCCTCGCGATGGCTGATCGCACCGATCTACTTCGGCTTGTCCCTCGGCATGCTGGCCCTGGCGATGAAATTTTTCCAGGAAGTg                                                           >  7:147117/1‑149 (MQ=255)                                                                                          gCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCTCTGGCGCTGAAATTTTTCCAGGAAGTGGt                                                         >  5:297753/1‑149 (MQ=255)                                                                                          gCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGt                                                         >  2:145316/1‑149 (MQ=255)                                                                                           cTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGACGTggac                                                        <  7:232396/149‑3 (MQ=255)                                                                                               tttAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACg                                                    >  3:26297/1‑149 (MQ=255)                                                                                                  aGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTcc                                                 >  3:267855/1‑149 (MQ=255)                                                                                                   gTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTcc                                                 <  4:26297/148‑1 (MQ=255)                                                                                                    tCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGACGTGGACCACTTCCTg                                               <  5:154326/149‑1 (MQ=255)                                                                                                      aCCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGAGCTGAAATTTTTCCAGGAAGTGGTCGACGTCCTGcc                                             >  7:13216/1‑149 (MQ=255)                                                                                                          gatgaGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGc                                                                                               >  2:96611/1‑95 (MQ=255)                                                                                                          gatgaGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGc                                                                                               <  1:96611/95‑1 (MQ=255)                                                                                                                tCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCttt                                   >  8:140104/1‑149 (MQ=255)                                                                                                                 ccTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCTCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCttt                                   <  6:297753/148‑1 (MQ=255)                                                                                                                  cTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGc                                 <  5:162013/149‑1 (MQ=255)                                                                                                                  cTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGCCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGc                                 >  3:282974/1‑149 (MQ=255)                                                                                                                     aaCTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGtt                              <  2:198383/149‑1 (MQ=255)                                                                                                                       cTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGa                            <  2:264238/149‑1 (MQ=255)                                                                                                                           aTGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCg                         >  4:148953/1‑148 (MQ=255)                                                                                                                              gAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAg                      <  3:197550/148‑1 (MQ=255)                                                                                                                              gAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAg                      >  4:197550/1‑148 (MQ=255)                                                                                                                               aaCGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTcc                                                       >  5:207522/1‑114 (MQ=255)                                                                                                                               aaCGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTcc                                                       <  6:207522/114‑1 (MQ=255)                                                                                                                                       cTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGAt            >  6:263748/1‑149 (MQ=255)                                                                                                                                          gAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGATGAAATTTTTCCAGGAATTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGAtg           <  8:147117/147‑2 (MQ=255)                                                                                                                                            aaaCGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTgg       <  1:145316/149‑1 (MQ=255)                                                                                                                                              gcgcgAGGCATACATCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATcc          >  3:251237/2‑144 (MQ=255)                                                                                                                                              gcgcgAGGCATACATCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATcc          <  4:251237/143‑1 (MQ=255)                                                                                                                                              aCGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGt      >  3:303168/1‑148 (MQ=255)                                                                                                                                                 cgATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATc  <  4:267855/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATC  >  NC_002947/796866‑797158

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