Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 49,905 T→G L28R (CTC→CGC)  folA → dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091349,9050TG100.0% 64.9 / NA 20L28R (CTC→CGC) folAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/7);  total (13/7)

GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTG  >  NC_000913/49843‑49962
                                                              |                                                         
gTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATcgcg                                                         <  1:478240/65‑1 (MQ=255)
     cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTgg                                                    <  1:2958332/65‑1 (MQ=255)
       gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGtt                                                  >  1:1873399/1‑65 (MQ=255)
          gATCGCGTTATCGGCATGGAAAACGCCATGGTGTGGAACCTACCTGCCGATCGCGCCTGGTTTaa                                               <  1:49869/65‑1 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGGACCTGCCTGCCGATCGCGCCCGGTTTaaa                                              >  1:1843244/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:1781794/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:2054987/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:353757/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:324641/1‑65 (MQ=255)
            tCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAAc                                             <  1:2387671/65‑1 (MQ=255)
                  tATCGGCATGGAAAACGCCCTGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAaca                                       >  1:689458/1‑65 (MQ=255)
                   gtCGTCATGGAAAGCTCCTTGGTGTTGAACCGGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                                      <  1:1471133/64‑1 (MQ=255)
                   aTCGGCATGGAAAACCCCATGCTGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                                      <  1:939878/65‑1 (MQ=255)
                      ggCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCtt                                   >  1:971769/1‑65 (MQ=255)
                        cATGGAAACCGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaa                                 <  1:2204467/65‑1 (MQ=255)
                              aaaCGCCAGGCCGTGGAACCTGCCTGCCGACCGCGCCAGGTTTAAACGCAACACCTTAAATAAAc                           <  1:1461728/65‑1 (MQ=255)
                                       gctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                  >  1:2755399/4‑65 (MQ=255)
                                       gctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                  >  1:2176000/4‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg           >  1:1554292/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg           >  1:943110/1‑65 (MQ=255)
                                               aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc          >  1:2781491/1‑65 (MQ=255)
                                               aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc          >  1:2332426/1‑65 (MQ=255)
                                                cctgcctgCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc         >  1:730810/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa        >  1:2602761/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa        >  1:2006626/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa        >  1:617289/1‑65 (MQ=255)
                                                    cctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAc     >  1:1154722/1‑65 (MQ=255)
                                                       gccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg  >  1:1908852/1‑65 (MQ=255)
                                                              |                                                         
GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTG  >  NC_000913/49843‑49962

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