Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	864,723	C→A	15.3%	T155T (ACG→ACT)	uraA ←	uracil permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	864,723	0	C	A	15.3%	79.3 / 2.0	26	T155T (ACG→ACT)	uraA	uracil permease
Reads supporting (aligned to +/- strand):  ref base C (7/15);  new base A (2/2);  total (9/17)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

GCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAGATGATCACCGGGCCGATCACCACCGGC  >  NC_002947/864594‑864837                                                                                                                                  |                                                                                                                   gCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTg                                                                                                              <  2:447279/136‑1 (MQ=255) gCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTg                                                                                                              <  2:371649/136‑1 (MQ=255)      aGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACgg                                                                                                         <  2:459547/136‑1 (MQ=255)      aGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACgg                                                                                                         <  1:475513/136‑1 (MQ=255)      aGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACgg                                                                                                         <  2:199993/136‑1 (MQ=255)        gacagcGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCa                                                                                                       <  1:6016/131‑1 (MQ=255)           aaCGACAGGGCAAAGCCCACCAGCACGCCGGAGATGAAGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCCGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCa                                                                                                    <  1:279386/136‑1 (MQ=255)              gACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTg                                                                                                 >  1:408033/1‑136 (MQ=255)              gACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTg                                                                                                 >  2:356351/1‑136 (MQ=255)              gACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCAAAAGTg                                                                                                 >  1:393636/1‑136 (MQ=255)              gACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTAACAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTg                                                                                                 >  1:540306/1‑136 (MQ=255)              gACAGGGCAAAGCCAACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCAGTCTTGTACGGCATCAGTg                                                                                                 >  1:141914/1‑136 (MQ=255)              gACAGGGCAAAGCCAACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCAGTCTTGTACGGCATCAGTg                                                                                                 >  1:275303/1‑136 (MQ=255)                       aaGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGt                                                                                        <  1:144491/136‑1 (MQ=255)                                       cGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCAAt                                                                        <  2:99788/136‑1 (MQ=255)                                          agaTGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGc                                                                     >  2:247774/1‑136 (MQ=255)                                           gatgatGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGcc                                                                    <  1:152048/136‑1 (MQ=255)                                                    tACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTggcg                                                           >  1:512226/1‑136 (MQ=255)                                                            ggAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCAACGCCAATTTGGCGGCGATCgg                                                   >  1:496777/1‑136 (MQ=255)                                                                  ttcccTTGCCGATTACGGCAACGATGAAGGTTGTCAGCAGCGCTGTCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCAGTACCGTCCCCACCTTTGCCCATCGGCATGTTGGCGGCGATCGGTGcctt                                             <  2:393636/134‑3 (MQ=255)                                                                  tGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCAt                                             <  2:408033/136‑1 (MQ=255)                                                                  tGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCAGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCAt                                             <  2:141914/136‑1 (MQ=255)                                                                  tGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCAGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCAt                                             <  2:275303/136‑1 (MQ=255)                                                                            gAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCg                                   <  1:515294/136‑1 (MQ=255)                                                                                        gATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAGATGATc                       <  1:58794/136‑1 (MQ=255)                                                                                                             tGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAGATGATCACCGGGCCGATCACCACcggc  <  2:512226/136‑1 (MQ=255)                                                                                                                                  |                                                                                                                   GCCAAAGGCGAACGACAGGGCAAAGCCCACCAGCACGCCGGAGATGATGGGTACCAGGCGGAAAATGCCCTTGCCGAATACGGCAACGATCAAGGTTGTCAGCAGCGCTGGCATCGAGATCAGCATGGCCGTCTTGTACGGCATCAGTGCCGTACCGTCCCCACCTTTGCCCATCGCCATGTTGGCGGCGATCGGTGCCATGGCCAGGCCGATGGAGATGATCACCGGGCCGATCACCACCGGC  >  NC_002947/864594‑864837

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