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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,8840CT100.0% 157.2 / NA 48P21L (CCG→CTG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/27);  total (21/27)

TCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCA  >  NZ_CP009273/49820‑49969
                                                                |                                                                                     
tcAATGATCAGTCTGATTGCGGCGTTGGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                              <  2:706305/90‑1 (MQ=255)
taaaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                              <  2:837070/88‑1 (MQ=255)
 cAATGATCAGTCTGATTGCGGCTTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCt                                                             <  2:1099471/90‑1 (MQ=255)
  aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTg                                                            <  2:10527/90‑1 (MQ=255)
       tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGGCGATCTCGCCTGGTTTa                                                       >  1:1339137/1‑90 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGa                                                                                  >  2:995828/1‑62 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGa                                                                                  <  1:995828/62‑1 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                      <  1:863125/90‑1 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                      >  1:872599/1‑90 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                      >  1:1052442/1‑90 (MQ=255)
         aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                     >  2:219179/1‑90 (MQ=255)
         aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                     >  1:1164335/1‑90 (MQ=255)
         aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                     >  2:893805/1‑90 (MQ=255)
          gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                                    >  1:345251/1‑90 (MQ=255)
           tCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                      >  1:183134/1‑87 (MQ=255)
           tCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                      <  2:183134/87‑1 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                <  2:493590/90‑1 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                <  2:1052442/90‑1 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                <  1:1293330/90‑1 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                >  1:368580/1‑90 (MQ=255)
               aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                               <  2:1122234/90‑1 (MQ=255)
                   cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                           <  2:943274/90‑1 (MQ=255)
                   cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                           <  2:345251/90‑1 (MQ=255)
                   cggcgTTAGCGGTAGACCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                           <  1:1164157/90‑1 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                               >  1:39461/1‑85 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                               <  2:39461/85‑1 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                          >  2:273677/1‑90 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                          >  1:1051487/1‑90 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                          >  2:984019/1‑90 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                           <  1:879318/89‑1 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                           >  2:879318/1‑89 (MQ=255)
                     gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTta                                         <  2:616922/90‑1 (MQ=255)
                          aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                    <  1:1013836/90‑1 (MQ=255)
                          aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                    <  2:1051487/90‑1 (MQ=255)
                              gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                <  2:1218173/90‑1 (MQ=255)
                                aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                              <  2:35192/90‑1 (MQ=255)
                                aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                              <  2:872599/90‑1 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:700144/1‑90 (MQ=255)
                                        ttATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                      >  1:1345550/1‑90 (MQ=255)
                                        ttATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTCCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                      <  1:732449/90‑1 (MQ=255)
                                            cGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                  <  1:1206309/90‑1 (MQ=255)
                                                     aaaCCCCATGCTGTGGACCCTCCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg         <  2:1216069/90‑1 (MQ=255)
                                                      aaCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTgg        >  1:1164499/1‑90 (MQ=255)
                                                        cGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa      >  2:135490/1‑90 (MQ=255)
                                                        cGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa      <  1:273677/90‑1 (MQ=255)
                                                            aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                               >  2:577613/1‑45 (MQ=255)
                                                            aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                               <  1:577613/45‑1 (MQ=255)
                                                            aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca  >  1:225668/1‑90 (MQ=255)
                                                                |                                                                                     
TCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCA  >  NZ_CP009273/49820‑49969

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