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% 135.8 / NA 41I94L (ATT→CTT) folAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (29/12);  total (29/12)

CGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGATTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGCA  >  NC_000913/50040‑50163
                                                              |                                                             
cGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCtt                                                             >  1:442624/1‑65 (MQ=255)
cGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCtt                                                             >  1:699288/1‑65 (MQ=255)
 gTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTg                                                            >  1:3040292/1‑65 (MQ=255)
        aaGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGcggcgg                                                     <  1:317840/65‑1 (MQ=255)
          gTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTc                                                   >  1:1225778/1‑65 (MQ=255)
           tCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCTGTct                                                  >  1:828129/1‑64 (MQ=39)
           tCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTcg                                                  >  1:2533715/1‑65 (MQ=255)
           tCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACAAGAAATCATGGTGCTTGGCGGCGGTcg                                                  >  1:2010118/1‑65 (MQ=255)
           tCGGTGGATGAAGCCACCGCGGCGTGTGGTGACGTACCAGAAATCATGGCGCTTGGCGGCGGCcg                                                  >  1:1366500/1‑65 (MQ=37)
                 gATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTa                                            >  1:867444/1‑65 (MQ=255)
                   tGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATg                                          >  1:2735069/1‑65 (MQ=255)
                   tGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATg                                          >  1:1909003/1‑65 (MQ=255)
                    gAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGa                                         <  1:661703/65‑1 (MQ=255)
                    gAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGa                                         <  1:1904901/65‑1 (MQ=255)
                     aaGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGGGTTTATGaa                                        >  1:685928/1‑65 (MQ=255)
                       gCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACa                                      >  1:151078/1‑65 (MQ=255)
                       gCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACa                                      <  1:2235537/65‑1 (MQ=255)
                        ccATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGACCAg                                     >  1:1025000/1‑65 (MQ=255)
                        ccATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAg                                     >  1:956078/1‑65 (MQ=255)
                        ccATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAg                                     >  1:3231490/1‑65 (MQ=255)
                        ccATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGGCTATGAAAAg                                     >  1:1804784/1‑65 (MQ=37)
                        ccATCGCGGCGTGTGGTGACGTACCAGAAATCATGGCGCTTGGCGGCGGTCGCGTTTATGAACAg                                     >  1:2535514/1‑65 (MQ=255)
                         cATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGt                                    <  1:1871153/65‑1 (MQ=255)
                         cATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGt                                    <  1:1276273/65‑1 (MQ=255)
                         cATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGt                                    <  1:2229039/65‑1 (MQ=255)
                         cATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGt                                    <  1:1446832/65‑1 (MQ=255)
                            cgcgGCGTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCt                                 <  1:18953/65‑1 (MQ=255)
                             gcggcgTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCtt                                >  1:1990762/1‑65 (MQ=255)
                                 cgTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGGACAGTTCTTGCCa                            >  1:2678892/1‑65 (MQ=39)
                                 cgTGTGGTGACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCa                            >  1:320087/1‑65 (MQ=255)
                                    gtggtgACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCaaaa                         >  1:2550462/1‑65 (MQ=255)
                                    gtggtgACGTACCAGAAATCATGGTGCTTGGCGGCGGGCGCGTTTATGAACAGTTCTTGCCaaaa                         >  1:1660290/1‑65 (MQ=39)
                                         gACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCa                    >  1:159684/1‑65 (MQ=255)
                                         gACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCa                    >  1:2186613/1‑65 (MQ=255)
                                         gACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCa                    >  1:280734/1‑65 (MQ=255)
                                         gACGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCa                    >  1:2882116/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaa                  >  1:1052952/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaa                  >  1:261108/1‑65 (MQ=255)
                                           cGTACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaa                  >  1:1482098/1‑65 (MQ=255)
                                             tACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaaaa                <  1:138957/65‑1 (MQ=255)
                                             tACCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCaaaaa                <  1:1113335/65‑1 (MQ=255)
                                              aCCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAAc               >  1:2904703/1‑65 (MQ=255)
                                              aCCAGAAATCATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAAc               >  1:564281/1‑65 (MQ=255)
                                                  gAAATCATGGTGCTTGGCGGCGGTCGCGGTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTa           >  1:2343550/1‑65 (MQ=255)
                                                       cATGGTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGa      <  1:2034249/65‑1 (MQ=255)
                                                          ggTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGc   >  1:1380060/1‑65 (MQ=255)
                                                          ggTGATTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGc   <  1:2469887/65‑1 (MQ=255)
                                                           gTGCTTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGCa  >  1:3152245/1‑65 (MQ=255)
                                                              |                                                             
CGTGGGTGAAGTCGGTGGATGAAGCCATCGCGGCGTGTGGTGACGTACCAGAAATCATGGTGATTGGCGGCGGTCGCGTTTATGAACAGTTCTTGCCAAAAGCGCAAAAACTGTATCTGACGCA  >  NC_000913/50040‑50163

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