Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,897,454	A→G	23.5%	T175A (ACC→GCC)	sdaA →	L‑serine dehydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,897,454	0	A	G	23.5%	25.8 / 6.1	17	T175A (ACC→GCC)	sdaA	L‑serine dehydratase 1
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base G (2/2);  total (8/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 = 8.95e-01

TAACGGCGATGAAGTCGTCTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGA  >  NC_000913/1897327‑1897601                                                                                                                                |                                                                                                                                                    taACGGCGATGAAGTCGTCTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTaa                                                                                                                               <  1:54624/150‑1 (MQ=255)     ggCGATGAAGTCGTCTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCGCCGAACTGCTCGCGTACTGTAATGaa                                                                                                                           >  1:112474/1‑150 (MQ=255)     ggCGATGAAGTCGTCTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGaa                                                                                                                           >  1:36518/1‑150 (MQ=255)     ggCGATGAAGTCGTCTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATAGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCGCCGAACTGCTCGCGTACTGTAATGaa                                                                                                                           >  1:112473/1‑150 (MQ=255)         aTGAAGTCGTCTACAGCAAAAATTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGGTCAAATCTGCCACCGAACTGCGCGCGTACTGTAATGAAACCg                                                                                                                       >  2:241179/1‑150 (MQ=255)               tcgtcTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGGTGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTAtt                                                                                                                 >  1:73037/1‑150 (MQ=255)                cgtcTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTc                                                                                                                <  2:199165/150‑1 (MQ=255)                                                                   agaaCACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCGCCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCaagaaaga                                                             <  2:112473/150‑1 (MQ=255)                                                                   agaaCACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCGCCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCaagaaaga                                                             <  2:112474/150‑1 (MQ=255)                                                                   agaaCACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCaagaaaga                                                             <  2:107423/150‑1 (MQ=255)                                                                   agaaCACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCaagaaaga                                                             <  2:36518/150‑1 (MQ=255)                                                                      acacTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGAt                                                          <  1:49891/150‑1 (MQ=255)                                                                        acTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATcg                                                        >  1:115957/1‑150 (MQ=255)                                                                                                aaGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGc                                >  2:232607/1‑150 (MQ=255)                                                                                                                        aaTCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTatcgatc        >  1:198879/1‑150 (MQ=255)                                                                                                                              ccACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGa  <  1:12825/150‑1 (MQ=255)                                                                                                                              ccACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGa  <  2:57912/150‑1 (MQ=255)                                                                                                                                |                                                                                                                                                    TAACGGCGATGAAGTCGTCTACAGCAAAACTTATTATTCCATCGGCGGCGGTTTTATCGTCGATGAAGAACACTTTGGTCAGGATGCTGCCAACGAAGTAAGCGTGCCGTATCCGTTCAAATCTGCCACCGAACTGCTCGCGTACTGTAATGAAACCGGCTATTCGCTGTCTGGTCTCGCTATGCAGAACGAACTGGCGCTGCACAGCAAGAAAGAGATCGACGAGTATTTCGCGCATGTCTGGCAAACCATGCAGGCATGTATCGATCGCGGGA  >  NC_000913/1897327‑1897601

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