Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	2,601,596	A→G	T324A (ACC→GCC)	PP_2277 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,601,596	0	A	G	84.6%	28.6 / ‑0.7	13	T324A (ACC→GCC)	PP_2277	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (6/5);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CACTTGATTTGTATTTCTCCGTAGCTACTGGCGCACACCGTTACACAGATGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAACCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGATGGGTCCACGGCAGGACTCGCTATGGGAACGAATTGACAC  >  NC_002947/2601501‑2601725                                                                                                |                                                                                                                                  cACTTGATTTGTATTTCTCCGTAGCTACTGGCGCACACCGTTACACAGATGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCg                                                                                           >  1:583125/1‑136 (MQ=255)            tATTTCTCCGTAGCTACTGGCGCACACCGTTACACAGATGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTgg                                                                                <  1:534985/136‑1 (MQ=255)            tATTTCTCCGTAGCTACTGGCGCACACCGTTACACAGATGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTgg                                                                                <  1:814875/136‑1 (MQ=255)                                                gATGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTa                                            <  1:720591/136‑1 (MQ=255)                                                gATGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAACCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTa                                            >  1:121940/1‑136 (MQ=255)                                                  tGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGa                                          >  1:653684/1‑136 (MQ=255)                                                    gTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTgcgc                                                           >  1:151568/1‑117 (MQ=255)                                                    gTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTgcgc                                                           <  2:151568/117‑1 (MQ=255)                                                                        cAAGGTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGATGGGTCCACGGCAGGACTCGCt                    >  2:886854/1‑136 (MQ=255)                                                                            gTTTGGAGACCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGATGGGTCCACGGCAGGACTCGCTATgg                <  1:886854/136‑1 (MQ=255)                                                                                  agaCCTTCCATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGATGGGTCCACGGCAGGACTCGCTATGGGAACGa          >  1:235786/1‑136 (MQ=255)                                                                                         ccATAGAGCCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGATGGGTCCACGGCAGGACTCGCTATGGGAACGAATTGaca   >  2:17575/1‑136 (MQ=255)                                                                                          cATAGAACCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGATGGGTCCACGGCAGGACTCGCTATGGGAACGAATTGacac  >  1:312842/1‑136 (MQ=255)                                                                                                |                                                                                                                                  CACTTGATTTGTATTTCTCCGTAGCTACTGGCGCACACCGTTACACAGATGGTCGATTCCTTGCTTTGATTCAAGGTTTGGAGACCTTCCATAGAACCATTAGTGACGAAACCCTTAAAGATGAAGCTGTTTTTCGAAGAGATTTGGCTAGGGTGCTATGGAGTGCGCCAAAAGACTTGCGTAGATGGGTCCACGGCAGGACTCGCTATGGGAACGAATTGACAC  >  NC_002947/2601501‑2601725

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