Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,910 T→C W30R (TGG→CGG)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9100TC100.0% 138.0 / NA 38W30R (TGG→CGG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/24);  total (14/24)

CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49828‑49996
                                                                                  |                                                                                      
cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTaa                                                                                 >  2:362891/1‑90 (MQ=255)
 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTaaa                                                                                >  1:345710/1‑90 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCaa                                                                           <  1:281035/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCAGGGAAAACGCCAGGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAac                                                                          <  2:1039963/90‑1 (MQ=255)
          gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACAcc                                                                       >  2:680382/1‑90 (MQ=255)
           cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCt                                                                      <  2:970307/90‑1 (MQ=255)
           cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCt                                                                      <  1:305612/90‑1 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAac                                                                          >  2:856004/1‑85 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAac                                                                          <  1:856004/85‑1 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCtt                                                                     >  1:648114/1‑90 (MQ=255)
                  aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTtaaataa                                                               <  2:234787/90‑1 (MQ=255)
                           tCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                                      >  2:951791/1‑90 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  2:826419/1‑90 (MQ=255)
                              cgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                                   <  1:826419/90‑1 (MQ=255)
                                ttATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 <  1:456488/90‑1 (MQ=255)
                                 tATGGGCATGAAAACCCCCATGCTGGGAACCCCCCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAACCCCGTGATTATggg                                                <  2:788226/90‑1 (MQ=255)
                                        agggAAACCCCCATGCTGTGGACCCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  2:195808/88‑1 (MQ=255)
                                        aTGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  1:922037/90‑1 (MQ=255)
                                        aTGGAAAACCCCATGCTGTGGACCCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  2:385452/90‑1 (MQ=255)
                                         tgaaaaaAGCTATGCTGTGAAACCGCCCTGCCGATCTCGCCCGGTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCATa                                        <  2:1029017/87‑1 (MQ=255)
                                            aaaaCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                     <  1:282770/90‑1 (MQ=255)
                                                cGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                                 >  2:669623/1‑90 (MQ=255)
                                                   cATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTtaaataa                                                               <  1:1037260/57‑1 (MQ=255)
                                                   cATGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTtaaataa                                                               >  2:1037260/1‑57 (MQ=255)
                                                     tGCTGTGGACCCTCCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa                            <  2:792275/90‑1 (MQ=255)
                                                     tGCTGTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa                            <  1:362891/90‑1 (MQ=255)
                                                       ctgTGGACCCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc                          <  1:960783/88‑1 (MQ=255)
                                                       ctgTGGAACCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc                          <  2:345710/88‑1 (MQ=255)
                                                              aCCTGCCTGCCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc                   >  1:411673/1‑90 (MQ=255)
                                                                ctgcctgcCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg                 <  2:757383/90‑1 (MQ=255)
                                                                ctgcctgcCGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg                 >  1:601487/1‑90 (MQ=255)
                                                                 tgcctgccGATCTCGCCCGGTTTAAACGCAACGCCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt                <  2:648114/90‑1 (MQ=255)
                                                                 tgcctgccGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt                >  2:663513/1‑90 (MQ=255)
                                                                   cctgccGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                                        >  2:521363/1‑64 (MQ=255)
                                                                   cctgccGATCTCGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                                        <  1:521363/64‑1 (MQ=255)
                                                                           tctcGCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACg      >  2:113651/1‑90 (MQ=255)
                                                                               gCCCGGTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  2:140912/90‑1 (MQ=255)
                                                                               gCCCGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  1:65055/90‑1 (MQ=255)
                                                                                  |                                                                                      
CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49828‑49996

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: