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,8980GA94.7% 57.2 / ‑4.0 19A26T (GCC→ACC) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (8/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

AAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCG  >  NZ_CP009273/49815‑49981
                                                                                   |                                                                                   
aaaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGAtc                                                                               <  1:317691/90‑1 (MQ=255)
     tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGcc                                                                          <  1:617863/90‑1 (MQ=255)
             cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTaa                                                                  >  1:539490/1‑90 (MQ=255)
             cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTaa                                                                  >  2:103077/1‑90 (MQ=255)
                   gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCaa                                                            <  1:476962/90‑1 (MQ=255)
                    aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAac                                                           <  1:590095/90‑1 (MQ=255)
                    aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAac                                                           <  2:197632/90‑1 (MQ=255)
                        cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCt                                                       <  2:24709/90‑1 (MQ=255)
                        cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCt                                                       <  1:613343/90‑1 (MQ=255)
                         ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCtt                                                      >  1:507372/1‑90 (MQ=255)
                              tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaata                                                 <  1:244321/90‑1 (MQ=255)
                                          gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGAtc                                                                               >  1:224819/1‑48 (MQ=255)
                                          gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGAtc                                                                               <  2:224819/48‑1 (MQ=255)
                                          gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                     >  2:154307/1‑90 (MQ=255)
                                             ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                  >  1:137352/1‑90 (MQ=255)
                                             ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                  >  2:263024/1‑90 (MQ=255)
                                                    cATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                           <  1:118356/90‑1 (MQ=255)
                                                          aaaCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                     <  2:539490/90‑1 (MQ=255)
                                                                        ggAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCAgtc       >  1:165588/1‑90 (MQ=255)
                                                                             ctgcctacCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg  >  2:189970/1‑90 (MQ=255)
                                                                                   |                                                                                   
AAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCG  >  NZ_CP009273/49815‑49981

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