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,9030TA100.0% 152.5 / NA 44D27E (GAT→GAAfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (19/25);  total (19/25)

TCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTG  >  NZ_CP009273/49820‑49984
                                                                                   |                                                                                 
tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcc                                                                             <  1:148089/90‑1 (MQ=255)
tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcc                                                                             <  2:78304/90‑1 (MQ=255)
tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGcc                                                                             <  1:401505/90‑1 (MQ=255)
  aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTg                                                                           <  1:591787/90‑1 (MQ=255)
  aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTg                                                                           <  2:749010/90‑1 (MQ=255)
      aTCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGttt                                                                       >  2:183689/1‑90 (MQ=255)
        cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGGCATCCCGTGGAACCTGCCTGCCCCACTCGCCTGGcttca                                                                     >  1:401246/1‑88 (MQ=255)
         aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGGCTGGTTTaaa                                                                    >  2:404249/1‑90 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                               <  2:977523/90‑1 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                               <  1:538735/90‑1 (MQ=255)
              gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                               <  1:676247/90‑1 (MQ=255)
               aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                              <  2:1115172/90‑1 (MQ=255)
                   cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCt                                                          <  2:654577/90‑1 (MQ=255)
                   cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCt                                                          <  1:169645/90‑1 (MQ=255)
                    ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCtt                                                         >  1:574466/1‑90 (MQ=255)
                       gTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaa                                                      >  2:927139/1‑90 (MQ=255)
                            cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                                 <  2:886005/90‑1 (MQ=255)
                               tAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                                              <  1:1146846/90‑1 (MQ=255)
                                 gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                            >  1:50720/1‑90 (MQ=255)
                                 gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                            >  1:207817/1‑90 (MQ=255)
                                     gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                        >  1:922190/1‑90 (MQ=255)
                                        ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                     >  1:535221/1‑90 (MQ=255)
                                        ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                     <  2:740879/90‑1 (MQ=255)
                                         tATCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                    <  1:424042/90‑1 (MQ=255)
                                          aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                   <  2:426066/90‑1 (MQ=255)
                                            cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                                 <  2:1150345/90‑1 (MQ=255)
                                             ggCATGGAAACCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                                <  1:294972/90‑1 (MQ=255)
                                                 tGGAAAACGCCATGCCGTGGGACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                            >  1:1234039/1‑90 (MQ=255)
                                                            aTGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                 >  1:253630/1‑90 (MQ=255)
                                                            aTGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                 >  1:383498/1‑90 (MQ=255)
                                                              gCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat               <  2:535221/90‑1 (MQ=255)
                                                               ccGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc              <  2:135620/90‑1 (MQ=255)
                                                               ccGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc              <  2:377902/90‑1 (MQ=255)
                                                                       cctgcctgCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTcc      <  1:644655/90‑1 (MQ=255)
                                                                       cctgcctgCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTcc      <  2:63973/90‑1 (MQ=255)
                                                                       cctgccggCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTcc      <  1:432572/90‑1 (MQ=255)
                                                                        ctgcctgcCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg     >  1:865701/1‑90 (MQ=255)
                                                                         tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                     >  2:203792/1‑73 (MQ=255)
                                                                         tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa                     <  1:203792/73‑1 (MQ=255)
                                                                         tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt    >  2:181023/1‑90 (MQ=255)
                                                                         tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt    >  2:668185/1‑90 (MQ=255)
                                                                         tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt    >  1:813002/1‑90 (MQ=255)
                                                                         tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt    >  1:730460/1‑90 (MQ=255)
                                                                           cctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTg  >  1:1016535/1‑90 (MQ=255)
                                                                                   |                                                                                 
TCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTG  >  NZ_CP009273/49820‑49984

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