Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,998,494 C→T *329* (TAG→TAA dcyD ← D‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,998,4940CT100.0% 89.0 / NA 29*329* (TAG→TAAdcyDD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/16);  total (13/16)

AAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACA  >  NC_000913/1998384‑1998626
                                                                                                              |                                                                                                                                    
aaaaaaCATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAgcg                                                                                                          <  1:834842/139‑1 (MQ=255)
           cgccgATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccg                                                                                               <  2:1030798/139‑1 (MQ=255)
                       gTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATa                                                                                   >  2:438335/1‑139 (MQ=255)
                         tGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                 <  1:184740/139‑1 (MQ=255)
                          gCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAAc                                                                                <  1:142007/139‑1 (MQ=255)
                                  taCCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                 >  1:43493/1‑130 (MQ=255)
                                  taCCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                 <  2:43493/130‑1 (MQ=255)
                                        gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTc                                                                  >  1:1048065/1‑139 (MQ=255)
                                         ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAAt                                                                           <  2:1092444/129‑1 (MQ=255)
                                         ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAAt                                                                           >  1:1092444/1‑129 (MQ=255)
                                                  gcagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAg                                                        <  1:438335/139‑1 (MQ=255)
                                                       aaCGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGttt                                                   <  1:547328/139‑1 (MQ=255)
                                                                          ccAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGttt                                                   <  1:1072687/120‑1 (MQ=255)
                                                                          ccAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGttt                                                   >  2:1072687/1‑120 (MQ=255)
                                                                                gTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAgcg                                                                                                          >  2:467436/1‑59 (MQ=255)
                                                                                gTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAgcg                                                                                                          <  1:467436/59‑1 (MQ=255)
                                                                                           tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCttt               >  1:213977/1‑139 (MQ=255)
                                                                                           tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCACtt                >  2:915990/1‑138 (MQ=255)
                                                                                           tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCACtt                <  1:915990/138‑1 (MQ=255)
                                                                                             cATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                           <  1:498088/61‑1 (MQ=255)
                                                                                             cATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                           >  2:498088/1‑61 (MQ=255)
                                                                                                gAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGgcg                                                                                                    <  1:969049/49‑1 (MQ=255)
                                                                                                gAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGgcg                                                                                                    >  2:969049/1‑49 (MQ=255)
                                                                                                gAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCgg          >  2:991251/1‑139 (MQ=255)
                                                                                                       tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                 <  2:326435/61‑1 (MQ=255)
                                                                                                       tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                 >  1:326435/1‑61 (MQ=255)
                                                                                                       tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt              >  1:57259/1‑128 (MQ=255)
                                                                                                       tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt              <  2:57259/128‑1 (MQ=255)
                                                                                                        ccGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATAcaca  <  2:61508/139‑1 (MQ=255)
                                                                                                              |                                                                                                                                    
AAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACA  >  NC_000913/1998384‑1998626

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