Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 378,772 A→G V267A (GTC→GCC)  frmA ← S‑(hydroxymethyl)glutathione dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913378,7720AG100.0% 55.6 / NA 17V267A (GTC→GCC) frmAS‑(hydroxymethyl)glutathione dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/8);  total (9/8)

TTCCATACGCGACCAGTGACCAACTGGAATGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGA  >  NC_000913/378643‑378901
                                                                                                                                 |                                                                                                                                 
ttCCATACGCGACCAGTGACCAACTGGGATGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCGCCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGt                                                                                                                                >  1:328581/1‑133 (MQ=255)
       cgcgACCAGTGACCAGCTCGAATGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGAt                                                                                                                         <  1:205879/133‑1 (MQ=255)
          gACCAGTGACCAACTGGAATGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCa                                                                                                                      >  1:102632/1‑133 (MQ=255)
           aCCAGTGACCAACTGGAATGGGCGGGTGGAGATTTCCTGACAGGCAACCGCGACCCCGAAGATCACCGACTGACCCCACCCGCGCTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCCATGCAt                                                                                                                     >  2:317901/1‑133 (MQ=255)
               gTGACCAACTCGAATTGACGGGTGTAGACTTCTTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCaa                                                                                                                 <  2:130419/133‑1 (MQ=255)
                      aCTGGAATGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATg                                                                                                          <  2:3330/133‑1 (MQ=255)
                          gAATGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCg                                                                                                      >  2:127831/1‑133 (MQ=255)
                            aTGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGAt                                                                                                    >  1:2297/1‑133 (MQ=255)
                              ggACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATAc                                                                                                  >  1:308172/1‑133 (MQ=255)
                                  gggTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCa                                                                                              <  1:127832/133‑1 (MQ=255)
                                                               aCCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTa                                                                 <  1:127831/133‑1 (MQ=255)
                                                                    gatgatCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGt                                                            >  1:370335/1‑133 (MQ=255)
                                                                    gacgatCGCCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGt                                                            <  2:328581/130‑1 (MQ=255)
                                                                                   aCCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTc                                             >  2:366555/1‑133 (MQ=255)
                                                                                    ccccAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCg                                            >  2:76349/1‑133 (MQ=255)
                                                                                              ggTGCGCACTTTCCAGCGCCGCACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCa                                  <  2:102632/133‑1 (MQ=255)
                                                                                                                ccgcACGCATCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCg                <  2:308172/133‑1 (MQ=255)
                                                                                                                   cacgcaTCACGTTGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGAcg             >  2:49398/1‑133 (MQ=255)
                                                                                                                              ttGGCGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGa  <  2:2297/133‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                 
TTCCATACGCGACCAGTGACCAACTGGAATGGACGGGTGGAGATTTCCTGACCGGCAACCGCGACCCCGATGATCACCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGA  >  NC_000913/378643‑378901

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