Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,048,383 A→G V220A (GTA→GCA)  gcvP ← glycine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,048,3830AG100.0% 76.0 / NA 23V220A (GTA→GCA) gcvPglycine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/13);  total (10/13)

TATCCGCGCCCTGTTTACCCGGCGCAGTTAACAGCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTACCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTGCGGATGCACATCGGAAGCCA  >  NC_000913/3048249‑3048518
                                                                                                                                      |                                                                                                                                       
tatCCGCGCCCTGTTTACCCGGCGCAGTTAACAGCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCctgcctgc                                                                                                                                    >  2:444585/1‑140 (MQ=255)
            gTTTACCCGGCGCAGTTAACAGCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACg                                                                                                                        <  2:120085/140‑1 (MQ=255)
                 cccGGCGCAGTTAACAGCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGaa                                                                                                                   >  2:78340/1‑140 (MQ=255)
                     gcgcAGTTAACAGCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACg                                                                                                               >  1:416766/1‑140 (MQ=255)
                         aGTTAACAGCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCt                                                                                                           <  1:631292/140‑1 (MQ=255)
                               cAGCACCAGCGCCATAATATCGGCGGCAACGCTGATCACTATTTTGCGCGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGt                                                                                                     <  2:246989/140‑1 (MQ=255)
                                 gCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCg                                                                                                   <  2:36544/140‑1 (MQ=255)
                                     cAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCa                                                                                               >  1:36285/1‑140 (MQ=255)
                                                 atCGGCGGCAACGCTGAACACAATTTTGCGTGATTTCAGTTCGCTAATTAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGt                                                                                   <  2:338481/140‑1 (MQ=255)
                                                     gcggcAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAGCAGCACGCCGAAGACGCCCTGATGGCCGAGTACTTTTTGCGCGTCATc                                                                               >  2:278519/1‑140 (MQ=255)
                                                           aCGCTGACCTCAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAAt                                                                         <  2:473103/140‑1 (MQ=255)
                                                                         tttGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCaa                                                           <  1:78340/140‑1 (MQ=255)
                                                                                gATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTGGCGCGTCATCGACAATCACTTGAAAACCAAAGGTTTc                                                    >  1:263174/1‑140 (MQ=255)
                                                                                  tttCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGCGCACTTTTCGCGAGTCTTCGACAATCACTTCAAAACCAAAGGTTTCgg                                                  >  2:281621/1‑140 (MQ=255)
                                                                                          cGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTa                                          <  2:36285/140‑1 (MQ=255)
                                                                                          cGCTAATAAGCGCAGCGTAGTCGTGCATTCCACCGGTCGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTa                                          <  1:401370/140‑1 (MQ=255)
                                                                                               ataaGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGAc                                     <  1:590002/140‑1 (MQ=255)
                                                                                               ataaGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGAc                                     <  1:525811/140‑1 (MQ=255)
                                                                                                           tAGTCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTGTCGCCAAGAGTACCGACCACATCCAGCGt                         >  1:362991/1‑140 (MQ=255)
                                                                                                              tCGTGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTg                      >  2:183912/1‑140 (MQ=255)
                                                                                                              tCGTGAATTTCACCGATAGTGCCTGTCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTg                      <  1:281621/140‑1 (MQ=255)
                                                                                                                 tGAATTTCACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTGCgg                   >  2:625585/1‑140 (MQ=255)
                                                                                                                        cACCGGTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTGCGGATGcaca            >  2:831438/1‑140 (MQ=255)
                                                                                                                             gTAGTGCCTGCCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTGCGGATGCACATCGGa       <  1:444585/140‑1 (MQ=255)
                                                                                                                                  gcctgcctGCTGTAAAAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTGCGGATGCACATCGGAAGCCa  <  1:278519/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TATCCGCGCCCTGTTTACCCGGCGCAGTTAACAGCACCAGCGCCATAATATCGGCGGCAACGCTGACCACAATTTTGCGTGATTTCAGTTCGCTAATAAGCGCAGTGTAGTCGTGAATTTCACCGGTAGTGCCTACCTGCTGTAACAGCACGCCGAAGACGTCCTGATGGTCGAGCACTTTTTGCGCGTCATCGACAATCACTTCAAAACCAAAGGTTTCGGCACGAGTACGGACCACATCCAGCGTTTGCGGATGCACATCGGAAGCCA  >  NC_000913/3048249‑3048518

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: