Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,912 G→T W30C (TGG→TGT folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9120GT96.2% 172.4 / 0.8 53W30C (TGG→TGTfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (24/27);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49827‑49996
                                                                                     |                                                                                    
tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTa                                                                                  <  1:858736/90‑1 (MQ=255)
 cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaa                                                                                 >  1:198029/1‑90 (MQ=255)
 cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaa                                                                                 >  1:54211/1‑90 (MQ=255)
  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa                                                                                >  2:857200/1‑90 (MQ=255)
  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa                                                                                >  1:463814/1‑90 (MQ=255)
  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa                                                                                >  2:1227284/1‑90 (MQ=255)
      tGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  2:1434178/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  2:1420508/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  1:1429230/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  1:1357064/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCACGTTATCGGGATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  1:330972/90‑1 (MQ=255)
       ctttGCGCCGTTAGCGGTAGTTCGCGTTATCGGCATGGAAAACGCCATGCCGTGAACCCTTCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  1:1462495/88‑1 (MQ=255)
           gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACAcc                                                                       <  2:759880/90‑1 (MQ=255)
            cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  1:444801/90‑1 (MQ=255)
            cggcgTTAGCGGTAGATCGCGTTATCGGCATGAAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  1:1263890/90‑1 (MQ=255)
            cggcgTTAGCGGTAGATCGCGTTATCGGCATGAAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  1:1295166/90‑1 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           >  1:779785/1‑84 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  2:779785/84‑1 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  2:421864/1‑90 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  1:444231/1‑90 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  2:116676/89‑1 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      >  1:116676/1‑89 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAATCAACACCtt                                                                     >  1:145756/1‑90 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCCTGCCGGGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  1:1434615/1‑90 (MQ=255)
                  tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTtaaata                                                                <  2:96700/90‑1 (MQ=255)
                            tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGttt                                                                                    <  2:1010384/60‑1 (MQ=255)
                            tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGttt                                                                                    >  1:1010384/1‑60 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGGTCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGGGGTTa                                                    >  1:918094/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:1418400/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:488525/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:523881/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGGGATTa                                                    >  1:123720/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAAGCAACACCTTAAATAAACCCGTGGTTa                                                    >  1:920751/1‑90 (MQ=255)
                                 ttATCGGCATGGAAAACGCCTTGCCGTGGACCCTCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 <  1:857200/90‑1 (MQ=255)
                                 ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 <  2:198029/90‑1 (MQ=255)
                                  tATCGGCATGGAAATCGCCACGCCGGGGAGCCTGCCTGCGGATCTAGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                <  2:826904/90‑1 (MQ=255)
                                  tATCGGCATGAAAACCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                <  1:86761/90‑1 (MQ=255)
                                      ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                                            <  2:1419037/90‑1 (MQ=255)
                                      ggCATGAAAACCCCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                                            <  1:187617/90‑1 (MQ=255)
                                             aaacCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                     <  1:322438/90‑1 (MQ=255)
                                             aaaaCGCCATGCGGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                     <  1:1241477/90‑1 (MQ=255)
                                             aaaaCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                     <  2:321853/90‑1 (MQ=255)
                                             aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                     >  1:1301868/1‑90 (MQ=255)
                                                 cGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCTATACCTGGGa                                 <  2:35766/90‑1 (MQ=255)
                                                   ccATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                               <  1:296585/90‑1 (MQ=255)
                                                     aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                             >  1:1251501/1‑90 (MQ=255)
                                                     aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                             >  1:1454064/1‑90 (MQ=255)
                                                        ccGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc                          <  2:182631/90‑1 (MQ=255)
                                                          gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                        >  2:402859/1‑90 (MQ=255)
                                                          gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                        >  2:529518/1‑90 (MQ=255)
                                                           tGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt                       >  2:1448597/1‑90 (MQ=255)
                                                            ggAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  2:1301868/90‑1 (MQ=255)
                                                                 ctgcctgcCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATGGGTCGTCCg                 >  2:1313117/1‑90 (MQ=255)
                                                                  tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAACCGGTCGGCCGt                >  1:1444168/1‑90 (MQ=255)
                                                                                gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  2:1418400/90‑1 (MQ=255)
                                                                                     |                                                                                    
TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49827‑49996

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