Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,042,454	G→C	100%	A131P (GCC→CCC)	PP_5295 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,042,454	0	G	C	100.0%	148.0 / NA	44	A131P (GCC→CCC)	PP_5295	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (24/20);  total (24/20)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑GCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGGCGCCAGCGCGAA  >  NC_002947/6042313‑6042594                                                                                                                                               |                                                                                                                                              cAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCt                                                                                                                                        >  8:298124/1‑149 (MQ=255)        tGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCg                                                                                                                                 >  2:42699/1‑149 (MQ=255)            aCGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTg                                                                                                                              <  7:298124/148‑1 (MQ=255)             cgcgGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGa                                                                                                                            <  5:302138/149‑1 (MQ=255)              gcgGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGaa                                                                                                                           <  2:319475/149‑1 (MQ=255)                cgCAAGGTGGATTCCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGCCCATGAGCGCGCCGCGGAACTGGTTCACGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCGGAACCCGCTGGAAGt                                                                                                                         >  6:210651/2‑149 (MQ=255)                cgCAAGGTGCATTCCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGt                                                                                                                         >  8:301343/2‑149 (MQ=255)                 gCAAGGTGCATTCCACCTTGCGCCTCAACGAAGACAGAAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTCGTTGCCGCCGCCGAGGAAGTGGCTGCACTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTg                                                                                                                        >  1:132938/1‑149 (MQ=255)                  cAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGc                                                                                                                       >  1:249567/1‑149 (MQ=255)                       tGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCa                                                                                                                  >  1:277303/1‑149 (MQ=255)                               aCCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAgg                                                                                                          >  8:101390/1‑149 (MQ=255)                               aCCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAgg                                                                                                          >  6:287171/1‑149 (MQ=255)                               aCCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAgg                                                                                                          >  7:28181/1‑149 (MQ=255)                                ccTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGc                                                                                                         <  6:122939/149‑1 (MQ=255)                                     cgcCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGc                                                                                                     <  3:5741/148‑1 (MQ=255)                                     cgcCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCt                                                                                                    <  6:293053/149‑1 (MQ=255)                                      gcCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCtg                                                                                                   <  1:273138/149‑1 (MQ=255)                                         tCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCtggtg                                                                                                >  3:76515/1‑149 (MQ=255)                                                    aGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGCCCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCtg                                                                                                   <  5:14685/136‑1 (MQ=255)                                                    aGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGCCCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCtg                                                                                                   >  6:14685/1‑136 (MQ=255)                                                             aaaCCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGAcccc                                                                              >  4:188280/1‑147 (MQ=255)                                                               aCCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAgg                                                                          <  3:188280/149‑1 (MQ=255)                                                                     gAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGt                                                                                                                    <  5:163131/101‑1 (MQ=255)                                                                     gAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGt                                                                                                                    >  6:163131/1‑101 (MQ=255)                                                                       aaGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AAc                                                                 <  2:249567/149‑1 (MQ=255)                                                                                aTGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGcc                                                        <  7:101390/149‑1 (MQ=255)                                                                                 tGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACACGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGCCCCCCAGGCACTG‑AACGCCGAAGCCa                                                       >  2:87260/1‑149 (MQ=255)                                                                                                       gccgccgccGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTATCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCt                                 >  6:99397/1‑149 (MQ=255)                                                                                                         cgccgccgAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTcgac                                         >  7:296419/1‑139 (MQ=255)                                                                                                         cgccgccgAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTcgac                                         <  8:296419/139‑1 (MQ=255)                                                                                                               cgAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAAc                         <  2:277303/149‑1 (MQ=255)                                                                                                                     aGTGGCTGGCCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAgg                   <  7:301343/149‑1 (MQ=255)                                                                                                                       tggctggcCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGcc                 >  6:286872/1‑149 (MQ=255)                                                                                                                           tggcCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAgg                                                                          >  1:161912/1‑89 (MQ=255)                                                                                                                           tggcCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAgg                                                                          <  2:161912/89‑1 (MQ=255)                                                                                                                             gcCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTGCTACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGgc             <  1:87260/148‑1 (MQ=255)                                                                                                                             gcCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGgc             <  4:76515/147‑1 (MQ=255)                                                                                                                             gcCTGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGgc             <  8:28181/147‑1 (MQ=255)                                                                                                                                tGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGGCGCCAg        >  5:253463/1‑149 (MQ=255)                                                                                                                                tGATCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGGCGCCAg        >  8:195824/1‑149 (MQ=255)                                                                                                                                  aTCAAGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGCAGCGCTGGCCGAACTCAAGGCCGGGCGCCAgcg      >  5:246746/1‑149 (MQ=255)                                                                                                                                     aaGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGCCACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGTGCTGGCCGAACTCAAGGCCGGGCGCCAGcgcg    >  1:17454/1‑148 (MQ=255)                                                                                                                                     aaGCAGCCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGGCGCCAGcgcg    >  4:253570/1‑148 (MQ=255)                                                                                                                                      agcagcCG‑CCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGGCGCCAGCGCGaa  <  6:246746/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              CAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCG‑GCCCCCCTGAACCCGCTGGAAGTGCTGTCATGGCCAGGCGTGCTGGTGGCCGACACCAGCGACCCCCAGGCACTG‑AACGCCGAAGCCATGGCGCTGTTCGACGAAGCGCTGGCCGAACTCAAGGCCGGGCGCCAGCGCGAA  >  NC_002947/6042313‑6042594

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