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,9030TA91.7% 28.4 / ‑3.8 12D27E (GAT→GAAfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (4/7);  total (5/7)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

GATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCG  >  NZ_CP009273/49834‑49981
                                                                     |                                                                              
gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                            <  2:20916/90‑1 (MQ=255)
gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                            <  2:303997/90‑1 (MQ=255)
gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                            <  2:70665/90‑1 (MQ=255)
 aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                           <  1:213876/90‑1 (MQ=255)
      ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCtt                                                      >  1:327023/1‑90 (MQ=255)
           tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaata                                                 >  2:147593/1‑90 (MQ=255)
                                 cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                           <  2:327023/90‑1 (MQ=255)
                                  aTGGAAAAGGCCATCCCGGGGAACGTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                          <  1:115817/90‑1 (MQ=255)
                                          cGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                  >  2:233483/1‑90 (MQ=255)
                                                          ctgcctgcCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg  >  1:301964/1‑90 (MQ=255)
                                                          ctgcctgcCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                <  1:333291/60‑1 (MQ=255)
                                                          ctgcctgcCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                >  2:333291/1‑60 (MQ=255)
                                                                     |                                                                              
GATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCG  >  NZ_CP009273/49834‑49981

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