Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 369,076 A→G K156E (AAA→GAA)  mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913369,0760AG100.0% 107.1 / NA 31K156E (AAA→GAA) mhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/12);  total (19/12)

AACGCCTTTATTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGA  >  NC_000913/368944‑369202
                                                                                                                                    |                                                                                                                              
aaCGCCTTTATTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACgg                                                                                                                          >  1:543520/1‑139 (MQ=255)
         aTTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGc                                                                                                                 >  1:507897/1‑139 (MQ=255)
               ccGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCggg                                                                                                           >  1:696107/1‑139 (MQ=255)
                    ggTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATa                                                                                                      >  2:866571/1‑139 (MQ=255)
                                       gAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTgg                                                                                   <  2:309364/139‑1 (MQ=255)
                                         aGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTa                                                                                 >  1:215531/1‑139 (MQ=255)
                                         aGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTa                                                                                 >  2:22524/1‑139 (MQ=255)
                                            ggTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGcc                                                                              >  2:177441/1‑139 (MQ=255)
                                             gtgtCGCGTTTTACGAATGTGCGCTGATTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGAGCAGCGGGAAATAGACAAAGCCCAGTGGCTGGTAGCCt                                                                             >  2:479945/1‑139 (MQ=255)
                                               gtCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCtgt                                                                           <  2:543520/139‑1 (MQ=255)
                                                tCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCtgtg                                                                          >  1:177443/1‑139 (MQ=255)
                                                     ttttCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGt                                                                     >  1:21410/1‑139 (MQ=255)
                                                         ccGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAg                                                                 <  1:866571/139‑1 (MQ=255)
                                                          cGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGc                                                                <  1:177441/139‑1 (MQ=255)
                                                                  cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTg                                                        >  1:92203/1‑139 (MQ=255)
                                                                  cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTg                                                        >  1:572677/1‑139 (MQ=255)
                                                                      gctTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCg                                                    >  1:854167/1‑139 (MQ=255)
                                                                            tttCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCAc                                              <  2:21410/139‑1 (MQ=255)
                                                                                 cgcgAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGa                                                                                                                      <  1:459701/62‑1 (MQ=255)
                                                                                 cgcgAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGa                                                                                                                      >  2:459701/1‑62 (MQ=255)
                                                                                       cTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCa                                                                                                  <  1:316656/76‑1 (MQ=255)
                                                                                       cTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCa                                                                                                  >  2:316656/1‑76 (MQ=255)
                                                                                          gAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGc                                                                >  2:292369/1‑107 (MQ=255)
                                                                                          gAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGc                                                                <  1:292369/107‑1 (MQ=255)
                                                                                                ttCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAg                          >  1:969853/1‑139 (MQ=255)
                                                                                                                        aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGt                                                                                                    <  2:968708/41‑1 (MQ=255)
                                                                                                                        aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGt                                                                                                    >  1:968708/1‑41 (MQ=255)
                                                                                                                        aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa  <  2:507897/139‑1 (MQ=255)
                                                                                                                        aCCTTGCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGa  <  2:969853/139‑1 (MQ=255)
                                                                                                                             gCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAg                                                                                          >  2:813196/1‑46 (MQ=255)
                                                                                                                             gCACCTGGAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAg                                                                                          <  1:813196/46‑1 (MQ=255)
                                                                                                                                    |                                                                                                                              
AACGCCTTTATTCAGCCGCAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGA  >  NC_000913/368944‑369202

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