Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,237,369	G→A	12.9%	G274S (GGC→AGC)	PP_2835 →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,237,369	0	G	A	12.9%	72.3 / 3.1	31	G274S (GGC→AGC)	PP_2835	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (14/13);  new base A (2/2);  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 = 9.62e-01

CCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGACAGGCCTTGGAGCAGGGC  >  NC_002947/3237236‑3237506                                                                                                                                      |                                                                                                                                          ccACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCt                                                                                                                           <  2:216555/150‑1 (MQ=255)  cacaAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGc                                                                                                                                 <  2:521932/143‑1 (MQ=255)  cacaAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGc                                                                                                                                 >  1:521932/1‑143 (MQ=255)    caAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTcga                                                                                                                        <  2:510017/150‑1 (MQ=255)    caAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGACTCAGTGGGTCACTTCTGTCCCTGCTCTTTTTAGTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTcga                                                                                                                        <  1:203872/150‑1 (MQ=255)                       cTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTgggg                                                                                                     <  2:33443/150‑1 (MQ=255)                           tGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGAt                                                                                                 <  2:125057/150‑1 (MQ=255)                                gAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGtt                                                                                            >  2:716975/1‑150 (MQ=255)                                   gCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTg                                                                                         <  1:50441/150‑1 (MQ=255)                                     ccTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGc                                                                                       >  2:537124/1‑150 (MQ=255)                                       tCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCaa                                                                                     >  2:624687/1‑150 (MQ=255)                                           cTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCAGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGAt                                                                                                 >  2:397900/1‑134 (MQ=255)                                           cTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCAGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGAt                                                                                                 >  2:397858/1‑134 (MQ=255)                                           cTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCAGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGAt                                                                                                 <  1:397900/134‑1 (MQ=255)                                           cTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCAGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGAt                                                                                                 <  1:397858/134‑1 (MQ=255)                                                   tGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAg                                                                         >  1:453912/1‑150 (MQ=255)                                                                                  tCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCt                                                                                          >  2:15759/1‑102 (MQ=255)                                                                                  tCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCt                                                                                          <  1:15759/102‑1 (MQ=255)                                                                                      tGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGAt                                                                                                 >  1:608597/1‑91 (MQ=255)                                                                                      tGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGAt                                                                                                 <  2:608597/91‑1 (MQ=255)                                                                                      tGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTgg                                                                                              <  2:407066/94‑1 (MQ=255)                                                                                      tGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTgg                                                                                              >  1:407066/1‑94 (MQ=255)                                                                                      tGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCa                                      >  1:475712/1‑150 (MQ=255)                                                                                       ggggCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAg                                     <  1:624687/150‑1 (MQ=255)                                                                                          gcgcCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGcc                                  >  2:538386/1‑150 (MQ=255)                                                                                               tCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTAGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGcc                             >  1:117041/1‑150 (MQ=255)                                                                                                         cgcTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGAc                   >  2:577252/1‑150 (MQ=255)                                                                                                            tGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCGGCCCACTAGACCCCGAACGCCCAATACGACTGCCCGCGCAATCAGGCCTGTCCCTGCGTCGACatg                >  2:351158/1‑148 (MQ=255)                                                                                                            tGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGACAgg                >  2:229967/1‑150 (MQ=255)                                                                                                                   tGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGACAGGCCTTGGa         <  1:537124/150‑1 (MQ=255)                                                                                                                          cgcATGTTCGGTGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCGAGCAGGCCTTGCCCTGCGTCGACAGGCCTTGGAGCAGGGc  <  1:716975/150‑1 (MQ=255)                                                                                                                                      |                                                                                                                                          CCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGACAGGCCTTGGAGCAGGGC  >  NC_002947/3237236‑3237506

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