Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,826,546 T→C E126E (GAA→GAG astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,826,5460TC100.0% 101.8 / NA 29E126E (GAA→GAGastEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/15);  total (14/15)

TCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAC  >  NC_000913/1826418‑1826677
                                                                                                                                |                                                                                                                                   
tCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTcgc                                                                                                                           >  1:109848/1‑139 (MQ=255)
                      gTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCt                                                                                                     <  2:109848/139‑1 (MQ=255)
                                 aaCTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAg                                                                                         >  1:957277/1‑139 (MQ=255)
                                    tGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCt                                                                                       >  2:813853/1‑139 (MQ=255)
                                                 aGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGcc                                                                          <  2:342898/139‑1 (MQ=255)
                                                     aCGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCg                                                                                                                    >  2:892564/1‑93 (MQ=255)
                                                     aCGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCg                                                                                                                    <  1:892564/93‑1 (MQ=255)
                                                                       aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGa                                                                               >  1:137242/1‑112 (MQ=255)
                                                                       aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGa                                                                               <  2:137242/112‑1 (MQ=255)
                                                                            ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCt                                                                                                     >  2:361345/1‑85 (MQ=255)
                                                                            ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCt                                                                                                     <  1:361345/85‑1 (MQ=255)
                                                                            ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                               >  1:204827/1‑139 (MQ=255)
                                                                            ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                               >  1:1047747/1‑139 (MQ=255)
                                                                                ccAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTa                                                    >  2:697607/1‑130 (MQ=255)
                                                                                ccAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTa                                                    <  1:697607/130‑1 (MQ=255)
                                                                                 cAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTcgcgcg                                                                                                                        >  2:1153599/1‑61 (MQ=255)
                                                                                 cAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTcgcgcg                                                                                                                        <  1:1153599/61‑1 (MQ=255)
                                                                                        caGATTCTTTGCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCt                                  <  2:957277/139‑1 (MQ=255)
                                                                                                     cTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCa                      >  1:803101/1‑139 (MQ=255)
                                                                                                      tGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCag                     <  2:1047747/139‑1 (MQ=255)
                                                                                                       ggTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGATAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCagg                    >  2:915265/1‑139 (MQ=255)
                                                                                                         tCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCaggag                  <  1:813853/139‑1 (MQ=255)
                                                                                                                  tCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAAGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACAtt                                                                   <  2:1062510/81‑1 (MQ=255)
                                                                                                                  tCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAAGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACAtt                                                                   >  1:1062510/1‑81 (MQ=255)
                                                                                                                   cTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg        <  2:803101/139‑1 (MQ=255)
                                                                                                                      ccAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGAt     <  2:886760/139‑1 (MQ=255)
                                                                                                                         ggCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAc  <  2:1103025/139‑1 (MQ=255)
                                                                                                                            aCTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACg                                                                                >  2:379418/1‑58 (MQ=255)
                                                                                                                            aCTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACg                                                                                <  1:379418/58‑1 (MQ=255)
                                                                                                                                |                                                                                                                                   
TCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCAC  >  NC_000913/1826418‑1826677

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