Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 49,910 T→A W30R (TGG→AGG)  folA → dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091349,9100TA100.0% 79.5 / NA 23W30R (TGG→AGG) folAdihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (14/9);  total (14/9)

GTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49850‑49964
                                                            |                                                      
gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGtt                                                    >  1:269506/1‑65 (MQ=255)
   gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTaa                                                 <  1:12390/65‑1 (MQ=255)
    aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGGTTcaa                                                >  1:2450573/1‑65 (MQ=255)
    aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCCACCGCTCTCGCCAGGGTTaaa                                                >  1:2807118/1‑65 (MQ=255)
    aTCGCGTTATCGGCATGGAAAAAGCCATGCCGGGGAACCCGCCTACCGATCTCGCCAGGTTTaaa                                                >  1:1231428/1‑65 (MQ=255)
     tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAAc                                               >  1:2312358/1‑65 (MQ=255)
       gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGc                                             >  1:449771/1‑65 (MQ=255)
           tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAaca                                         >  1:3005932/1‑65 (MQ=255)
            aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAacac                                        <  1:300934/65‑1 (MQ=255)
            aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAacac                                        <  1:369262/65‑1 (MQ=255)
            aTCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                        <  1:554746/65‑1 (MQ=255)
             tccgCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACAcc                                       <  1:217777/62‑1 (MQ=255)
             tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACAcc                                       <  1:1684096/65‑1 (MQ=255)
             tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACAcc                                       <  1:2766016/65‑1 (MQ=255)
                gCATGGAAAACGCCATGCCGTGGACCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTta                                    <  1:1617691/65‑1 (MQ=255)
                 cATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTtaa                                   <  1:698083/65‑1 (MQ=255)
                      aaaaCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTtaaataaa                              >  1:297565/1‑65 (MQ=255)
                          cGCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCg                          >  1:1050943/1‑65 (MQ=255)
                           gCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGt                         >  1:1522890/1‑65 (MQ=255)
                           gCCATGCCGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGt                         >  1:851175/1‑65 (MQ=255)
                                 ccGTGGAACCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGCTTAt                   >  1:1426974/1‑65 (MQ=255)
                                       aaCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg             >  1:2275479/1‑65 (MQ=255)
                                        aCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            >  1:695492/1‑65 (MQ=255)
                                        aCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            >  1:2452010/1‑65 (MQ=255)
                                        aCCTGCCTACCGATCTCGCCAGGTTTAAACGCAACACCCTAAATAAACCCGTGATTATGGgccgc            >  1:1920070/1‑65 (MQ=255)
                                          ctgcctacCGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcta          >  1:1124728/1‑63 (MQ=255)
                                                accGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg    >  1:561566/2‑65 (MQ=255)
                                                  cGATCTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg  >  1:2609643/1‑65 (MQ=255)
                                                            |                                                      
GTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49850‑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: