Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	261,052	A→G	D222G (GAC→GGC)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	261,052	0	A	G	95.0%	62.1 / ‑3.9	20	D222G (GAC→GGC)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (10/9);  total (11/9)
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.54e-01

CTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAA  >  NC_000913/260917‑261178                                                                                                                                        |                                                                                                                               cTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGt                                                                                                                             >  1:322956/1‑139 (MQ=255)    aCCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGc                                                                                                                          >  1:205348/1‑139 (MQ=255)                   cGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTAt                                                                                                           <  1:343797/139‑1 (MQ=255)                                  ttAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGc                                                                                            >  1:79565/1‑139 (MQ=255)                                      aGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACT‑GCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCgg                                                                                       >  1:17776/1‑139 (MQ=255)                                                ggCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACT‑GCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGcc                                                                             <  2:17776/139‑1 (MQ=255)                                                 gCATTGATGACGCACTGCGCGCGATTGCAGGTGACAGCGTTTCAGGCATCGGAACTGGCGGCGTGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGTAAGcc                                                                             <  2:104776/139‑1 (MQ=255)                                                  cATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCg                                                                            >  2:215931/1‑139 (MQ=255)                                                    ttGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCggg                                                                          >  1:569246/1‑139 (MQ=255)                                                    ttGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCggg                                                                          <  2:569246/139‑1 (MQ=255)                                                       atgaCGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGac                                                                                                       <  1:265029/107‑1 (MQ=255)                                                       atgaCGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGac                                                                                                       >  2:265029/1‑107 (MQ=255)                                                               cTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGgtg                                                               <  1:215931/139‑1 (MQ=255)                                                                          tGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGCACCGAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGaa                                                    >  1:496593/1‑139 (MQ=255)                                                                                              gCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGt                                                                       >  1:127984/1‑100 (MQ=255)                                                                                              gCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGt                                                                       <  2:127984/100‑1 (MQ=255)                                                                                                                  aGTACCGAATTGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTc            <  2:496593/139‑1 (MQ=255)                                                                                                                            tGCAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGaaaa  >  1:91862/1‑139 (MQ=255)                                                                                                                              cAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGgt                                                                >  1:239190/1‑75 (MQ=255)                                                                                                                              cAGGCCGCTGGCGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGgt                                                                <  2:239190/75‑1 (MQ=255)                                                                                                                                        |                                                                                                                               CTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAA  >  NC_000913/260917‑261178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑