Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	797,015	T→C	8.9%	Y10H (TAT→CAT)	PP_0682 →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	797,015	0	T	C	8.9%	120.3 / 3.1	45	Y10H (TAT→CAT)	PP_0682	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base C (2/2);  total (24/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 = 9.96e-01

TGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATC  >  NC_002947/796886‑797158                                                                                                                                  |                                                                                                                                                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                      <  3:197550/148‑1 (MQ=255)                                                                                                          gAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAg                      >  4:197550/1‑148 (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)                                                                                                                                  |                                                                                                                                                TGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTTTCCAGGAAGTGGTCCACGTCCTGCCCAATGTCTTTGCGTTGAGCGAAGCTGACCTGATCCTGGTGATC  >  NC_002947/796886‑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: ‑