Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 371,379 T→C G47G (GGT→GGC mhpC → 2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913371,3790TC100.0% 132.9 / NA 39G47G (GGT→GGCmhpC2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/24);  total (15/24)

TTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAA  >  NC_000913/371244‑371513
                                                                                                                                       |                                                                                                                                      
ttATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTa                                                                                                                                    >  1:783980/1‑140 (MQ=255)
      gCCACAAACCGAAGCAGCCACCCGCCGTTATCTGAATCTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTActggct                                                                                                                              <  1:279681/140‑1 (MQ=255)
                 aaGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTc                                                                                                                   <  2:783980/140‑1 (MQ=255)
                      gccACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCg                                                                                                              <  1:151999/140‑1 (MQ=255)
                                  tttCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATcc                                                                                                  <  2:812794/140‑1 (MQ=255)
                                    tCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCTGCCCGGGCCCTAATGGGTGGGCGAAGTTCAGCCGCAATATCGATCCGc                                                                                                >  2:220136/1‑140 (MQ=255)
                                        aaTGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTTGCTGGGCCAACTTCAGCCGCAATATCGATCCGCTGGt                                                                                            >  1:288212/1‑140 (MQ=255)
                                                    gCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTa                                                                                >  2:764070/1‑140 (MQ=255)
                                                            aaCGCTGCGCATCCATTTTAATTACTGCGGAAAAGTCGATGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGa                                                                        <  2:288212/140‑1 (MQ=255)
                                                            aaCGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGa                                                                        <  1:764070/140‑1 (MQ=255)
                                                                 tGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCctg                                                                   <  1:658432/140‑1 (MQ=255)
                                                                           ttttAATGACTGCGGTCAAGGCTACCAAATCGTTGTCATGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTc                                                         <  1:220136/140‑1 (MQ=255)
                                                                           ttttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTc                                                         >  1:135574/1‑140 (MQ=255)
                                                                             ttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATTGATCCGCTGGTATAGGGGGGCTATCGGGGGGTCCTGCTGGATTGTCCg                                                       >  2:255745/1‑140 (MQ=255)
                                                                             ttAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCg                                                       >  1:501209/1‑140 (MQ=255)
                                                                                     tGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGGGATCCTGCTGGATTGTCCGGGTTgggg                                               >  1:389/1‑140 (MQ=255)
                                                                                     tGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCAGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGt                                                    <  1:90602/135‑1 (MQ=255)
                                                                                     tGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCAGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGt                                                    >  2:90602/1‑135 (MQ=255)
                                                                                      gCGGACAAGGCGACGCAACCGTCGTTCTGATGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGc                                              <  1:234115/140‑1 (MQ=255)
                                                                                      gCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGTCTATCGGATGATCCAGCTGGACTGTCCGGGTTGGGGc                                              >  1:370753/1‑140 (MQ=255)
                                                                                      gCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGc                                              >  2:661621/1‑140 (MQ=255)
                                                                                       cGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCAGCTGGTAGAGGCGGGCTTTCGGGAGATTCTGTTGGTTTGGCCGGGGTGGGGCa                                             >  2:34290/1‑140 (MQ=255)
                                                                                        ggACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCaa                                            <  2:446343/140‑1 (MQ=255)
                                                                                                cgacgaAACCGTTGTTCTGCTGCATCGTTCCGGCCCGTGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGAtt                                    <  1:255745/140‑1 (MQ=255)
                                                                                                cgacgaAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGAtt                                    <  1:661621/140‑1 (MQ=255)
                                                                                                  acgaAACCGTTGTCCTGTTGCATGGTTCCGGTCCGGGCGCTACTGGCTTGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCg                                  <  2:284746/140‑1 (MQ=255)
                                                                                                  acgaAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCg                                  <  2:135574/140‑1 (MQ=255)
                                                                                                  acgaAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGAACCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCg                                  <  1:554540/140‑1 (MQ=255)
                                                                                                  acgaAACCGTTGTCCTGCTGCAAGGTTCCGGACCGTGCGCTACTGGCTTGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCg                                  <  2:263466/140‑1 (MQ=255)
                                                                                                   cgaAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCgg                                 <  2:389/140‑1 (MQ=255)
                                                                                                   cgaAACCGTTGTCCTGCGGCATTGTTCCGGCTCGGCCGCCACTGGCTGAGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCgg                                 <  2:370753/140‑1 (MQ=255)
                                                                                                    gaAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGt                                <  1:34290/140‑1 (MQ=255)
                                                                                                     aaaCCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTc                               <  1:103795/140‑1 (MQ=255)
                                                                                                      aaCCGTTGTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCg                              <  2:501209/140‑1 (MQ=255)
                                                                                                             gTCCTGCTGCATGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAg                       >  2:200198/1‑140 (MQ=255)
                                                                                                                        tGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGAt            <  2:594447/140‑1 (MQ=255)
                                                                                                                        tGGTTCCGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGAt            >  1:594447/1‑140 (MQ=255)
                                                                                                                             ccGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGAt       >  2:205169/1‑140 (MQ=255)
                                                                                                                             ccGGCCCGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGAt       >  1:642848/1‑140 (MQ=255)
                                                                                                                                  ccGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTaa  <  2:642848/140‑1 (MQ=255)
                                                                                                                                  ccGGGCGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTaa  <  1:205169/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
TTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAA  >  NC_000913/371244‑371513

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: