Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,898 G→A A26T (GCC→ACC)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,8980GA97.6% 134.8 / NA 41A26T (GCC→ACC) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/25);  minor base C (0/1);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGT  >  NZ_CP009273/49815‑49982
                                                                                   |                                                                                    
aaaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGAtc                                                                                <  1:1622/90‑1 (MQ=255)
     tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGcc                                                                           <  2:1074778/90‑1 (MQ=255)
     tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGcc                                                                           <  1:319721/90‑1 (MQ=255)
      cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGccg                                                                          <  2:1013036/90‑2 (MQ=255)
            tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTATCGATCTCGCCGGGTTTa                                                                    >  1:796507/1‑90 (MQ=255)
            tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGGTTa                                                                    >  1:58645/1‑90 (MQ=255)
              aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTaaa                                                                  >  2:682856/1‑90 (MQ=255)
              aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTACCGATCTCGCCGGGTTTaaa                                                                  >  1:335151/1‑90 (MQ=255)
               gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAAc                                                                 <  2:412990/90‑1 (MQ=255)
                   gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCaa                                                             <  1:398258/90‑1 (MQ=255)
                   gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCaa                                                             <  2:1037934/90‑1 (MQ=255)
                   gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCaa                                                             <  1:137058/90‑1 (MQ=255)
                   gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCaa                                                             <  2:83193/90‑1 (MQ=255)
                    aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAac                                                            <  2:231480/90‑1 (MQ=255)
                        cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCt                                                        <  1:325731/90‑1 (MQ=255)
                        cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCt                                                        <  1:18715/90‑1 (MQ=255)
                        cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCt                                                        <  1:621108/90‑1 (MQ=255)
                         ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCtt                                                       >  2:392493/1‑90 (MQ=255)
                         ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCtt                                                       >  2:991183/1‑90 (MQ=255)
                         ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCtt                                                       >  2:26222/1‑90 (MQ=255)
                             ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTtaaat                                                   <  1:1115149/90‑1 (MQ=255)
                               aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTtaaataa                                                 <  2:1024189/90‑1 (MQ=255)
                                    tAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCg                                            >  2:1032567/1‑90 (MQ=255)
                                    tAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCg                                            <  1:154949/90‑1 (MQ=255)
                                          gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAac                                                            >  1:695555/1‑68 (MQ=255)
                                          gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAac                                                            <  2:695555/68‑1 (MQ=255)
                                          gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCt                                                        <  1:84702/72‑1 (MQ=255)
                                          gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCt                                                        >  2:84702/1‑72 (MQ=255)
                                           cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                     <  1:829303/90‑1 (MQ=255)
                                             ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                   <  2:363623/90‑1 (MQ=255)
                                               aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                 >  2:833535/1‑90 (MQ=255)
                                                tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCCATCTCGCCGGGGTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                                >  1:306079/1‑90 (MQ=255)
                                                     aTGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                           <  1:639531/90‑1 (MQ=255)
                                                      tGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                          <  2:77981/90‑1 (MQ=255)
                                                         aaaaCGCCATGCCGTGGACCCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAACCCCGTGATTATGGGCCGCCATACCt                       <  1:639976/90‑1 (MQ=255)
                                                         aaaaCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                       <  1:450391/90‑1 (MQ=255)
                                                         aaaaCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                       <  2:883434/90‑1 (MQ=255)
                                                                 aTGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGAagtca               >  1:164031/1‑90 (MQ=255)
                                                                  tGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa              >  2:724111/1‑90 (MQ=255)
                                                                       tGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt         >  2:960048/1‑90 (MQ=255)
                                                                        ggAACCTGCCTCCCGCTCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc        <  2:53470/90‑1 (MQ=255)
                                                                            cctgcctaCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTcc    <  2:796507/90‑1 (MQ=255)
                                                                             ctgcctacCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg   >  2:1156357/1‑90 (MQ=255)
                                                                              tgcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt  <  1:724111/90‑1 (MQ=255)
                                                                              tgcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCACCATACCTGGGAATCAATCGGTCGTCCGt  >  1:739850/1‑90 (MQ=255)
                                                                               gcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                           <  2:629331/64‑1 (MQ=255)
                                                                               gcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                           >  1:629331/1‑64 (MQ=255)
                                                                               gcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                           <  2:1018225/64‑1 (MQ=255)
                                                                               gcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                           >  1:1018225/1‑64 (MQ=255)
                                                                                   |                                                                                    
AAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGT  >  NZ_CP009273/49815‑49982

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