Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,434,815	G→T	21.0%	G253C (GGC→TGC)	PP_2129 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,434,815	0	G	T	21.0%	46.4 / 2.8	19	G253C (GGC→TGC)	PP_2129	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/10);  new base T (2/2);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GTGCGGTGCATCTTCACGACCATCGATCCAGTCACCGGGGAGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAACGACCGCTCAAACCCT  >  NC_002947/2434683‑2434925                                                                                                                                     |                                                                                                               gtgCGGTGCATCTTCAGGACCATCGATCCAGTCACCTGGGTGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAATGAAGGCGACGTGCTGTTCGGCCAGAAGCTGGCTGTGGATGGCa                                                                                                             <  2:404073/136‑1 (MQ=255)  tgCGGTGCATCTTCACGACCATCGATCCAGTCACCGGGGAGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAg                                                                                                            >  1:236536/1‑136 (MQ=255)                          aTCCAGTCACCGTGGAGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTTGCTGGAGGTGGGGATGGa                                                                                    >  2:14339/1‑136 (MQ=255)                                  cccGGGGAGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGaaa                                                                            <  2:500711/135‑1 (MQ=255)                                  cccGGGGAGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGaaa                                                                            <  1:460646/135‑1 (MQ=255)                                  aCCGGGGAGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGaaa                                                                            <  1:478518/136‑1 (MQ=255)                                                      ccGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCATAACCTGGCTGTGGATTGCAGTGGGTGGCTGGAGGTGGTGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGc                                                        >  2:36336/1‑136 (MQ=255)                                                      ccGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCATAACCTGGCTGTGGATTGCAGTGGGTGGCTGGAGGTGGTGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGc                                                        >  2:140435/1‑136 (MQ=255)                                                        gTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCt                                                      >  2:41285/1‑136 (MQ=255)                                                                  ggCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCggggc                                            >  1:296981/1‑133 (MQ=255)                                                                                    ccGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTTGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGc                          <  1:14339/136‑1 (MQ=255)                                                                                     cgcgAAAAAGAAGGCGACTTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGTGTGAACCCGCt                         <  1:141283/136‑1 (MQ=255)                                                                                     cgcgAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCt                         <  1:461943/136‑1 (MQ=255)                                                                                      gcgAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTc                        >  2:376020/1‑136 (MQ=255)                                                                                            aagaagGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACaa                  <  1:376020/136‑1 (MQ=255)                                                                                                   cGACGTGCTGTTCGGCCATAACCTGGCTGTGGATTGCAGTGGGTGGCTGGAGGTGGTGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAACGACCGc           <  1:36336/136‑1 (MQ=255)                                                                                                   cGACGTGCTGTTCGGCCATAACCTGGCTGTGGATTGCAGTGGGTGGCTGGAGGTGGTGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAACGACCGc           <  1:140435/136‑1 (MQ=255)                                                                                                      cGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAACGACCGCTCa        <  2:468572/136‑1 (MQ=255)                                                                                                            gTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAACGACCGCTCAAACCCt  <  1:198197/136‑1 (MQ=255)                                                                                                                                     |                                                                                                               GTGCGGTGCATCTTCACGACCATCGATCCAGTCACCGGGGAGCGTAGCGCTGACCGTGAGCCGATGGCTACGCTGAAGACTTTCCGCGAAAAAGAAGGCGACGTGCTGTTCGGCCAGAACCTGGCTGTGGATGGCAGTGGGTGGCTGGAGGTGGGGATGGAAGTAGAAATCCTGGAATAGGTGTGTTGCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAACGACCGCTCAAACCCT  >  NC_002947/2434683‑2434925

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