Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,851,740	G→T	63.6%	G262C (GGC→TGC)	relA →	ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,851,740	0	G	T	63.6%	12.2 / 10.4	19	G262C (GGC→TGC)	relA	ATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (6/6);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATC  >  NC_002947/1851596‑1851874                                                                                                                                                 |                                                                                                                                       ttCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGGCGACATCAGCTGCCgg                                                                                                                                   >  2:201793/1‑150 (MQ=255)           tGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACata                                                                                                                         >  1:190601/1‑150 (MQ=255)             gAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATc                                                                                                                       <  1:290438/150‑1 (MQ=255)                                     gCCAAGCTGACGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCAACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAg                                                                                               >  2:20692/1‑150 (MQ=255)                                                  aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTAt                                                                                                                    <  2:157526/116‑1 (MQ=255)                                                  aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTAt                                                                                                                    >  1:157526/1‑116 (MQ=255)                                                  aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTgg                                                                                  >  1:202177/1‑150 (MQ=255)                                                                          gcTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtgc                                                          <  1:20692/150‑1 (MQ=255)                                                                            ttCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtgcgt                                                        >  2:2219/1‑150 (MQ=255)                                                                            ttCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtgcgt                                                        >  1:291680/1‑150 (MQ=255)                                                                                           atgaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTg                                         <  2:202177/150‑1 (MQ=255)                                                                                             gaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGt                                       <  2:190601/150‑1 (MQ=255)                                                                                                ccAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGcc                                    >  2:383619/1‑150 (MQ=255)                                                                                                                  gCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTa                  <  1:383619/150‑1 (MQ=255)                                                                                                                      gCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGa                                                                                 >  1:48070/1‑83 (MQ=255)                                                                                                                      gCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGa                                                                                 <  2:48070/83‑1 (MQ=255)                                                                                                                             gCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTgg       <  2:230509/150‑1 (MQ=255)                                                                                                                             gCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTgg       <  2:291680/150‑1 (MQ=255)                                                                                                                                  aaGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATc  <  1:461464/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                       TTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATC  >  NC_002947/1851596‑1851874

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