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.6% 132.7 / ‑3.2 41D27E (GAT→GAAfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (21/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

CAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCA  >  NZ_CP009273/49821‑49987
                                                                                  |                                                                                    
cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCt                                                                               >  2:1102276/1‑90 (MQ=255)
cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCt                                                                               <  2:1001311/90‑1 (MQ=255)
 aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTg                                                                              <  2:1504293/90‑1 (MQ=255)
       cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGGAACTTGCTGCCGAACTCGCGTCGTTTaa                                                                        >  1:875305/1‑90 (MQ=255)
       cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTaa                                                                        >  2:1299275/1‑90 (MQ=255)
             gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                                  <  2:1513620/90‑1 (MQ=255)
             gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                                  <  1:89227/90‑1 (MQ=255)
              aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                                 <  1:1147036/90‑1 (MQ=255)
              aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                                 <  2:484936/90‑1 (MQ=255)
              aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                                 >  1:631255/1‑90 (MQ=255)
                   ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCtt                                                            >  2:1231045/1‑90 (MQ=255)
                         aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaataa                                                      <  1:185663/90‑1 (MQ=255)
                           cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                                    >  1:42168/1‑90 (MQ=255)
                             gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                                  <  2:42168/90‑1 (MQ=255)
                               aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                                                <  2:14712/90‑1 (MQ=255)
                                gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                               >  2:725124/1‑90 (MQ=255)
                                gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                               >  2:533187/1‑90 (MQ=255)
                                gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                               >  2:1073342/1‑90 (MQ=255)
                                    gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAAAGCAACACCTTAAATAAACCCGTGATTa                                           >  1:1095937/1‑90 (MQ=255)
                                     cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGGTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                          >  1:770898/1‑90 (MQ=255)
                                        tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                       <  2:1127123/90‑1 (MQ=255)
                                        tATCGGCATGAAAAACGCCATGCCGTGGACCCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                       <  1:1576083/90‑1 (MQ=255)
                                         aTCGGCATGCAAACCGCCTCGCCGTGGACCCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                      <  1:1518924/90‑1 (MQ=255)
                                          tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                                  >  1:354929/1‑61 (MQ=255)
                                          tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                                  <  2:354929/61‑1 (MQ=255)
                                               aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCAt                                <  1:725537/90‑1 (MQ=255)
                                                   aaaaCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                            <  1:725124/90‑1 (MQ=255)
                                                    aaaCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                           <  1:885057/90‑1 (MQ=255)
                                                       cGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                        >  1:522068/1‑90 (MQ=255)
                                                       cGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                        >  1:281402/1‑90 (MQ=255)
                                                       cGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                        <  2:349844/90‑1 (MQ=255)
                                                          cATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                     >  2:621744/1‑90 (MQ=255)
                                                           aTGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                    >  2:1472477/1‑90 (MQ=255)
                                                                gTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                                  >  1:281382/1‑39 (MQ=38)
                                                                gTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                                  <  2:281382/39‑1 (MQ=38)
                                                                  ggAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             <  2:37386/90‑1 (MQ=255)
                                                                      cctgcctgCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTcc         <  2:1095937/90‑1 (MQ=255)
                                                                        tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt       >  1:1163989/1‑90 (MQ=255)
                                                                        tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt       >  1:1244247/1‑90 (MQ=255)
                                                                        tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt       >  1:364331/1‑90 (MQ=255)
                                                                             gccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCa  >  2:237902/1‑90 (MQ=255)
                                                                                  |                                                                                    
CAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCA  >  NZ_CP009273/49821‑49987

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