Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,884 C→T P21L (CCG→CTG)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,8840CT97.0% 102.5 / NA 33P21L (CCG→CTG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (12/20);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGA  >  NZ_CP009273/49797‑49965
                                                                                       |                                                                                 
aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGt                                                                                 <  2:226568/90‑1 (MQ=255)
     tttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAc                                                                            <  2:286505/90‑1 (MQ=255)
              cGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCtgcctgcc                                                                   <  1:610963/90‑1 (MQ=255)
                 gAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGAt                                                                <  2:804280/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  1:143869/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  1:662276/90‑1 (MQ=255)
                         aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACCGCATGCCGTGGAACCTGCCAGCCGCGCTCACCTg                                                        >  1:491974/1‑90 (MQ=255)
                            gATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGtt                                                     <  1:579896/90‑1 (MQ=255)
                               cAGTCTGATTGCGGCGTTAACGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                  >  2:470895/1‑90 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  1:287952/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  2:173186/90‑1 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAc                                                                            <  2:511248/52‑1 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAc                                                                            >  1:511248/1‑52 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                      >  1:728695/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                      >  1:935855/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                      >  2:534278/1‑90 (MQ=255)
                                                tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                                 <  2:886665/90‑1 (MQ=255)
                                                 aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                <  1:172747/90‑1 (MQ=255)
                                                     gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                            <  2:797420/90‑1 (MQ=255)
                                                       aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                          <  1:534278/90‑1 (MQ=255)
                                                       aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                          <  1:463886/90‑1 (MQ=255)
                                                            gcgTTATCGGCATGTAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                     >  1:285585/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGGACCTGCCTGCCGATCTCGCCTGGGTTAAACGCAACACCTTAAATAAACCCGGGATTa                     >  1:871215/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          >  1:89247/1‑53 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  2:89247/53‑1 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                     >  1:384924/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAAAACCTTAAATAAACCCCTGATTa                     >  1:856028/1‑90 (MQ=255)
                                                                       aTGTAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt          <  2:285585/90‑1 (MQ=255)
                                                                       aTGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt          <  1:165370/90‑1 (MQ=255)
                                                                                   aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                       >  2:301321/1‑49 (MQ=255)
                                                                                   aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                       <  1:301321/49‑1 (MQ=255)
                                                                                   aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa  <  1:840582/86‑1 (MQ=255)
                                                                                   aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa  >  2:840582/1‑86 (MQ=255)
                                                                                       |                                                                                 
ACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGA  >  NZ_CP009273/49797‑49965

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