Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,408,553	G→A	A331A (GCG→GCA)	alaA →	glutamate‑pyruvate aminotransferase, glutamic‑pyruvic transaminase (GPT), alanine transaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,408,553	0	G	A	100.0%	70.3 / NA	21	A331A (GCG→GCA)	alaA	glutamate‑pyruvate aminotransferase, glutamic‑pyruvic transaminase (GPT), alanine transaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/9);  total (12/9)

GCTAGGTGGTTATCAGAGCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCGCTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACTGGCCGTGGCCGGA  >  NC_000913/2408427‑2408684                                                                                                                               |                                                                                                                                    gcTAGGTGGTTATCAGAGCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTc                                                                                                                         >  2:536777/1‑139 (MQ=255)                 agCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCa                                                                                                         >  1:114145/1‑139 (MQ=255)                 agCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCa                                                                                                         <  2:216654/139‑1 (MQ=255)                   cATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATg                                                                                                                            <  1:151807/118‑1 (MQ=255)                   cATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATg                                                                                                                            >  2:151807/1‑118 (MQ=255)                   cATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCaaa                                                                                                       <  1:536777/139‑1 (MQ=255)                                        tGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCa                                                                                                         <  2:106032/116‑1 (MQ=255)                                        tGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCa                                                                                                         >  1:106032/1‑116 (MQ=255)                                                        gAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTgg                                                                  >  2:475275/1‑139 (MQ=255)                                                         agcagcGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCaa                                                                                                        >  1:401548/1‑100 (MQ=255)                                                         agcagcGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCaa                                                                                                        <  2:401548/100‑1 (MQ=255)                                                          gcagcGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGAt                                                                >  2:210033/1‑139 (MQ=255)                                                            agcGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGAttt                                                              >  1:534162/1‑139 (MQ=255)                                                                        cgTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGaaaa                                                                           <  1:478082/114‑1 (MQ=255)                                                                        cgTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGaaaa                                                                           >  2:478082/1‑114 (MQ=255)                                                                                             ttCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAgg                             >  2:61846/1‑139 (MQ=255)                                                                                                    cGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCa                      <  1:210033/139‑1 (MQ=255)                                                                                                               tGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGa                                                                 <  1:313887/85‑1 (MQ=255)                                                                                                               tGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGa                                                                 >  2:313887/1‑85 (MQ=255)                                                                                                                       cGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACTGGCCGTGGCCgg   >  1:445298/1‑139 (MQ=255)                                                                                                                        gtggtgCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACTGGCCGTGGCCGGa  <  2:114145/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                    GCTAGGTGGTTATCAGAGCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCGCTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACTGGCCGTGGCCGGA  >  NC_000913/2408427‑2408684

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