Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,573,892	C→A	13.3%	P486Q (CCG→CAG)	PP_3150 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,573,892	0	C	A	13.3%	89.5 / 2.4	30	P486Q (CCG→CAG)	PP_3150	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base A (2/2);  total (16/14)
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.60e-01

CCGGGCCGGCAATGGACAATGTGCTCAACCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACTTTGAACTGTTCGTCGAGGACT  >  NC_002947/3573766‑3574026                                                                                                                               |                                                                                                                                       ccgggccggCAATGGACAATGTGCTCAACCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGcc                                                                                                                               <  2:341787/136‑1 (MQ=255)        ggCAATGGACAATGTGCTCAACCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGgcg                                                                                                                        >  1:332962/1‑136 (MQ=255)                       gCTCAACCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCgatgat                                                                                                         <  1:136977/136‑1 (MQ=255)                           aaCCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCAGGAACAGCCCACGGCGCTTATCGTCGATGATGAAt                                                                                                     >  1:260607/1‑136 (MQ=255)                           aaCCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCAGGAACAGCCCACGGCGCTTATCGTCGATGATGAAt                                                                                                     >  1:403804/1‑136 (MQ=255)                           aaCCCGTTGATCCAGGCGCACCCGGCAACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAAt                                                                                                     <  1:53836/136‑1 (MQ=255)                            aCCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATg                                                                                                    >  2:387337/1‑136 (MQ=255)                            aCCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATg                                                                                                    >  2:360550/1‑136 (MQ=255)                                   gATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTg                                                                                             >  2:418568/1‑136 (MQ=255)                                             cACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATAAGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGt                                                                                   >  1:191487/1‑136 (MQ=255)                                                    ccACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGATCAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAgcg                                                                            <  1:213230/136‑1 (MQ=255)                                                    ccACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAgcg                                                                            <  2:39270/136‑1 (MQ=255)                                                      aCGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGcc                                                                          <  1:88016/136‑1 (MQ=255)                                                                     ttCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCAAGCGCCTCACGCATCAggggg                                                           >  1:148393/1‑136 (MQ=255)                                                                       cAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGCATGTTACCACCGCAACACCAAACAGGGGGTg                                                         >  1:419467/1‑136 (MQ=255)                                                                             cAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTTATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTc                                                   <  1:74013/136‑1 (MQ=255)                                                                                          tACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGaaaa                                      <  1:141831/136‑1 (MQ=255)                                                                                               ccGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCCCGCATCAGGGTGTGCCGTTCTTGTTGCTGAAAAGCGCt                                 >  2:493144/1‑136 (MQ=255)                                                                                               ccGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTAACGTCGATGATGATTGGGGCCTGTTTGGCTGGTTACAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAgcgca                                 >  1:458841/1‑135 (MQ=255)                                                                                                 gggTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCg                               <  1:70260/136‑1 (MQ=255)                                                                                                  ggTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGc                              <  1:387337/136‑1 (MQ=255)                                                                                                    tCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCACGCaa                            >  1:74680/1‑136 (MQ=255)                                                                                                     cAGCTTCTTCGATACGGAGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGAAGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAg                           >  2:300836/1‑136 (MQ=255)                                                                                                      aGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAgg                          <  1:158790/136‑1 (MQ=255)                                                                                                       gCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGa                         >  2:95527/1‑136 (MQ=255)                                                                                                       gCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTAGCAAggg                         >  2:404075/1‑135 (MQ=255)                                                                                                           gatctaTACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGTTGAATGGGGCCTGTTTGGCTGGTTCCAGCGGCTCACGCATCAGGGGGTGGCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACttt                     <  1:236918/131‑1 (MQ=255)                                                                                                                    ggTGATGGTGCAGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACTTTGAACTGTtc            <  2:260607/136‑1 (MQ=255)                                                                                                                    ggTGATGGTGCAGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACTTTGAACTGTtc            <  2:403804/136‑1 (MQ=255)                                                                                                                         tGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACTTTGAACTGTTCGTCGa       >  1:244528/1‑136 (MQ=255)                                                                                                                              ccGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACTTTGAACTGTTCGTCGAGGACt  <  1:364142/136‑1 (MQ=255)                                                                                                                              ccGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTACCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACTTTGAACTGTTCGTCGAGGACt  <  1:418568/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                       CCGGGCCGGCAATGGACAATGTGCTCAACCCGTTGATCCAGGCGCACCCGGCCACGATCCAGGCGGCCTTCAGTGGCAGCGCCGAGCACTACACCCGGGTCAGCTTCTTCGATACGGTGATGGTGCCGGAACAGCCCACGGCGCTTATCGTCGATGATGAATGGGGCCTGTTTGGCTGGTTCCAGCGCCTCACGCATCAGGGGGTGCCGTTCTTGTTGCTGAAAAGCGCTCGCAAGGACTTTGAACTGTTCGTCGAGGACT  >  NC_002947/3573766‑3574026

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