Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,905 T→G L28R (CTC→CGC)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9050TG96.0% 79.1 / ‑3.2 25L28R (CTC→CGC) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (10/14);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAC  >  NZ_CP009273/49828‑49991
                                                                             |                                                                                      
cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaa                                                                            <  1:342130/90‑1 (MQ=255)
 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGGTTaaa                                                                           >  2:273265/1‑90 (MQ=255)
   tCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCTATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACg                                                                         >  2:533405/1‑90 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCaa                                                                      <  1:23057/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAac                                                                     <  2:619997/90‑1 (MQ=255)
          gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACAcc                                                                  >  1:418650/1‑90 (MQ=255)
           cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                                 <  2:229092/90‑1 (MQ=255)
           cggcgTTAGCGGTAGACCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                                 <  2:468522/90‑1 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCtt                                                                >  2:632238/1‑90 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCtt                                                                >  2:561883/1‑90 (MQ=255)
             gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTta                                                               <  1:588879/90‑1 (MQ=255)
                  aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaaataa                                                          <  1:416598/90‑1 (MQ=255)
                        aGATCGCGTTATCGGCATGGAAAACGCTATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                                                    <  1:533405/90‑1 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTa                                                                             >  2:38664/1‑60 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTa                                                                             <  1:38664/60‑1 (MQ=255)
                                  aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                          <  1:561883/90‑1 (MQ=255)
                                    cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                                        >  2:47684/1‑90 (MQ=255)
                                       cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                     <  2:418650/90‑1 (MQ=255)
                                        aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                    <  1:228616/90‑1 (MQ=255)
                                                cGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                                 <  2:382086/53‑1 (MQ=255)
                                                cGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                                 >  1:382086/1‑53 (MQ=255)
                                                  ccATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                          >  2:524579/1‑90 (MQ=255)
                                                                   cctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                                                    <  2:187750/47‑1 (MQ=255)
                                                                   cctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                                                    >  1:187750/1‑47 (MQ=255)
                                                                          aTCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc  <  2:591818/90‑1 (MQ=255)
                                                                             |                                                                                      
CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAC  >  NZ_CP009273/49828‑49991

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