Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,905 T→G L28R (CTC→CGC)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9050TG97.8% 139.5 / NA 46L28R (CTC→CGC) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (22/23);  minor base C (1/0);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCA  >  NZ_CP009273/49822‑49987
                                                                                   |                                                                                  
aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTg                                                                              <  2:659299/90‑1 (MQ=255)
aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTg                                                                              <  2:212807/90‑1 (MQ=255)
aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTAGCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTg                                                                              <  1:190058/90‑1 (MQ=255)
     tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTa                                                                         <  2:102572/90‑1 (MQ=255)
      cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaa                                                                        >  2:158251/1‑90 (MQ=255)
      cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaa                                                                        <  1:29886/90‑1 (MQ=255)
      cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaa                                                                        >  1:520007/1‑90 (MQ=255)
            gATTGCGGCTTTAGCGGTAGATCGCGCTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCaa                                                                  <  2:443095/90‑1 (MQ=255)
            gATTGCGGCGTTAGCGGTAGCTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCaa                                                                  <  1:528365/90‑1 (MQ=255)
            gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCaa                                                                  <  2:218120/90‑1 (MQ=255)
            gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCaa                                                                  <  1:506885/90‑1 (MQ=255)
            gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAACGCCGTGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAac                                                                 <  1:688067/90‑1 (MQ=255)
             aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAac                                                                 >  2:15754/1‑90 (MQ=255)
             aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAac                                                                 >  2:260305/1‑90 (MQ=255)
                 cggctTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                             <  2:385255/90‑1 (MQ=255)
                 cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                             >  2:32319/1‑90 (MQ=255)
                  ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCtt                                                            >  1:402550/1‑90 (MQ=255)
                  ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCtt                                                            >  2:196574/1‑90 (MQ=255)
                        aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaaataa                                                      <  1:245254/90‑1 (MQ=255)
                               gATCGCGTTATCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGATCCCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                               >  1:424867/1‑90 (MQ=255)
                                 tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                                             >  2:631853/1‑90 (MQ=255)
                                 tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGGTTAAACGCAACACCTTAAATAAACCCGTGAt                                             >  2:645896/1‑90 (MQ=255)
                                   gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                                                       >  1:659695/1‑62 (MQ=255)
                                   gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                                                       <  2:659695/62‑1 (MQ=255)
                                   gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                           >  2:95999/1‑90 (MQ=255)
                                   gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGGGAttt                                           >  2:463653/1‑89 (MQ=255)
                                   gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGGGATTa                                           >  2:409707/1‑90 (MQ=255)
                                      ttATCGGTATGGAAAACGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                        <  1:536406/90‑1 (MQ=255)
                                      ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                             <  1:544231/69‑1 (MQ=255)
                                      ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCt                                                             >  2:544231/1‑69 (MQ=255)
                                             cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                 <  1:167073/90‑1 (MQ=255)
                                              aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                <  2:250074/90‑1 (MQ=255)
                                                  aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                            <  2:80967/90‑1 (MQ=255)
                                                  aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                            <  1:177216/90‑1 (MQ=255)
                                                               gTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg               <  1:118719/90‑1 (MQ=255)
                                                               gTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg               >  1:12045/1‑90 (MQ=255)
                                                                 ggAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             <  2:520007/90‑1 (MQ=255)
                                                                 ggAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             <  2:229211/90‑1 (MQ=255)
                                                                       tgcctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             <  1:312171/84‑1 (MQ=255)
                                                                       tgcctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc             >  2:312171/1‑84 (MQ=255)
                                                                       tgcctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt       >  1:349125/1‑90 (MQ=255)
                                                                       tgcctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt       >  2:675113/1‑90 (MQ=255)
                                                                       tgcctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGGCGTCCGt       >  1:99293/1‑90 (MQ=255)
                                                                         cctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTg     >  1:458188/1‑90 (MQ=255)
                                                                         cctgccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTg     >  1:532810/1‑90 (MQ=255)
                                                                            gccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCa  >  1:224616/1‑90 (MQ=255)
                                                                                   |                                                                                  
AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCA  >  NZ_CP009273/49822‑49987

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