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

GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACC  >  NC_000913/49843‑49960
                                                              |                                                       
tttAGCGTTAGACCGCTTTATCGGCTTGAAAAACGCCATGCTGGGGACCCTGCCTGCCGATcgcg                                                       <  1:860477/64‑1 (MQ=255)
       gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGtt                                                >  1:923060/1‑65 (MQ=255)
       gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGtt                                                >  1:2614479/1‑65 (MQ=255)
           aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                            >  1:3466163/1‑65 (MQ=255)
             ctcgTTTTCCGCATGAAAACCGCCATCCTGTGGACCCTGCCTGACGATCGCGCCTGGTTTAAACg                                          <  1:973081/63‑1 (MQ=255)
              gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGGTTAAAcct                                         >  1:2901870/1‑63 (MQ=255)
              gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCCGGTTTAAACGc                                         >  1:2979699/1‑65 (MQ=255)
                  tATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAaca                                     >  1:3140217/1‑65 (MQ=255)
                   aTCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                                    <  1:1173521/65‑1 (MQ=255)
                                   ccATGCTGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                    <  1:3261204/65‑1 (MQ=255)
                                       gctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                >  1:2783518/4‑65 (MQ=255)
                                       gctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                >  1:436146/4‑65 (MQ=255)
                                        ctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt               >  1:2454342/3‑65 (MQ=255)
                                            ggAACCTGCCTGCCGATCGCGCCTGGTTTAACCGCAACACCTTAAATAAACCCGTGATTATGGgc           >  1:1315331/1‑65 (MQ=255)
                                            ggAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc           >  1:1081467/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg         >  1:934544/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg         >  1:1052839/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg         >  1:1427975/1‑65 (MQ=255)
                                              aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg         >  1:3642011/1‑65 (MQ=255)
                                               aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc        >  1:3337117/1‑65 (MQ=255)
                                               aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc        >  1:212355/1‑65 (MQ=255)
                                               aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc        >  1:2198266/1‑65 (MQ=255)
                                                cctgcctgCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc       >  1:3557072/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa      >  1:3649486/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa      >  1:2910395/1‑65 (MQ=255)
                                                 ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa      >  1:2580735/1‑65 (MQ=255)
                                                   gcctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa    >  1:2709678/1‑65 (MQ=255)
                                                     ctgccGATCCCGCCTGTTTTAACCGCAACACCTAAAAAAACCCCGTGATTATGGGCCGCCATAcc  <  1:490975/65‑1 (MQ=255)
                                                              |                                                       
GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACC  >  NC_000913/49843‑49960

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