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,8840CT96.3% 74.2 / ‑4.9 27P21L (CCG→CTG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (9/17);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

AATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGA  >  NZ_CP009273/49799‑49965
                                                                                     |                                                                                 
aaTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTgg                                                                               >  2:616589/1‑90 (MQ=255)
                       aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTg                                                        <  2:661349/90‑1 (MQ=255)
                       aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTg                                                        <  2:453974/90‑1 (MQ=255)
                            tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                                   <  1:616589/90‑1 (MQ=255)
                             cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTg                                                        <  1:37148/84‑1 (MQ=255)
                             cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTg                                                        >  2:37148/1‑84 (MQ=255)
                             cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                  >  2:77749/1‑90 (MQ=255)
                             cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                  >  2:691668/1‑90 (MQ=255)
                                  tGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                             <  2:211272/90‑1 (MQ=255)
                                   gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  1:85879/90‑1 (MQ=255)
                                   gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  1:50514/90‑1 (MQ=255)
                                    aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           <  1:481540/90‑1 (MQ=255)
                                         ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGGTTAAACGCAACACCtt                                      >  1:575694/1‑90 (MQ=255)
                                         ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          >  1:110854/1‑70 (MQ=255)
                                         ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  2:110854/70‑1 (MQ=255)
                                              tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGACCTCGCCTGGTTTAAACACAACACCTtaaata                                 <  2:545593/90‑1 (MQ=255)
                                               aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                <  2:659170/90‑1 (MQ=255)
                                                     aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                          <  2:610622/90‑1 (MQ=255)
                                                     aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                          <  2:637165/90‑1 (MQ=255)
                                                          gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                     >  1:556103/1‑90 (MQ=255)
                                                          gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                     >  2:693356/1‑90 (MQ=255)
                                                             ttATCGGCATGGGAAACGCCATGCTGTGGAACCTGCCCGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                  >  2:627667/1‑90 (MQ=255)
                                                                     aTGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt          <  2:677919/90‑1 (MQ=255)
                                                                      gggAAACCGCCATGCTGTGGAACCTCCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa         <  1:147373/89‑1 (MQ=255)
                                                                         aaaaCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt      <  2:556103/90‑1 (MQ=255)
                                                                             cGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa  >  1:381651/1‑90 (MQ=255)
                                                                             cGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa  <  1:77749/90‑1 (MQ=255)
                                                                                     |                                                                                 
AATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGA  >  NZ_CP009273/49799‑49965

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