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.7% 251.5 / ‑3.9 79P21L (CCG→CTG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/50);  total (28/51)
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.04e-01

ACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAA  >  NZ_CP009273/49797‑49970
                                                                                       |                                                                                      
aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGt                                                                                      <  1:1015564/90‑1 (MQ=255)
                  aaaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGAtc                                                                    <  2:1462974/90‑1 (MQ=255)
                  aaaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGAtc                                                                    <  1:1403996/90‑1 (MQ=255)
                  aaaTCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGAtc                                                                    <  1:1492599/90‑1 (MQ=255)
                  aaaTCTAAATGATCAGTCTGATTGCGGCGTTAGCGGTAGACCGCGTTATCGGCATGGAAACCGCCATGCTGTGGACCCTGCCTGCCGAtc                                                                    <  1:1427255/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                               <  2:28754/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                               <  2:1278665/90‑1 (MQ=255)
                       tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGcc                                                               <  2:1239602/90‑1 (MQ=255)
                        cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCt                                                              <  1:345363/90‑1 (MQ=255)
                        cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCt                                                              <  1:1203294/90‑1 (MQ=255)
                        cAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCt                                                              <  2:868448/90‑1 (MQ=255)
                         aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTg                                                             <  1:710118/90‑1 (MQ=255)
                               cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                       >  2:557692/1‑90 (MQ=255)
                               cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                       >  2:733098/1‑90 (MQ=255)
                                aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                      >  2:93382/1‑90 (MQ=255)
                                aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                      >  1:804582/1‑90 (MQ=255)
                                aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaaa                                                      >  2:523362/1‑90 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCTGATCTCGCCTGGTTTAAACGCaa                                                 <  2:729883/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                 <  2:1084355/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                 <  1:733098/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                 <  1:523362/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                 <  2:681764/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                                 <  1:1426924/90‑1 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                                  >  2:569591/1‑89 (MQ=255)
                                     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                                  <  1:569591/89‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                <  2:551746/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                <  1:1324363/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                <  2:917063/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                <  2:1157480/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                <  1:808976/90‑1 (MQ=255)
                                      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                                <  1:1426852/90‑1 (MQ=255)
                                          cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                            <  2:729180/90‑1 (MQ=255)
                                          cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                            <  1:1616042/90‑1 (MQ=255)
                                          cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                                            <  1:1336084/90‑1 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                           >  2:367157/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                           >  2:285687/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                           >  2:87085/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                           >  1:761639/1‑90 (MQ=255)
                                           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                           >  1:1414854/1‑90 (MQ=255)
                                               ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                                       >  2:88965/1‑90 (MQ=255)
                                               ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGGGGGACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                                       >  1:707912/1‑90 (MQ=255)
                                                   cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAc                                   <  1:1519586/90‑1 (MQ=255)
                                                     gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAccc                                 <  2:1449602/90‑1 (MQ=255)
                                                      tAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                                <  2:733252/90‑1 (MQ=255)
                                                      tAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAAGCAACACCTTAAATAAACCCg                                >  1:1115075/1‑90 (MQ=255)
                                                       aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                               <  1:391785/90‑1 (MQ=255)
                                                       aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt                               <  2:1308177/90‑1 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGAtc                                                                    <  2:378924/52‑1 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGAtc                                                                    >  1:378924/1‑52 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                              >  2:1548542/1‑90 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                              >  1:258335/1‑90 (MQ=255)
                                                        gATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAAACCGTg                              >  2:1111092/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                          >  2:1156243/1‑90 (MQ=255)
                                                            gcgTTATCGGCATGGAAAACGCCATGCTGGGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                          >  1:1325955/1‑90 (MQ=255)
                                                               ttATCGGCATGGAAAACGCCATGCTGTGGAACctgcctgc                                                                         <  2:910202/40‑1 (MQ=255)
                                                               ttATCGGCATGGAAAACGCCATGCTGTGGAACctgcctgc                                                                         >  1:910202/1‑40 (MQ=255)
                                                               ttATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                       <  2:1325955/90‑1 (MQ=255)
                                                               ttATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                       <  2:6058/90‑1 (MQ=255)
                                                               ttATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                       >  2:658076/1‑90 (MQ=255)
                                                               ttATCGGCAGGAAAAACGCCATGCTGTGAACCCTCCCTGCCGATCCCGCCTGGTTTAAACGCAACACCTTAAATAAGCCCGTGATTATgg                       <  1:1111092/90‑1 (MQ=255)
                                                                tATCGGCATGGAAACCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                      <  1:667123/90‑1 (MQ=255)
                                                                 aTCGGCATGGAAAACGCCATGCTGTGGAACCtgcctgcc                                                                        >  2:1505170/1‑39 (MQ=255)
                                                                 aTCGGCATGGAAAACGCCATGCTGTGGAACCtgcctgcc                                                                        <  1:1505170/39‑1 (MQ=255)
                                                                 aTCGGAATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                     <  2:1420828/90‑1 (MQ=255)
                                                                   cGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                   <  2:1115075/90‑1 (MQ=255)
                                                                     gCATGGAAACCCCAATGCTGGGGACCCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                 <  1:147581/90‑1 (MQ=255)
                                                                      cATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCt                                                              <  2:918852/44‑1 (MQ=255)
                                                                      cATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCt                                                              >  1:918852/1‑44 (MQ=255)
                                                                      cATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                <  1:367157/90‑1 (MQ=255)
                                                                           aaaaCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt           <  2:1414854/90‑1 (MQ=255)
                                                                           aaaaCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt           <  2:1239151/90‑1 (MQ=255)
                                                                            aaaCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg          <  1:566679/90‑1 (MQ=255)
                                                                            aaaCGCCATGCTGGGGACCCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg          <  1:354384/90‑1 (MQ=255)
                                                                              aCGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg        >  1:850794/1‑90 (MQ=255)
                                                                               cGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                     <  2:1065889/60‑1 (MQ=255)
                                                                               cGCCATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataa                                     >  1:1065889/1‑60 (MQ=255)
                                                                                 ccATGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat     <  1:1156243/90‑1 (MQ=255)
                                                                                   aTGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca   >  1:1356066/1‑90 (MQ=255)
                                                                                    tGCTGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa  >  1:738595/1‑90 (MQ=255)
                                                                                       |                                                                                      
ACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAA  >  NZ_CP009273/49797‑49970

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