Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,098,121	A→C	E79A (GAA→GCA)	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,121	0	A	C	100.0%	87.0 / NA	25	E79A (GAA→GCA)	ghrA	glyoxylate/hydroxypyruvate reductase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/7);  total (18/7)

CATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGTTTTGACGATTATCGCATCCAGCAAAATAGTT  >  NC_000913/1097974‑1098255                                                                                                                                                    |                                                                                                                                       cATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGcaa                                                                                                                                      <  1:65002/150‑1 (MQ=255)   tGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGc                                                                                                                                   <  2:8031/151‑1 (MQ=255)          gCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAAcc                                                                                                                             >  2:60818/1‑150 (MQ=255)          gCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAAcc                                                                                                                             >  2:10620/1‑150 (MQ=255)                      aCTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCAc                                                                                                                 <  1:319330/150‑1 (MQ=255)                                    gCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAg                                                                                                  >  2:46/1‑151 (MQ=255)                                       ttAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATa                                                                                               >  1:201548/1‑151 (MQ=255)                                        tAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATAc                                                                                              <  2:6101/151‑1 (MQ=255)                                         aGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATAc                                                                                              >  1:272880/1‑150 (MQ=255)                                             tGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTa                                                                                         <  1:51056/151‑1 (MQ=255)                                                                tGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGa                                                                      >  1:106943/1‑151 (MQ=255)                                                                    ggCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAatat                                                                  >  2:257650/1‑151 (MQ=255)                                                                         ggCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTg                                                              <  1:300315/150‑1 (MQ=255)                                                                               aTCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTgtcagtca                                                       >  2:309519/1‑151 (MQ=255)                                                                               aTCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTgtcagtca                                                       >  2:56989/1‑151 (MQ=255)                                                                               aTCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTgtcagtca                                                       >  1:305316/1‑151 (MQ=255)                                                                                    aaaGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGtg                                                   >  1:46125/1‑150 (MQ=255)                                                                                            gtTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTgg                                          >  1:126982/1‑151 (MQ=255)                                                                                             ttCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTgg                                          <  1:218745/150‑1 (MQ=255)                                                                                                   cTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGAc                                   <  2:126981/151‑1 (MQ=255)                                                                                                       gggCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGttt                               >  1:67067/1‑151 (MQ=255)                                                                                                       gggCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGttt                               >  1:190838/1‑151 (MQ=255)                                                                                                             gtgtTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGTTTTGACGa                         >  2:133516/1‑151 (MQ=255)                                                                                                                                  aGCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGTTTTGACGATTATCGCATCCAGCAAAATAg    >  2:294424/1‑151 (MQ=255)                                                                                                                                   gCTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGCATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGTTTTGACGATTATCGCATCCAGCAAAATAGt   >  1:102905/1‑151 (MQ=255)                                                                                                                                    cTACAGGCACACCCTGCAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGTTTTGACGATTATCGCATCCAGCAAAATAGtt  >  2:32085/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                       CATGGAAAAGCGGAGATAATGACTCTGCTGATTATGCTTTAGTCTGGCATCCTCCTGTTGAAATGCTGGCAGGGCGCGATCTTAAAGCGGTGTTCGCACTCGGGGCCGGTGTTGATTCTATTTTGAGCAAGCTACAGGCACACCCTGAAATGCTGAACCCTTCTGTTCCACTTTTTCGCCTGGAAGATACCGGTATGGGCGAGCAAATGCAGGAATATGCTGTCAGTCAGGTGCTGCATTGGTTTCGACGTTTTGACGATTATCGCATCCAGCAAAATAGTT  >  NC_000913/1097974‑1098255

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