Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,846,297 Δ1 bp 9.1% coding (588/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,2970T. 9.1% 129.5 / 8.1 44coding (588/1377 nt)pfeS‑IIhistidine kinase
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base . (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

ATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAA  >  NC_002947/1846149‑1846442
                                                                                                                                                    |                                                                                                                                                   
aTCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGAt                                                                                                                                                     <  5:72687/149‑1 (MQ=255)
 tCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGAt                                                                                                                                                     >  6:10798/1‑148 (MQ=255)
 tCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGAt                                                                                                                                                     >  2:157457/1‑148 (MQ=255)
      gCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGTCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTg                                                                                                                                               <  8:53067/149‑1 (MQ=255)
       ccAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGt                                                                                                                                             <  7:73806/149‑1 (MQ=255)
          gCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCgg                                                                                                                                          >  5:270270/1‑149 (MQ=255)
             ccGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGGACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGACAGCTGATCGGccc                                                                                                                                       <  5:198547/149‑1 (MQ=255)
               gCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGAGCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGc                                                                                                                                     <  3:256811/149‑1 (MQ=255)
                   gCTGTGGCGCATTGTCACCAACAGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCTTGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCt                                                                                                                                   <  6:9272/147‑1 (MQ=255)
                    cTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGa                                                                                                                                >  7:6433/1‑149 (MQ=255)
                           gcATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTcgcg                                                                                                                          >  7:76257/1‑148 (MQ=255)
                             aTTGTCACCAACGGCATCATCCCGGGGCTGTGCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGAGCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCAc                                                                                                                       <  5:191501/149‑1 (MQ=255)
                                     cAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGccc                                                                                                               >  8:50194/1‑149 (MQ=255)
                                      aaCGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCa                                                                                                                                                   >  5:9421/1‑113 (MQ=255)
                                      aaCGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCa                                                                                                                                                   <  6:9421/113‑1 (MQ=255)
                                           catcatCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCAt                                                                                                          >  2:39828/1‑148 (MQ=255)
                                            atcatcCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGTCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATAGCCCGGCATGa                                                                                                       >  7:108312/1‑149 (MQ=255)
                                            atcatcCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGa                                                                                                       >  6:165619/1‑149 (MQ=255)
                                                 cccGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAAt                                                                                                  >  6:16876/1‑149 (MQ=255)
                                                      ggCTGTTCAGCTTGCTGTTGTGTGTAGGCGTGTACCGCATGCTGATCTTGCCGGTGAACCAAGTTCGCGAGGAGGCCAATGCGTGGCGGGCCGATGAGCTTTCGGCCCGCCTGGACTCGCGCACCATCGCGCGGGATGACGAATTGGGc                                                                                             <  8:6433/149‑1 (MQ=255)
                                                         tGTTCACCTTGCTGTTGTGTGTAGGTCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATAGCCCGGCATGACGAATTGGGCGaa                                                                                          <  8:108312/149‑1 (MQ=255)
                                                           ttCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACt                                                                                        >  6:270848/1‑149 (MQ=255)
                                                              aCCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGgc                                                                                     <  5:165619/149‑1 (MQ=255)
                                                                cTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGgcgc                                                                                   <  5:270848/149‑1 (MQ=255)
                                                                  tGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGt                                                                                 >  7:104811/1‑149 (MQ=255)
                                                                  tGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGt                                                                                 >  7:94950/1‑149 (MQ=255)
                                                                   gCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGt                                                                                 <  1:39828/148‑1 (MQ=255)
                                                                       ttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCtt                                                                            >  2:154673/1‑149 (MQ=255)
                                                                          tgtgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGCACGAATTGGGCGAACTGGCGCGTGCCttt                                                                           <  1:157457/148‑1 (MQ=255)
                                                                            tgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                                        <  8:104811/148‑1 (MQ=255)
                                                                                ggCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGAt                                                                   <  5:16876/149‑1 (MQ=255)
                                                                                                  aTCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCgg                                                         >  1:178787/1‑141 (MQ=255)
                                                                                                  aTCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCgg                                                         <  2:178787/141‑1 (MQ=255)
                                                                                                              aaCCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTAt                                     >  1:175334/1‑149 (MQ=255)
                                                                                                                         cgAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                          >  4:174560/1‑149 (MQ=255)
                                                                                                                          gAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGtt                         >  4:155100/1‑149 (MQ=255)
                                                                                                                          gAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCCCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGtt                         >  7:181958/1‑149 (MQ=255)
                                                                                                                                gCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACTATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgctgc                    >  4:240196/1‑148 (MQ=255)
                                                                                                                                gCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACTATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgctgc                    <  3:240196/148‑1 (MQ=255)
                                                                                                                                  cAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                          >  7:212622/1‑140 (MQ=255)
                                                                                                                                  cAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                          <  8:212622/140‑1 (MQ=255)
                                                                                                                                    aTGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGAc               <  7:50194/149‑1 (MQ=255)
                                                                                                                                         tGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTc          >  6:24197/1‑149 (MQ=255)
                                                                                                                                           gCGGGCCGACAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCt             <  6:146195/144‑1 (MQ=255)
                                                                                                                                           gCGGGCCGACAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCt             >  5:146195/1‑144 (MQ=255)
                                                                                                                                           gCGGGCCGACAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTcc         >  1:7223/1‑148 (MQ=255)
                                                                                                                                                 cGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGaaa  <  8:94950/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAA  >  NC_002947/1846149‑1846442

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