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_0009132,408,5530GA100.0% 30.1 / NA 10A331A (GCG→GCAalaAglutamate‑pyruvate aminotransferase, glutamic‑pyruvic transaminase (GPT), alanine transaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/5);  total (5/5)

GACCGCGCTAGGTGGTTATCAGAGCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCGCTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACTGGCCGTGGCCGGA  >  NC_000913/2408421‑2408684
                                                                                                                                    |                                                                                                                                   
gACCGCGCTAGGTGGTTATCAGAGCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGTGTGAAACCTCGTGGTGCACTGtat                                                                                                                               >  1:399258/1‑139 (MQ=255)
                        cATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTGTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCaaa                                                                                                       <  2:21870/139‑1 (MQ=255)
                                      ttaCCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCa                                                                                         <  2:399258/139‑1 (MQ=255)
                                              ggCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATc                                                                                 >  2:547973/1‑139 (MQ=255)
                                                                            gcgTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCa                                                                                                         <  1:552091/85‑1 (MQ=255)
                                                                            gcgTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCa                                                                                                         >  2:552091/1‑85 (MQ=255)
                                                                                                                      aaaCCTCGTGGTGCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACtggccgt         >  1:680557/1‑139 (MQ=255)
                                                                                                                             gtggtgCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCAcgacga                                                                                   <  1:659231/58‑1 (MQ=255)
                                                                                                                             gtggtgCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCAcgacga                                                                                   >  2:659231/1‑58 (MQ=255)
                                                                                                                             gtggtgCACTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACTGGCCGTGGCCGGa  <  2:680557/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GACCGCGCTAGGTGGTTATCAGAGCATCAGTGAATTTATTACCCCTGGCGGTCGTCTTTATGAGCAGCGTAACCGCGCGTGGGAACTGATCAACGATATTCCGGGCGTTTCCTGCGTGAAACCTCGTGGTGCGCTGTATATGTTCCCGAAAATCGACGCCAAACGCTTTAACATTCACGACGATCAGAAAATGGTGTTGGATTTCCTGTTGCAGGAAAAAGTTCTGTTGGTGCAAGGGACGGCATTCAACTGGCCGTGGCCGGA  >  NC_000913/2408421‑2408684

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