Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	833,898	G→A	15.4%	A191T (GCA→ACA)	PP_0717 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	833,898	0	G	A	15.4%	64.3 / 4.4	26	A191T (GCA→ACA)	PP_0717	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (8/14);  new base A (2/2);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

CGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCT  >  NC_002947/833767‑834016                                                                                                                                    |                                                                                                                       gtgtggtgCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCAtt                                                                                                                    <  1:486084/134‑1 (MQ=255) gggtggtGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCATCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCAtt                                                                                                                    <  1:20460/135‑1 (MQ=255) cGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCAtt                                                                                                                    <  2:11784/136‑1 (MQ=255) cGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCAtt                                                                                                                    <  2:104743/136‑1 (MQ=255)           tACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCt                                                                                                          <  1:242694/136‑1 (MQ=255)           tACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCt                                                                                                          <  1:328174/136‑1 (MQ=255)                 gttCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAg                                                                                                    >  2:94994/1‑136 (MQ=255)                     cccGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACa                                                                                                >  1:472830/1‑136 (MQ=255)                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATg                                                                                              <  2:37442/136‑1 (MQ=255)                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGACAtt                                                                                                                    <  1:396103/114‑1 (MQ=255)                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGACAtt                                                                                                                    <  1:393695/114‑1 (MQ=255)                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGACAtt                                                                                                                    >  2:396103/1‑114 (MQ=255)                       cGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGACAtt                                                                                                                    >  2:393695/1‑114 (MQ=255)                             tGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTg                                                                                        >  2:10202/1‑136 (MQ=255)                                cAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTTAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATc                                                                                     <  1:433740/136‑1 (MQ=255)                                cAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATc                                                                                     <  1:367297/136‑1 (MQ=255)                                 aaTGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCa                                                                                    >  1:446086/1‑136 (MQ=255)                                        cAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGc                                                                             <  2:57750/136‑1 (MQ=255)                                        cAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGc                                                                             <  2:446086/136‑1 (MQ=255)                                                                                           cgcTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCg                          >  1:444461/1‑136 (MQ=255)                                                                                              tGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGa                       >  2:19589/1‑136 (MQ=255)                                                                                                 gTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCt                    <  1:94994/136‑1 (MQ=255)                                                                                                 gTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCt                    <  2:472830/136‑1 (MQ=255)                                                                                                     cACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATc                >  1:477070/1‑136 (MQ=255)                                                                                                                                gctgGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACcttgct  >  2:310021/1‑123 (MQ=255)                                                                                                                                gctgGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACcttgct  <  1:310021/123‑1 (MQ=255)                                                                                                                                    |                                                                                                                       CGGTGGTGCATACGTTGTTCCCCGAACCTGCCAATGCGCCAGCACCGCCGGTGCCGCCGCTGCTGGGTAGCGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCT  >  NC_002947/833767‑834016

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