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,9030TA97.5% 133.8 / ‑5.4 40D27E (GAT→GAAfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (18/21);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCA  >  NZ_CP009273/49820‑49987
                                                                                   |                                                                                    
tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcc                                                                                <  2:252530/90‑1 (MQ=255)
 cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcca                                                                               <  1:456711/90‑2 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                                        >  1:218705/1‑90 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTaa                                                                        >  2:795128/1‑90 (MQ=255)
         aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcc                                                                                <  2:436021/81‑1 (MQ=255)
         aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcc                                                                                >  1:436021/1‑81 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                  <  1:827998/90‑1 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                  <  2:47646/90‑1 (MQ=255)
               aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                 <  1:238508/90‑1 (MQ=255)
               aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                 >  1:92375/1‑90 (MQ=255)
               aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                 <  1:318354/90‑1 (MQ=255)
               aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                 <  1:667745/90‑1 (MQ=255)
               aTTGCGGCGTTAGCGGTAGATCGCGTTATCGCCATGGAAAACACCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                 <  1:572397/90‑1 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcca                                                                               >  1:270691/1‑70 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcca                                                                               <  2:270691/71‑2 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCtt                                                            >  2:704058/1‑90 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCtt                                                            >  2:704592/1‑90 (MQ=255)
                         tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTtaaata                                                       <  2:807634/90‑1 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCCTGCCGGGGAACCTGCCTGCCCAACTCGCCAGGGTTAAAAGCAAAACCTTAAATAAACCCGGGATTa                                           >  2:702727/1‑90 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGGACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                           >  2:872875/1‑90 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                  <  2:123371/67‑1 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                  >  1:123371/1‑67 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGGTTAAACGCAACACCTTAAATAAACCCGTGATTa                                           >  1:614736/1‑90 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCGGCCTGCCGAACTCGCCCGGGTTAAACGCAAAACCTTAAATAAACCCCTGATTa                                           >  1:862355/1‑90 (MQ=255)
                                      cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                          >  1:179450/1‑90 (MQ=255)
                                        ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                  >  2:5491/1‑64 (MQ=255)
                                        ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                  <  1:5491/64‑1 (MQ=255)
                                          aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                      <  2:92375/90‑1 (MQ=255)
                                               cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                 <  2:460768/90‑1 (MQ=255)
                                                aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                <  1:702727/90‑1 (MQ=255)
                                                aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                <  1:587290/90‑1 (MQ=255)
                                                aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                >  1:479641/1‑90 (MQ=255)
                                                 tGGAAAACGCCATGCCGTGGACCCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                               <  1:639210/90‑1 (MQ=255)
                                                           cATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                     >  2:553807/1‑90 (MQ=255)
                                                                 gTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg               <  1:872875/90‑1 (MQ=255)
                                                                   ggAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             <  1:245875/90‑1 (MQ=255)
                                                                         tgcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             <  2:346723/84‑1 (MQ=255)
                                                                         tgcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             >  1:346723/1‑84 (MQ=255)
                                                                          gcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt      >  1:654516/1‑90 (MQ=255)
                                                                              gccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCa  >  2:183819/1‑90 (MQ=255)
                                                                                   |                                                                                    
TCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCA  >  NZ_CP009273/49820‑49987

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