Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,039,703	A→C	19.0%	P33P (CCT→CCG) ‡	PP_1812 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,039,703	0	A	C	19.0%	61.3 / 7.8	31	P33P (CCT→CCG) ‡	PP_1812	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base C (3/3);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGG  >  NC_002947/2039573‑2039845                                                                                                                                   |                                                                                                                                                ggAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAgg                                                                                                                               <  2:23850/149‑1 (MQ=255)               tGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCgg                                                                                                                 >  5:144501/1‑149 (MQ=255)                                                 ggcgggcgTTGTAAGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAg                                                                               >  7:221771/1‑149 (MQ=255)                                                            tAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGa                                                                           <  1:10280/142‑1 (MQ=255)                                                            tAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGa                                                                           >  2:10280/1‑142 (MQ=255)                                                              ggAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGc                                                                 >  5:302236/1‑149 (MQ=255)                                                                       ccGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGAtt                                                        >  2:222869/1‑149 (MQ=255)                                                                         gggCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCg                                                      <  4:159882/149‑1 (MQ=255)                                                                           gCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTgggtgggt                                                                                   >  6:207158/1‑119 (MQ=255)                                                                           gCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTgggtgggt                                                                                   <  5:207158/119‑1 (MQ=255)                                                                                  gttccAGTGGGTATCGGACAAACGGTGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtgg                                                                                     <  8:491/106‑1 (MQ=255)                                                                                      cAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGcc                                         >  3:33075/1‑149 (MQ=255)                                                                                        gTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGACGCGTTTTGCTGACCGATTCGGTAACCATTCGCCg                                        <  1:222869/149‑1 (MQ=255)                                                                                         tGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGGTTTGCTGACCGATTCGGTAACCATTCGCCgg                                       <  6:144501/148‑1 (MQ=255)                                                                                               tCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTCGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTgggtgggt                                                                                   <  3:203653/99‑1 (MQ=255)                                                                                               tCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTCGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTgggtgggt                                                                                   >  4:203653/1‑99 (MQ=255)                                                                                                 ggACAAACGGGGCTTGGGCTACAGGCACAAGGCCCGGGCCTTGTGCGGTCCGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGGGGGTGGGGGAAGTTGA‑GCGTTTTGCTGACCGCTTCGGTAACCATTCGCCGGGGTAGTTg                               <  6:340897/148‑1 (MQ=255)                                                                                                  gACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtgg                                                                                     >  7:222918/1‑94 (MQ=255)                                                                                                  gACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtgg                                                                                     <  8:222918/94‑1 (MQ=255)                                                                                                          gggATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGc                     >  2:46422/1‑149 (MQ=255)                                                                                                           gtttttGGGTTCCGGGACATGGCCTTGGCTTTGTGGGTTTTGGTCCGGGCCGTGCGGGGGGGAATTGGCCATGGTTTGGGGGTGGGTGAAGTTGT‑TCGTTTTGCTGACTGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCg                    <  8:221771/146‑1 (MQ=255)                                                                                                                   ctgaacgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATtg                        <  1:113529/131‑1 (MQ=255)                                                                                                                   ctgaacgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATtg                        >  2:113529/7‑137 (MQ=255)                                                                                                                      cAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCAGGCCGAGCGGGGCGGAAACGGCCACGGTTtgt                                                                                           <  7:87642/68‑1 (MQ=255)                                                                                                                      cAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCAGGCCGAGCGGGGCGGAAACGGCCACGGTTtgt                                                                                           >  8:87642/1‑68 (MQ=255)                                                                                                                         gCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGc                     >  8:273189/1‑134 (MQ=255)                                                                                                                           acaAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCggg                                      <  6:140239/115‑1 (MQ=255)                                                                                                                           acaAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCggg                                      >  5:140239/1‑115 (MQ=255)                                                                                                                           acaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAg               <  2:229750/138‑1 (MQ=255)                                                                                                                           acaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAg               >  1:229750/1‑138 (MQ=255)                                                                                                                             aaGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAgg  >  5:185573/1‑149 (MQ=255)                                                                                                                                   |                                                                                                                                                GGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGA‑GCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGG  >  NC_002947/2039573‑2039845

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