Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	797,018	G→A	8.7%	A11T (GCC→ACC) ‡	PP_0682 →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	797,018	0	G	A	8.7%	129.4 / 2.4	46	A11T (GCC→ACC) ‡	PP_0682	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base A (2/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

TGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATCCTGTCA  >  NC_002947/796886‑797164                                                                                                                                     |                                                                                                                                                   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                                                                                                     <  1:96611/95‑1 (MQ=255)                                                                                      gatgaGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGc                                                                                                     >  2:96611/1‑95 (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                            >  4:197550/1‑148 (MQ=255)                                                                                                          gAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAg                            <  3:197550/148‑1 (MQ=255)                                                                                                           aaCGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTcc                                                             <  6:207522/114‑1 (MQ=255)                                                                                                           aaCGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTcc                                                             >  5:207522/1‑114 (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)                                                                                                                                   aTGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATCCTGTCa  >  6:170489/1‑149 (MQ=255)                                                                                                                                     gCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTc                                                        >  3:171176/1‑93 (MQ=255)                                                                                                                                     gCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTc                                                        <  4:171176/93‑1 (MQ=255)                                                                                                                                     |                                                                                                                                                   TGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATCCTGTCA  >  NC_002947/796886‑797164

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