Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 49,884 C→T P21L (CCG→CTG)  folA → dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091349,8840CT100.0% 51.7 / NA 16P21L (CCG→CTG) folAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/7);  total (9/7)

AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAC  >  NC_000913/49822‑49937
                                                              |                                                     
aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGt                                                     <  1:198840/65‑1 (MQ=255)
 aTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCtgtg                                                    <  1:2022194/65‑1 (MQ=255)
   gATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGa                                                  <  1:856970/65‑1 (MQ=255)
   gATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGa                                                  >  1:1550616/1‑65 (MQ=255)
   gATCAGTCTGATTGCGGCGTGAGCGGTAGATCGCGTTATCGGGATGTAAAACGCCATGCTGTGGa                                                  <  1:83341/65‑1 (MQ=255)
     tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAc                                                <  1:705481/65‑1 (MQ=255)
      cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGGAAACGCCATGCCGTGGGAcc                                               >  1:240005/1‑65 (MQ=255)
        gTCTGATGGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAcctg                                             <  1:2848713/65‑1 (MQ=255)
        gTCGGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAcctg                                             <  1:242932/65‑1 (MQ=255)
          cTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCCTGCTGGGGAAcgtgcc                                           >  1:786605/1‑65 (MQ=255)
               tGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCg                                      >  1:72285/1‑65 (MQ=255)
                gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATTCTGCGGAACCTCCCTGCCGa                                     >  1:731953/1‑65 (MQ=255)
                     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:2054987/1‑65 (MQ=255)
                                aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                     >  1:324641/1‑65 (MQ=255)
                                aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                     >  1:353757/1‑65 (MQ=255)
                                aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                     >  1:1781794/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)
                                                              |                                                     
AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAC  >  NC_000913/49822‑49937

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