Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,903 T→A D27E (GAT→GAA folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9030TA96.7% 101.3 / NA 30D27E (GAT→GAAfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (10/19);  minor base G (0/1);  total (10/20)
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

TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGT  >  NZ_CP009273/49827‑49982
                                                                            |                                                                               
tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTa                                                                    <  2:438351/90‑1 (MQ=255)
tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTa                                                                    <  2:648085/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGACCCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                             <  2:336914/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                             <  2:831024/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                             <  1:458512/90‑1 (MQ=255)
        aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                            <  2:423292/90‑1 (MQ=255)
        aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                            >  2:727259/1‑90 (MQ=255)
            cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCt                                                        <  2:486661/90‑1 (MQ=255)
                gTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaa                                                    >  2:697317/1‑90 (MQ=255)
                 ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaat                                                   <  1:638779/90‑1 (MQ=255)
                    gCGGTAGACCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaataaa                                                <  1:467877/90‑1 (MQ=255)
                      ggTAGTTCGCGTTATCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAAcc                                              <  2:132701/90‑1 (MQ=255)
                                 ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                             >  2:173994/1‑64 (MQ=255)
                                 ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                             <  1:173994/64‑1 (MQ=255)
                                 ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                   <  1:360416/90‑1 (MQ=255)
                                  tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                  <  1:690732/90‑1 (MQ=255)
                                  tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                  <  1:664487/90‑1 (MQ=255)
                                      ggCATGGAAAACGCCATGCCGTGGACCCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                              <  2:311118/90‑1 (MQ=255)
                                        cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                            <  2:579763/90‑1 (MQ=255)
                                         aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAGACGCaa                                                             >  2:901996/1‑56 (MQ=255)
                                         aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAGACGCaa                                                             <  1:901996/56‑1 (MQ=255)
                                          tGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                          <  1:827533/90‑1 (MQ=255)
                                                aCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCt                                                        >  2:445019/1‑54 (MQ=255)
                                                aCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCt                                                        <  1:445019/54‑1 (MQ=255)
                                                      tGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa              >  2:509875/1‑90 (MQ=255)
                                                      tGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa              >  1:242406/1‑90 (MQ=255)
                                                            ggAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc        >  1:887435/1‑90 (MQ=255)
                                                            ggAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc        <  1:727259/90‑1 (MQ=255)
                                                                 ctgcctgcCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg   >  2:541236/1‑90 (MQ=255)
                                                                  tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt  >  2:887656/1‑90 (MQ=255)
                                                                            |                                                                               
TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGT  >  NZ_CP009273/49827‑49982

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