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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9030TG96.2% 81.4 / ‑3.6 26D27E (GAT→GAGfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (13/12);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAC  >  NZ_CP009273/49828‑49991
                                                                           |                                                                                        
cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTaa                                                                            >  1:75911/1‑90 (MQ=255)
 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTaaa                                                                           >  2:746729/1‑90 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCaa                                                                      <  1:440914/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAac                                                                     <  1:724147/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAac                                                                     <  2:657602/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAac                                                                     <  2:104441/90‑1 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCtt                                                                >  2:769689/1‑90 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGGTTAAACGCAACACCtt                                                                >  1:364530/1‑90 (MQ=255)
                  aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTtaaataa                                                          <  2:40819/90‑1 (MQ=255)
                      gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                                      <  1:126827/90‑1 (MQ=255)
                         gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCggg                                                   >  1:445051/1‑88 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTtaaataaa                                                         <  1:258914/80‑1 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTtaaataaa                                                         >  2:258914/1‑80 (MQ=255)
                                ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAGCACCTTAAATAAACCCGTGATTATgg                                            >  2:531216/1‑90 (MQ=255)
                                   tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                                         <  2:178190/90‑1 (MQ=255)
                                       cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                     <  1:769689/90‑1 (MQ=255)
                                        aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                    <  2:75911/90‑1 (MQ=255)
                                        aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                    >  1:378452/1‑90 (MQ=255)
                                                cGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                    <  1:390512/82‑1 (MQ=255)
                                                cGCCATGCCGTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                    >  2:390512/1‑82 (MQ=255)
                                                         gTGGAACCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                   >  2:786310/1‑90 (MQ=255)
                                                              aCCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCt                                                                 <  2:110968/39‑1 (MQ=39)
                                                              aCCTGCCTGCCGAGCTCGCCTGGTTTAAACGCAACACCt                                                                 >  1:110968/1‑39 (MQ=39)
                                                                 tgcctgccGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt           >  1:77846/1‑90 (MQ=255)
                                                                   cctgccGAGCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTg         >  2:218235/1‑90 (MQ=255)
                                                                          aTCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc  <  2:247098/90‑1 (MQ=255)
                                                                           |                                                                                        
CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAC  >  NZ_CP009273/49828‑49991

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