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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9100TG96.2% 160.4 / ‑1.9 52W30G (TGG→GGG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (20/30);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49825‑49996
                                                                                     |                                                                                      
gATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATAGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGtt                                                                                    <  2:511731/90‑1 (MQ=255)
   cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTaa                                                                                 >  2:968890/1‑90 (MQ=255)
   cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGGTTaa                                                                                 >  2:548706/1‑90 (MQ=255)
     gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAAc                                                                               >  2:710007/1‑90 (MQ=255)
         gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                           <  1:108514/90‑1 (MQ=255)
         gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                           <  1:218462/90‑1 (MQ=255)
         gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                           <  1:631320/90‑1 (MQ=255)
         gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                           <  1:877082/90‑1 (MQ=255)
         gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGAAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                                                           <  1:464139/90‑1 (MQ=255)
          aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac                                                                          <  2:52582/90‑1 (MQ=255)
          aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac                                                                          <  2:397805/90‑1 (MQ=255)
          aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac                                                                          <  2:297961/90‑1 (MQ=255)
          aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac                                                                          <  1:777775/90‑1 (MQ=255)
             gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                                                       >  1:699087/1‑90 (MQ=255)
              cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCt                                                                      >  1:146660/1‑90 (MQ=255)
              cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCt                                                                      <  2:724350/90‑1 (MQ=255)
               ggcgTTAGCGGTGGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt                                                                     >  1:612585/1‑90 (MQ=255)
               ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac                                                                          >  1:548652/1‑85 (MQ=255)
               ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac                                                                          <  2:548652/85‑1 (MQ=255)
               ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt                                                                     >  2:381839/1‑90 (MQ=255)
               ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt                                                                     >  2:771039/1‑90 (MQ=255)
               ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt                                                                     >  1:215832/1‑90 (MQ=255)
               ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt                                                                     >  2:509476/1‑90 (MQ=255)
                   ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaat                                                                 <  2:839534/90‑1 (MQ=255)
                   ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaat                                                                 >  2:312763/1‑90 (MQ=255)
                    tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaata                                                                <  1:777832/90‑1 (MQ=255)
                           aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGt                                                         <  2:639490/90‑1 (MQ=255)
                                gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCCGCCGATCTCGCCGGGGTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:668038/1‑90 (MQ=255)
                                   ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                                                       <  1:243368/68‑1 (MQ=255)
                                   ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                                                       >  2:243368/1‑68 (MQ=255)
                                   ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 >  2:133621/1‑90 (MQ=255)
                                   ttATCGGAATGAAAACCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 <  1:392859/90‑1 (MQ=255)
                                    tATCGGCATGAAAACCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                <  2:651042/90‑1 (MQ=255)
                                           aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  2:668038/90‑1 (MQ=255)
                                           aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  2:273945/90‑1 (MQ=255)
                                              gAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAcc                                      <  2:359473/90‑1 (MQ=255)
                                                aaaCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                                    <  2:682015/90‑1 (MQ=255)
                                                   cGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                                 >  2:153902/1‑90 (MQ=255)
                                                       aTGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaataaa                                                              <  2:15445/57‑1 (MQ=255)
                                                       aTGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaataaa                                                              >  1:15445/1‑57 (MQ=255)
                                                       aTGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                             >  1:365652/1‑90 (MQ=255)
                                                          ccGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc                          <  2:253412/90‑1 (MQ=255)
                                                              ggAACCTCCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  1:787529/90‑1 (MQ=255)
                                                              ggAACCTCCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  2:918059/90‑1 (MQ=255)
                                                                   ctgcctgcCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg                 >  2:306722/1‑90 (MQ=255)
                                                                    tgcctgccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt                >  1:781703/1‑90 (MQ=255)
                                                                             aTCTCGCCTGGTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc       <  1:949611/90‑1 (MQ=255)
                                                                             aTCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc       <  2:16593/90‑1 (MQ=255)
                                                                             aTCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc       <  1:92676/90‑1 (MQ=255)
                                                                                tcGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCa    <  1:710007/90‑1 (MQ=255)
                                                                                  gCCGGGTTTAACCGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  1:463530/90‑1 (MQ=255)
                                                                                  gCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  1:147111/90‑1 (MQ=255)
                                                                                     |                                                                                      
GATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49825‑49996

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