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,8840CT98.5% 208.3 / ‑4.2 65P21L (CCG→CTG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/44);  total (20/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

ACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGA  >  NZ_CP009273/49797‑49965
                                                                                       |                                                                                 
aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGt                                                                                 <  2:1191609/90‑1 (MQ=255)
 cAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCtgtg                                                                                >  1:621702/1‑90 (MQ=255)
   aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGa                                                                              <  2:165868/90‑1 (MQ=255)
        ttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAcctg                                                                         <  1:1406521/90‑1 (MQ=255)
         tttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAcctgc                                                                        <  2:856687/90‑1 (MQ=255)
             tCGGGAAATCTCAATGATCAGTCTGATTGCGGGGTTAGCGGTAGATCGCGTTACCGGCATGGAAAACGCCATGCTGTGGAACctgcctgc                                                                    <  2:1843208/90‑1 (MQ=255)
             tCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACctgcctgc                                                                    <  2:90174/90‑1 (MQ=255)
             tCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACctgcctgc                                                                    >  1:1573732/1‑90 (MQ=255)
                  aaaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGAtc                                                               <  1:1480076/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTTCCTGCCGATCTCGcc                                                          <  2:1112985/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  2:137684/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  2:1886027/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  1:527183/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  2:214528/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          <  2:24659/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAGCTCGcc                                                          <  1:149149/90‑1 (MQ=255)
                        cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCt                                                         <  2:974673/90‑1 (MQ=255)
                               cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGGTTaa                                                  >  2:1720761/1‑90 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  2:1573272/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  1:1514071/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            >  1:1179417/1‑90 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  1:654081/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCCCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  2:975908/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGACCGCGTTTTCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                            <  1:913676/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           <  1:651159/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           <  1:892462/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           <  2:590148/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           <  1:1175125/90‑1 (MQ=255)
                                         gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCTCGcc                                                          <  1:674376/72‑1 (MQ=255)
                                         gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                          >  2:674376/1‑72 (MQ=255)
                                          cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                       <  2:1242451/90‑1 (MQ=255)
                                          cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                       <  2:455510/90‑1 (MQ=255)
                                          cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                       <  1:1116824/90‑1 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGt                                                                                 <  2:786999/47‑1 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGt                                                                                 >  1:786999/1‑47 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                      >  1:1886183/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                      >  2:677508/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGGGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                      >  1:1177443/1‑90 (MQ=255)
                                                tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                                 <  1:176581/90‑1 (MQ=255)
                                                tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                                 >  2:563800/1‑90 (MQ=255)
                                                tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                                 >  2:603814/1‑90 (MQ=255)
                                                   cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                              <  1:1149997/90‑1 (MQ=255)
                                                     gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                            <  2:691274/90‑1 (MQ=255)
                                                       ataTCGCTTTGCCGGCATGGAAAACCCAATGCTGTGGAACCCGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                          <  1:1723400/88‑1 (MQ=255)
                                                       aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                          <  2:248153/90‑1 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           <  1:1604939/72‑1 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           >  2:1604939/1‑72 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCggg                         >  2:1789213/1‑88 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                         >  1:1106278/1‑90 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                         >  2:885737/1‑90 (MQ=255)
                                                          tCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                       >  1:1577421/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           >  1:1513738/1‑68 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                           <  2:1513738/68‑1 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                     >  2:465366/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                     >  2:1671474/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAAGCAACACCTTAAATAAACCCGTGATTa                     >  1:1718437/1‑90 (MQ=255)
                                                               ttATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                  <  2:1577421/90‑1 (MQ=255)
                                                                    ggCATGAAAACCCCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc             <  1:34498/90‑1 (MQ=255)
                                                                      cATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa           <  2:1886183/90‑1 (MQ=255)
                                                                       aTGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt          <  2:143663/90‑1 (MQ=255)
                                                                       aTGAAAACCGCCATGTTGTGAAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt          <  2:1046009/90‑1 (MQ=255)
                                                                       aTGAAAACCCCCATGCTGGGAAACCTCCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt          <  1:837076/90‑1 (MQ=255)
                                                                           aaaaCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt      <  2:24572/90‑1 (MQ=255)
                                                                               cGCCATGCTGGGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa  <  2:284748/90‑1 (MQ=255)
                                                                                   aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa           >  2:633097/1‑77 (MQ=255)
                                                                                   aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa           <  1:633097/77‑1 (MQ=255)
                                                                                       |                                                                                 
ACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGA  >  NZ_CP009273/49797‑49965

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