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% 123.8 / NA 36L28R (CTC→CGC) folAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (29/7);  total (29/7)

GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49843‑49964
                                                              |                                                           
tttAGCGGTAGACCGCGTTATCGGCATGGAAACCCCCATGCTGTGGAACCTGCCTGCCGATcgcg                                                           <  1:1004693/64‑1 (MQ=255)
gTTAGCGGTAGATCGCGTTATCGGCATGGAAACCGCCATGCTGTGGACCCTGCCTGCCGATcgcg                                                           <  1:3866678/65‑1 (MQ=255)
gTTAGCGGTAGATCGCGTTAGCGGCATGGAAAACGCCTTGCTGTGGAACCTGCCTGCCGATcgcg                                                           <  1:796317/65‑1 (MQ=255)
      gctAGACCGGGTTAGGGGAATGGAAAACGCGTTGCGGTGGAGCCTGCCTGCCGATCGCGCCTGGt                                                     <  1:2693585/63‑1 (MQ=255)
       gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGtt                                                    >  1:817160/1‑65 (MQ=255)
         aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTa                                                  >  1:1223018/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                                >  1:2104965/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                                >  1:3109264/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                                >  1:2865174/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                                >  1:285277/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                                >  1:487093/1‑65 (MQ=255)
             cgcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACg                                              >  1:3237357/1‑65 (MQ=255)
                gTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACgtct                                           >  1:2375853/1‑62 (MQ=255)
                   aTCGGCATGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                                        <  1:461149/65‑1 (MQ=255)
                   aTCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                                        <  1:2237580/65‑1 (MQ=255)
                    tCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACAcc                                       <  1:3231433/65‑1 (MQ=255)
                      ggCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCtt                                     >  1:3378340/1‑65 (MQ=255)
                        cATGGAAACCGCCATGCTTTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaa                                   <  1:13646/65‑1 (MQ=255)
                        cATGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaa                                   <  1:2793909/65‑1 (MQ=255)
                             aaaaCGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAAACCCTtcaataac                              >  1:3683377/1‑64 (MQ=255)
                              aaccGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAAc                             <  1:1322103/62‑1 (MQ=255)
                                  gCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                         >  1:1613565/1‑65 (MQ=255)
                                  gCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                         >  1:3185976/1‑65 (MQ=255)
                                   ccATGCTGTGGAACCTGCCTGCCGATCGCCCCTGGTTTAAACGCAACACCTTACATAAACTCGTg                        >  1:1795755/1‑65 (MQ=255)
                                       gctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                    >  1:3766411/4‑65 (MQ=255)
                                       gctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                    >  1:173269/4‑65 (MQ=255)
                                        ctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                   >  1:1523883/3‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg             >  1:2785181/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg             >  1:3202017/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg             >  1:222320/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGCGATTATGGgccg             >  1:2533160/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACACAACACCTTAAATAAACCCGTGATTATGGgccg             >  1:951267/1‑65 (MQ=255)
                                               aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            >  1:3949756/1‑65 (MQ=255)
                                               aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            >  1:1686477/1‑65 (MQ=255)
                                               aCCTGCATGCCGATCGCGCCTGGGTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            >  1:3008159/1‑65 (MQ=255)
                                                cctgcctgCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc           >  1:3228689/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa          >  1:1974720/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa          >  1:3901814/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa          >  1:1571629/1‑65 (MQ=255)
                                                     cggccGATCGCGCTTGTTTAAAACGCAACCCTAAAAATAACCCCGTGATTATGGGCCGCCATAcc      <  1:1829757/63‑1 (MQ=255)
                                                       gccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg    >  1:2605879/1‑65 (MQ=255)
                                                       gccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg    >  1:1896092/1‑65 (MQ=255)
                                                         cGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg  >  1:2075513/1‑65 (MQ=255)
                                                              |                                                           
GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49843‑49964

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