Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867	C→T	S177S (AGC→AGT)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867	0	C	T	95.5%	54.5 / ‑4.1	23	S177S (AGC→AGT)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/1);  total (22/1)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01

AGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCT  >  NC_000913/623‑1111                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     aGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCAtt                                                                                                                                                                                                                                                >  1:320239/1‑251 (MQ=255)    aaTTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCg                                                                                                                                                                                                                                             >  1:44684/1‑251 (MQ=255)    aaTTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCg                                                                                                                                                                                                                                             >  1:295145/1‑251 (MQ=255)    aaTTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGATGCGCTGATTTGACGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTCCCGTTATCGATCCGGTCGAAAAACTGCTGGCCGTGGGGCATTACCTCGACTCTACCGTCGATATTGATGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCg                                                                                                                                                                                                                                             >  1:252743/1‑251 (MQ=255)                          ccTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGCAACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGCGTCCACCCGCCGCATTCCGGCAAGTCGCATCCCGGCTGGTCACATGGGGCTGATgg                                                                                                                                                                                                                       >  2:93164/1‑251 (MQ=255)                                                            gATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGGATTGCGGCAAGGCGCATTCCGGCTGATCACATGGTACTGATGGCAGGTTTCCCCGCCGGTAATGAAAAAGGCGac                                                                                                                                                                                       >  2:276755/1‑248 (MQ=255)                                                               aGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACtggg                                                                                                                                                                                  >  1:117551/1‑250 (MQ=255)                                                                                        gTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGac                                                                                                                                                         >  1:77444/1‑251 (MQ=255)                                                                                            cGAGAAAATGACGATCGCCATTATGGCAGGCGTATTAGAAGCGCGCGGTAACAACGTTACCGTTATCGATACGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGCTATTGCTGCGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCAAATGGTGCTGATGGCAGGTTTACCCTCCGGTAATGAAAAAGGCACACTGGTGGTGATTGGGCGGAAAGGTTCCGactacc                                                                                                                                                     >  2:44618/1‑250 (MQ=255)                                                                                                   aTGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTActccgcggc                                                                                                                                              >  1:96321/1‑245 (MQ=255)                                                                                                       cGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACACCGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGAAATGAAAAAGGCGAACTGGTGGTACTTGGGCGCAACGGTTCCGACTCCTCCGCGGCGGTg                                                                                                                                          >  2:150776/1‑251 (MQ=255)                                                                                                        gATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTACTGAAAAAGGCGAACTGGGGGAACTTGGACGAAACGGGTCCGaatactacgcggcggggg                                                                                                                                         >  2:208118/1‑239 (MQ=255)                                                                                                            gCCATTATGGCCGGCGTATTAGAAGCGCGCGGGCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTACGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGCCTACTCCGCGGCGGAgctggg                                                                                                                                     >  1:318659/1‑250 (MQ=255)                                                                                                             ccATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAACAACTGCTGGCAGAGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGATGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTgctggct                                                                                                                                    >  1:44685/1‑251 (MQ=255)                                                                                                                     gCCGGCGTATTAGAAGAGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTggtt                                                                                                                            >  1:329575/1‑248 (MQ=255)                                                                                                                             aTTAGAAGCGAGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAACATGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGTGTCCACCCGCCGTATTGAGGAAAGTCGCATTCCGGCTGATCCCATGGTGCTGATGGCAGGTTTAACCGCCGGTAATGAAAAAGGCGTACTGGTGGTCCTGGGACGCAACGGTTCAGACTACGCCGCGGCGGTGAGGGCTACCTGTTTACGCGCCg                                                                                                                    >  2:226888/1‑251 (MQ=255)                                                                                                                                                  cGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCAGACTACTCCGCGGCGGTGCTGGCTGCCTGATTGAGCGGCGATTTTTGCGAGATTTggacgg                                                                                               >  2:272852/1‑251 (MQ=255)                                                                                                                                                                                  tgGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGATTCCACCAGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGATTCACCGCCGGTAATGAAACAGGCGAACTGGTGGTACTTGGAAGAAACGGTTCCGACTAATCCGCGGCGGTGCTCGCTGGCTGGTTACTCGACGATTGTTACGGGATGTGGACGTACGTTGAGGGGGGCTTAACCTGCGAccccgtg                                                               >  2:90780/1‑247 (MQ=255)                                                                                                                                                                                              aTTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGAt                                                   >  1:105616/1‑251 (MQ=255)                                                                                                                                                                                                          cTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTTTTGCGGGATTTGGGCGGACGTTGACGGGGTCTATACCTGCGACCCACGGCCGGTGCCCGATGCGAGGttgttg                                       >  1:338988/1‑251 (MQ=255)                                                                                                                                                                                                                  gATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCCCCGATTGTTGCGCGATTTTGCCGCACGGTGCCGGGGTCTATACATGAGCCCCCCGTCAGGGGCCCGATTCGAAGTTGATGGAGTCggt                               >  1:247810/1‑249 (MQ=255)                                                                                                                                                                                                                                          aTTGCGGCAAGTCGCATCCCGGCTGATCACAGGGTGCTGATGGCAGTTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGc        <  1:93164/250‑1 (MQ=255)                                                                                                                                                                                                                                                gCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTACTGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGCTGGAGCTTTCCt  >  1:129529/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     AGGAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCT  >  NC_000913/623‑1111

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