Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,846,339 Δ1 bp 10.8% 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,3390T.10.8% 107.2 / 7.8 37coding (588/1377 nt)pfeS‑IIhistidine kinase
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base . (2/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGT  >  NC_002947/1846191‑1846475
                                                                                                                                                    |                                                                                                                                         
aTCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATc                                                                                                                                          >  1:518175/1‑150 (MQ=255)
      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                                                                               >  2:274212/1‑146 (MQ=255)
                                                              aCCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACt                                                                               <  1:274212/146‑1 (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                         >  2:368314/1‑150 (MQ=255)
                                                                                                                ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgca                         <  1:180411/150‑1 (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)
                                                                                                                                                    |                                                                                                                                         
ATCCGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTGTCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGT  >  NC_002947/1846191‑1846475

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