Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,846,345:1 +A 11.4% coding (594/1377 nt) pfeS‑II → histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,846,3451.A11.4% 122.8 / 12.1 35coding (594/1377 nt)pfeS‑IIhistidine kinase
Reads supporting (aligned to +/- strand):  ref base . (14/17);  new base A (2/2);  total (16/19)
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.97e-01

GCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGT  >  NC_002947/1846197‑1846475
                                                                                                                                                     |                                                                                                                                  
gCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGTCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGt                                                                                                                                   <  1:36361/150‑1 (MQ=255)
   aGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCgg                                                                                                                                <  1:329880/150‑1 (MQ=255)
   aGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCgg                                                                                                                                >  1:335260/1‑150 (MQ=255)
           ttGCTGTGGCGCATTGTCACCAACGGCATCACCCTGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGCGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCAGGCGGGCCGACAGCTGATCGGCCCGccag                                                                                                                        <  2:556348/150‑3 (MQ=255)
                                  cGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCACGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCa                                                                                                 <  2:180918/150‑1 (MQ=255)
                                     catcatCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAATCAACTTCGCAAGCAGGCGAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCTCGGCATGa                                                                                              >  1:428482/1‑150 (MQ=255)
                                     catcatCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGa                                                                                              <  2:281884/150‑1 (MQ=255)
                                          tcCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGGCGGGACCGCATGCTGATCGTGCCGCTGAACCAACTTCGAGAGGATGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAAt                                                                                         >  1:624234/1‑150 (MQ=255)
                                           cccGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAAtt                                                                                        <  1:441901/150‑1 (MQ=255)
                                                        aCCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACt                                                                               <  1:274212/146‑1 (MQ=255)
                                                        aCCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACt                                                                               >  2:274212/1‑146 (MQ=255)
                                                                        tAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGTCCCTCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGa                                                           <  2:341332/150‑1 (MQ=255)
                                                                                     cATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCt                                              >  1:511218/1‑150 (MQ=255)
                                                                                          tGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCGTTGACCAGATGGCCGAGCGGCTGCAAg                                         <  2:592287/150‑1 (MQ=255)
                                                                                                 gCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGc                                                                                    >  2:278155/1‑100 (MQ=255)
                                                                                                 gCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGc                                                                                    <  1:278155/100‑1 (MQ=255)
                                                                                                    gCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGACAGCTGATCGGCCCGcca                                                                                                                         <  1:176793/60‑2 (MQ=21)
                                                                                                    gCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGACAGCTGATCGGCCCGcca                                                                                                                         >  2:176793/1‑59 (MQ=21)
                                                                                                      tGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTa                             <  2:428482/150‑1 (MQ=255)
                                                                                                          ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgca                         <  1:180411/150‑1 (MQ=255)
                                                                                                          ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgca                         >  2:368314/1‑150 (MQ=255)
                                                                                                           cAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgcag                        >  1:453445/1‑150 (MQ=255)
                                                                                                            aaCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgcagc                       >  2:697766/1‑150 (MQ=255)
                                                                                                                     aGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgc              <  1:4866/150‑1 (MQ=255)
                                                                                                                        aGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgctgc           <  2:335260/150‑1 (MQ=255)
                                                                                                                          gCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTgcgc         >  1:322470/1‑150 (MQ=255)
                                                                                                                            cAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATg                           <  1:602825/130‑1 (MQ=255)
                                                                                                                            cAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATg                           >  2:602825/1‑130 (MQ=255)
                                                                                                                                 ccaggCGGGCCGACAGCTGATCGGCCCGCCAGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGt  >  1:551005/4‑150 (MQ=255)
                                                                                                                                 ccTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGt  <  2:518175/150‑1 (MQ=255)
                                                                                                                                  cTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTAt                            >  1:650204/1‑123 (MQ=255)
                                                                                                                                  cTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTAt                            <  2:650204/123‑1 (MQ=255)
                                                                                                                                   tGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                               <  1:645404/87‑1 (MQ=255)
                                                                                                                                   tGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                               >  2:645404/1‑87 (MQ=255)
                                                                                                                                        ggCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGAc      <  2:376427/139‑1 (MQ=255)
                                                                                                                                        ggCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGAc      >  1:376427/1‑139 (MQ=255)
                                                                                                                                                     |                                                                                                                                  
GCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGT  >  NC_002947/1846197‑1846475

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