Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,889,885 T→C Y386H (TAT→CAT)  tnaA → tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,889,8850TC100.0% 166.9 / NA 44Y386H (TAT→CAT) tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/25);  total (19/25)

ACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTCAAACACATATGGACTTC  >  NC_000913/3889743‑3890031
                                                                                                                                              |                                                                                                                                                  
aCAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGt                                                                                                                                             <  2:278396/150‑1 (MQ=255)
aCAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGt                                                                                                                                             <  1:507992/150‑1 (MQ=255)
aCAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTc                                                                                                                                            <  2:254713/151‑1 (MQ=255)
    tATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCg                                                                                                                                        >  2:629134/1‑151 (MQ=255)
     aTCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCg                                                                                                                                        <  1:475416/150‑1 (MQ=255)
       cTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGt                                                                                                                                      >  1:528970/1‑150 (MQ=255)
            cGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGt                                                                                                                                >  2:473400/1‑151 (MQ=255)
            cGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGt                                                                                                                                >  1:498028/1‑151 (MQ=255)
              aTGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGc                                                                                                                              >  2:99396/1‑151 (MQ=255)
                     ggAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGa                                                                                                                        >  2:264570/1‑150 (MQ=255)
                      gAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGaaa                                                                                                                      <  1:272167/151‑1 (MQ=255)
                       aaGAGATTGGCGTCGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAAt                                                                                                                     >  1:313919/1‑151 (MQ=255)
                        agagATTGGCGTTGTCTTCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAAtt                                                                                                                    >  1:390974/1‑151 (MQ=255)
                                gCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCttt                                                                                                            <  1:362843/151‑1 (MQ=255)
                                         gCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTa                                                                                                     <  2:621471/149‑1 (MQ=255)
                                          ccaccaGGCTGGCGGTCACGCGGCATTCGTTTATGCCTGTAAACTGTTGCCGCATATCCCTCCAGACCAGTTCCCTGCACATTCACTATCCTTCGAGCTGCATAAAGTCGCCCTTATCCGTTCGCTAGAAATTCTCTCTTTCCTTTTAgac                                                                                                  <  2:124776/151‑3 (MQ=255)
                                           cagcagGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGcc                                                                                                 <  2:117108/151‑1 (MQ=255)
                                               agGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCcgcg                                                                                              >  1:199921/1‑150 (MQ=255)
                                               agGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGa                                                                                             <  2:508601/151‑1 (MQ=255)
                                                   gggcggTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCg                                                                                         >  2:90454/1‑151 (MQ=255)
                                                          cACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAAcc                                                                                   <  2:281826/150‑1 (MQ=255)
                                                                gCATTCGTTGATGCCGGTAAACTGTTGCCGCGTATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAaac                                                                            <  1:571433/151‑1 (MQ=255)
                                                                  aTTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAaacaa                                                                          <  2:110288/151‑1 (MQ=255)
                                                                  aTTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAaacaa                                                                          <  1:629153/151‑1 (MQ=255)
                                                                    tCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCGGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAaacaac                                                                         <  1:561655/150‑1 (MQ=255)
                                                                         gATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCAt                                                                   >  2:511703/1‑151 (MQ=255)
                                                                             ccGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGccc                                                               <  2:98808/151‑1 (MQ=255)
                                                                             ccGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGccc                                                               >  2:275247/1‑151 (MQ=255)
                                                                                 tAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCt                                                           >  2:226689/1‑151 (MQ=255)
                                                                                            cGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGGTCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCg                                                >  2:462636/1‑151 (MQ=255)
                                                                                             gCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGt                                               <  2:136724/151‑1 (MQ=255)
                                                                                               atatCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGttt                                             <  2:261312/151‑1 (MQ=255)
                                                                                                      ggCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCAtt                                      <  1:643671/151‑1 (MQ=255)
                                                                                                        cAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTcc                                    <  2:618197/151‑1 (MQ=255)
                                                                                                             cAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCcgcgcg                               >  2:468023/1‑151 (MQ=255)
                                                                                                                   ccGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACATCTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACat                         <  2:159732/151‑1 (MQ=255)
                                                                                                                   ccGGCACAGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACat                         >  1:239244/1‑151 (MQ=255)
                                                                                                                         caGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACt                    <  2:301048/150‑1 (MQ=255)
                                                                                                                         caGGCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTc                   <  1:560406/151‑1 (MQ=255)
                                                                                                                           ggCGCTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTCa                  >  1:509802/1‑150 (MQ=255)
                                                                                                                             cgcTGGCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTCAAac               >  2:562673/1‑151 (MQ=255)
                                                                                                                                 ggCCTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTCAAacac             <  2:179793/149‑1 (MQ=255)
                                                                                                                                   ccTGCGAGCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTCAAACACatat         <  1:110796/151‑1 (MQ=255)
                                                                                                                                          gCTGCATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTCAAACACATATGGACTTc  >  1:154958/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
ACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTGAACTGCTGCGTTTAACCATTCCGCGCGCAACATATACTCAAACACATATGGACTTC  >  NC_000913/3889743‑3890031

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑