Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,005,412	C→A	18.2%	A314A (GCC→GCA)	PP_5256 →	cyclic nucleotide‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,005,412	0	C	A	18.2%	57.0 / 3.8	22	A314A (GCC→GCA)	PP_5256	cyclic nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (3/15);  new base A (2/2);  total (5/17)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

AGTCGGTGGAGTACGAGGCCAGCGGCTTCGTGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGCGTTGAGCCAGCGCATGACTGCAGTCGAGTACCTGGCGGACCA  >  NC_002947/6005282‑6005534                                                                                                                                   |                                                                                                                           cgTCGGTGGAGTACGAGGCCAGCGGCTTCGGGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTgcc                                                                                                                       <  1:234024/135‑1 (MQ=255) aGTCGGTGGAGTACGAGGCCAGCGGCTTCGTGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTgcc                                                                                                                       >  1:169522/1‑136 (MQ=255)          aGTACGAGGCCAGCGGCTTCGTGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAg                                                                                                              <  2:1834/136‑1 (MQ=255)           gTACGAGGCCAGCGGCTTCGTGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAgc                                                                                                             <  2:402901/136‑1 (MQ=255)                  gCCAGCGGCTTCGTGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAgc                                                                                                      <  1:591780/136‑1 (MQ=255)                         gCTTCGTGGCCGATGCGGGCGCCTAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCAGTGCCGCCGCGCAGCCGCCAGCGCGAGGt                                                                                               >  1:321259/1‑136 (MQ=255)                         gCTTCGTGGCCGATGCGGGCGCCTAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCAGTGCCGCCGCGCAGCCGCCAGCGCGAGGt                                                                                               >  1:324477/1‑136 (MQ=255)                         gCTTCGTGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGt                                                                                               >  1:176/1‑136 (MQ=255)                                 gCCGATGCGGGCGCCAATGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCATCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATg                                                                                       <  1:529371/136‑1 (MQ=255)                                 gCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATg                                                                                       <  1:151907/136‑1 (MQ=255)                                                  gggCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCcgc                                                                      <  1:211217/136‑1 (MQ=255)                                                    gCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCcgctc                                                                    <  1:21165/136‑1 (MQ=255)                                                         gCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGa                                                               <  1:162951/136‑1 (MQ=255)                                                                      tGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAgcg                                                  <  2:176/136‑1 (MQ=255)                                                                      tGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCAGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAgcg                                                  <  2:324477/136‑1 (MQ=255)                                                                      tGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCAGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAgcg                                                  <  2:321259/136‑1 (MQ=255)                                                                            aCCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGc                                            <  2:226451/136‑1 (MQ=255)                                                                                          cACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGCGTTGAGCCAGCGCa                              >  2:597739/1‑136 (MQ=255)                                                                                                                tgtgGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGCGTTGAGCCAGCGCATGACTGCAGTCGAGTACCTGGc        <  2:137197/136‑1 (MQ=255)                                                                                                                      ccggcGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGCGTTGAGCCAGCGCATGACTGCAGTCGAGTACCTGGCGGACCa  <  1:39581/132‑1 (MQ=255)                                                                                                                      aCGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGCGTTGAGCCAGCGCATGACTGCAGTCGAGTACCTGGCGGACCa  <  1:161082/136‑1 (MQ=255)                                                                                                                      aCGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGGGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGCGTTGAGCCAGCGCATGACTGCAGTCGAGTACCTGGCGGACCa  <  1:541389/136‑1 (MQ=255)                                                                                                                                   |                                                                                                                           AGTCGGTGGAGTACGAGGCCAGCGGCTTCGTGGCCGATGCGGGCGCCAAGGGCGATGCACGTAACCAGCTGTTCGACCTGGCGCACCGGCACCTGGAGGCCAGCGGGGTGATGTGGAACGTCGACCTGGCCGTGCCGCCGCGCAGCCGCCAGCGCGAGGTGCTGGATGAGGTGCGAGTGTTCCGCTCGCTGAGTGACGAAGAGCGCGATGCGTTGAGCCAGCGCATGACTGCAGTCGAGTACCTGGCGGACCA  >  NC_002947/6005282‑6005534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑