Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,098,620	T→A	V245V (GTT→GTA)	ghrA →	glyoxylate/hydroxypyruvate reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,098,620	0	T	A	88.0%	67.3 / 2.9	25	V245V (GTT→GTA)	ghrA	glyoxylate/hydroxypyruvate reductase A
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (13/9);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

TGGAGTCGAACCCGTAAATCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGG  >  NC_000913/1098378‑1098857                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              tGGAGTCGAACCCGTACATCGTGGCCTGGCGTGCAACGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGgcgc                                                                                                                                                                                                                                       >  1:89591/1‑251 (MQ=255)                 aaTCGTGGCCTGGCGTGCAAAGCTTTGCCGGAAGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCTGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATTGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATATCGGCAAAGTAAAAGGCGCAATGTTTGATGttttt                                                                                                                                                                                                                       <  1:51635/251‑1 (MQ=255)                  aTCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGGATTGATTAATTTGTTACAGAATACCACTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGAGTTCATGTTGAGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGAAAAGTAAAGGGAGCAATGTTGGATGttttt                                                                                                                                                                                                                       >  1:188806/1‑250 (MQ=255)                        gCCTGGCGTGCAAAGCTTTGCCGGACTGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATcttg                                                                                                                                                                                                                <  2:89591/251‑4 (MQ=255)                                                                           ccAATGTCGGGTATTGAGTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATaa                                                                                                                                                             <  1:104723/251‑1 (MQ=255)                                                                                      tATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTc                                                                                                                                                  >  1:290174/1‑251 (MQ=255)                                                                                            ttAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCcgcg                                                                                                                                            >  1:14664/1‑251 (MQ=255)                                                                                              aaTTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGAt                                                                                                                                                                                                                             <  2:137087/168‑1 (MQ=255)                                                                                              aaTTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGAt                                                                                                                                                                                                                             >  1:137087/1‑168 (MQ=255)                                                                                                  tGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTATGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTAcc                                                                                                                                      >  1:168953/1‑251 (MQ=255)                                                                                                   gTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTAccc                                                                                                                                     <  1:102868/251‑1 (MQ=255)                                                                                                                      aCCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTggg                                                                                                                  >  1:107747/1‑250 (MQ=255)                                                                                                                        cGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGt                                                                                                                >  1:58583/1‑251 (MQ=255)                                                                                                                                            aaTTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGcc                                                                                            <  1:4627/251‑1 (MQ=255)                                                                                                                                                      aaaaTTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTgca                                                                                                                          >  2:78015/1‑209 (MQ=255)                                                                                                                                                      aaaaTTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGaa                                                                                                                          <  1:78015/211‑1 (MQ=255)                                                                                                                                                               ggATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAgggg                                                                            >  1:40057/1‑248 (MQ=255)                                                                                                                                                                          aTCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGCAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCggg                                                              >  2:225591/1‑251 (MQ=255)                                                                                                                                                                                       ggCGCGTGGGGTACATGTTGTGGAAGATGACCTGCTCGCGGCGGTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTGTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACcgcg                                                 <  2:132237/251‑1 (MQ=255)                                                                                                                                                                                            gtggtgTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCAACCGCGCAcgc                                            >  1:169953/1‑251 (MQ=255)                                                                                                                                                                                                   tCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTTCCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCAt                                                                                                <  1:179667/192‑1 (MQ=255)                                                                                                                                                                                                       tttgtGGAAGCTGACCAGCTCACGGTGCTTGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGGTTTTAATCGTGAACGCTTCCCGCCTTAAAGTCCGCTCTGGCAACATCCACGCATGACGATAACACCACATGTCGCCGCGATTGCCCTTCCCGCTGAAGCTGTGGAGTACATTGCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCTCGCGTCTACtaat                                  <  2:63526/249‑1 (MQ=255)                                                                                                                                                                                                               aGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCACGTCGACCGCGCACGCGGCTACTAATAAAGCATc                          >  1:55279/1‑250 (MQ=255)                                                                                                                                                                                                                            gcggcgCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTAt            <  1:39078/251‑1 (MQ=255)                                                                                                                                                                                                                                      aTAGCGGCAAAGTAAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGACAAAGGGGGGAGGGTCTGCGGGCAATTCGACCGCGCACGCGGCTACTAAGAAAGCATCAGTATTCCTGCTATCCTTGGCggg  >  2:69425/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              TGGAGTCGAACCCGTAAATCGTGGCCTGGCGTGCAAAGCTTTGCCGGACGGGAAGAACTGTCTGCATTTCTGAGCCAATGTCGGGTATTGATTAATTTGTTACCGAATACCCCTGAAACCGTCGGCATTATTAATCAACAATTACTCGAAAAATTACCGGATGGCGCGTATCTCCTCAACCTGGCGCGTGGTGTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAAGCTGTGGAGTACATTTCTCGCACCATTGCCCAGCTCGAAAAAGGGGAGAGGGTCTGCGGGCAAGTCGACCGCGCACGCGGCTACTAATAAAGCATCAGGATTCCTGCTATCCTTGGCGGG  >  NC_000913/1098378‑1098857

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