Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 50,102 A→C I94L (ATT→CTT)  folA → dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091350,1020AC100.0% 116.3 / NA 35I94L (ATT→CTT) folAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (24/11);  total (24/11)

CGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGATTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGC  >  NC_000913/50040‑50162
                                                              |                                                            
cGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCtt                                                            >  1:3328228/1‑65 (MQ=255)
cGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCtt                                                            >  1:3053308/1‑65 (MQ=255)
cGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCtt                                                            >  1:3355649/1‑65 (MQ=255)
cGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCtt                                                            >  1:2065049/1‑65 (MQ=255)
cGTGGGTGAAGTCGGGGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATAGTGCtt                                                            >  1:2164250/1‑65 (MQ=255)
  tGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTgg                                                          >  1:3125216/1‑65 (MQ=255)
  tGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTgg                                                          >  1:3342911/1‑65 (MQ=255)
       gAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGgcggcg                                                     >  1:3177677/1‑65 (MQ=255)
           tCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGcggcggctg                                                 >  1:2123858/1‑62 (MQ=255)
                 gATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGCCGCGTTTa                                           >  1:2088465/1‑65 (MQ=255)
                  aTGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTAt                                          >  1:3069182/1‑65 (MQ=255)
                    gAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGa                                        <  1:1514394/65‑1 (MQ=255)
                    gAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGa                                        <  1:816867/65‑1 (MQ=255)
                       gCCATCGCGGCGCGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGaata                                     >  1:2629220/1‑63 (MQ=38)
                        ccATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAg                                    >  1:1623743/1‑65 (MQ=255)
                         cATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGt                                   <  1:1752906/65‑1 (MQ=255)
                            cgcgGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCt                                <  1:3052655/65‑1 (MQ=255)
                                 cgTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCa                           >  1:2851583/1‑65 (MQ=255)
                                         gACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCa                   >  1:1464903/1‑65 (MQ=255)
                                         gACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCa                   >  1:3043795/1‑65 (MQ=255)
                                         gACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCa                   <  1:1457513/65‑1 (MQ=255)
                                          aCGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaa                  >  1:725323/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaa                 >  1:1517895/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaa                 >  1:591318/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaa                 >  1:2139063/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaa                 >  1:1482336/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAAAAGTTCTTGCCAAAAGCGCaaa                 >  1:1970303/1‑65 (MQ=39)
                                             tACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaaaa               <  1:2748552/65‑1 (MQ=255)
                                             tACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaaaa               <  1:51811/65‑1 (MQ=255)
                                             tACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaaaa               <  1:1058621/65‑1 (MQ=255)
                                               ccAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACt             <  1:3298182/65‑1 (MQ=255)
                                                cAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTg            >  1:2749597/1‑65 (MQ=255)
                                                 aGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGt           <  1:934638/65‑1 (MQ=255)
                                                    aaTCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATc        >  1:3298629/1‑65 (MQ=255)
                                                      tCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTg      >  1:52867/1‑65 (MQ=255)
                                                       cATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGa     <  1:2764148/65‑1 (MQ=255)
                                                          ggTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGc  >  1:34996/1‑65 (MQ=255)
                                                          ggTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGc  >  1:1147460/1‑65 (MQ=255)
                                                          ggTGCTTGGCGGCGGTCGCGTTCGTGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACgt  >  1:387920/1‑64 (MQ=255)
                                                              |                                                            
CGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGATTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGC  >  NC_000913/50040‑50162

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