Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,845,796 G→A 7.5% W29* (TGG→TGA)  pfeS‑II → histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,845,7960GA 7.5% 152.0 / 2.7 53W29* (TGG→TGA) pfeS‑IIhistidine kinase
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base A (2/2);  total (24/29)
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.89e-01

AGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAGCGC  >  NC_002947/1845651‑1845940
                                                                                                                                                 |                                                                                                                                                
aGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGc                                                                                                                                               <  1:24241/149‑1 (MQ=255)
aGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGc                                                                                                                                               >  2:24241/1‑149 (MQ=255)
aGATCCGCACCGTGGGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGGGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGc                                                                                                                                               <  7:85118/149‑1 (MQ=255)
         ccGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                                      <  2:260187/149‑1 (MQ=255)
          cGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                                      >  1:38203/1‑148 (MQ=255)
              tgGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAAc                                                                                                                                 <  2:38203/149‑1 (MQ=255)
                 ggCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCa                                                                                                                              <  5:142740/149‑1 (MQ=255)
                   cAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAg                                                                                                                             >  8:40837/1‑148 (MQ=255)
                       ggTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGc                                                                                                                        <  8:41321/149‑1 (MQ=255)
                        gTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGcc                                                                                                                       <  6:151738/149‑1 (MQ=255)
                            cGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTc                                                                                                                    <  7:40837/148‑1 (MQ=255)
                            cGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTc                                                                                                                    <  1:43275/148‑1 (MQ=255)
                                   aGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGa                                                                                                            <  5:330845/149‑1 (MQ=255)
                                    gTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACAt                                                                                                                                       >  2:134752/1‑121 (MQ=255)
                                    gTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACAt                                                                                                                                       <  1:134752/121‑1 (MQ=255)
                                         aGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCgg                                                                                                       >  3:30923/1‑148 (MQ=255)
                                                  gAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGAc                                                                                             >  7:152799/1‑149 (MQ=255)
                                                       aaCCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGc                                                                                        >  3:127128/1‑149 (MQ=255)
                                                          cATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGCTCGGCCTCAGGTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCggg                                                                                     <  8:306839/149‑1 (MQ=255)
                                                                      ccATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccgcc                                                                          >  2:204573/1‑148 (MQ=255)
                                                                       cATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACGCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGa                                                                        <  8:273961/149‑1 (MQ=255)
                                                                         ttCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAg                                                                      >  5:336459/1‑149 (MQ=255)
                                                                                ttCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACa                                                               <  7:207384/149‑1 (MQ=255)
                                                                                  cTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAgg                                                             >  5:313522/1‑149 (MQ=255)
                                                                                   tGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGc                                                            >  2:208484/1‑149 (MQ=255)
                                                                                     gAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGc                                                            <  3:96964/147‑1 (MQ=255)
                                                                                     gAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGc                                                            <  6:336459/147‑1 (MQ=255)
                                                                                       aGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTgg                                                        <  6:313522/149‑1 (MQ=255)
                                                                                          tGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGgcg                                                      >  1:124363/1‑148 (MQ=255)
                                                                                            gCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAg                                                   >  7:161049/1‑149 (MQ=255)
                                                                                            gCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAg                                                   <  8:152799/149‑1 (MQ=255)
                                                                                            gCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAg                                                   >  1:125806/1‑149 (MQ=255)
                                                                                             ccATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGt                                                  <  1:208484/149‑1 (MQ=255)
                                                                                                     cTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCt                                                          <  4:285312/133‑1 (MQ=255)
                                                                                                     cTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCt                                                          >  3:285312/1‑133 (MQ=255)
                                                                                                                   cctgcTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGAc                             <  4:30923/148‑1 (MQ=255)
                                                                                                                       ctgATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAATCTGAACAGGCCTGGCGCAGTGGCGGGCAGGCAGGGCTGGACCAgt                         <  2:124363/148‑1 (MQ=255)
                                                                                                                        tgatgaTCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTAt                                                                                >  7:147851/1‑92 (MQ=255)
                                                                                                                        tgatgaTCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTAt                                                                                <  8:147851/92‑1 (MQ=255)
                                                                                                                            tgccGGCCCCAGGTATAACTGAGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGTAGGGCTGGACCAgtgggtgg                   >  7:111610/4‑149 (MQ=255)
                                                                                                                                 gCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGcc              <  4:127128/149‑1 (MQ=255)
                                                                                                                                   cTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCAt            <  7:45341/149‑1 (MQ=255)
                                                                                                                                    tCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTCGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATg           >  6:245970/1‑149 (MQ=255)
                                                                                                                                      aGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTGGGAGCCGGCACGCCAGAGCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGTGCTGGACCAGTGGGTGGCGGCCTTGc          <  5:245970/148‑1 (MQ=255)
                                                                                                                                        ttATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggc                                               <  3:74441/109‑1 (MQ=255)
                                                                                                                                        ttATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggc                                               >  4:74441/1‑109 (MQ=255)
                                                                                                                                          ataCCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAg     >  4:119906/1‑149 (MQ=255)
                                                                                                                                             ccTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgt                         >  7:19135/1‑126 (MQ=255)
                                                                                                                                             ccTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgt                         <  8:19135/126‑1 (MQ=255)
                                                                                                                                             ccTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTggcgg                <  1:167662/135‑1 (MQ=255)
                                                                                                                                             ccTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTggcgg                >  2:167662/1‑135 (MQ=255)
                                                                                                                                              cTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAgcgc  >  3:17956/1‑148 (MQ=255)
                                                                                                                                              cTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAgcgc  >  7:290876/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAGCGC  >  NC_002947/1845651‑1845940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑