Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 722,788 A→G F543L (TTT→CTT)  kdpD ← sensory histidine kinase KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913722,7880AG97.1% 116.3 / NA 35F543L (TTT→CTT) kdpDsensory histidine kinase KdpD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/17);  minor base T (0/1);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTCGCGGTAACGGGCTACCCGCGCCTGATAACGCAC  >  NC_000913/722669‑722923
                                                                                                                       |                                                                                                                                       
gcCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGtt                                                                                                                     >  1:499013/1‑140 (MQ=255)
      aTCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCtgg                                                                                                               >  2:168545/1‑140 (MQ=255)
      aTCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCtgg                                                                                                               >  1:675849/1‑140 (MQ=255)
       tCAAAACTCCACTGCGCCATGGCATCGTCCCTCGGAGTCATTCCTTGCCGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCtggt                                                                                                              <  2:349512/140‑1 (MQ=255)
         aaaaCTCCACTGCGCGAGGGCATCGTCCCACGGCGACATTCCTTGCGGCTGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCtggtgg                                                                                                            <  2:333026/140‑1 (MQ=255)
                                  tCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGTGATATCCTGCGGACTGCGGcc                                                                                   >  1:57283/1‑140 (MQ=255)
                                  tCCCACGGCGTCATTACTTGCGGATTTGTCAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGcc                                                                                   <  1:310759/140‑1 (MQ=255)
                                   cccACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGTGATATCCTGCGGACTGCGGccc                                                                                  <  2:57283/140‑1 (MQ=255)
                                      aCGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTg                                                                                               <  2:509118/124‑1 (MQ=255)
                                      aCGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTg                                                                                               >  1:509118/1‑124 (MQ=255)
                                      aCGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCAcc                                                                               >  2:172433/1‑140 (MQ=255)
                                      aCGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGCTCGGTGGCGGCAGCGATATCCAGCGGACTGCGGCCCAcc                                                                               >  2:310759/1‑140 (MQ=255)
                                       cGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCg                                                                              >  1:681763/1‑140 (MQ=255)
                                            tCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCaga                                                                         >  1:834747/1‑140 (MQ=255)
                                            tCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCaga                                                                         <  2:675849/140‑1 (MQ=255)
                                                   ttGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCg                                                                  <  2:499013/140‑1 (MQ=255)
                                                    tGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGa                                                                 <  2:197293/140‑1 (MQ=255)
                                                    tGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGa                                                                 <  1:45878/140‑1 (MQ=255)
                                                    tGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGa                                                                 <  1:172433/140‑1 (MQ=255)
                                                      cGGATGTGTTAACAGCTGCAATCTACCGTCGTCATCGTGCAACAACGCCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACa                                                               <  1:240378/140‑1 (MQ=255)
                                                       ggATGTGTTAACGGCCGCAATTCACCGTTCTCATCGGCCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACAt                                                              <  1:229179/140‑1 (MQ=255)
                                                           gtgtTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCa                                                          >  2:457486/1‑140 (MQ=255)
                                                             gtTAACGGCTGCAATTTACCGTTGTGATCTGGCAACATCCCCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCata                                                        <  1:290229/140‑1 (MQ=255)
                                                                   ggCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTTGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACAt                                                              >  1:479200/1‑128 (MQ=255)
                                                                   ggCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTTGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACAt                                                              <  2:479200/128‑1 (MQ=255)
                                                                      tGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATTAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCAAGAACTGTCGACAGTTCATATAAGTGTCg                                               >  1:302426/1‑140 (MQ=255)
                                                                        cAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCgtg                                             >  2:839595/1‑140 (MQ=255)
                                                                                    gTCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTcgc                                 >  1:498085/1‑140 (MQ=255)
                                                                                     tCATCGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTATCGACATTTCATATTAGTGCCCTGTGCGTTGCTcgcg                                >  1:282850/1‑140 (MQ=255)
                                                                                         cGGGCAACAACACCTGACTGCGTGCATTAAGTGTGGAGGGAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGTCGTGTGCGTTGCTCGCGGTaa                            <  2:302426/140‑1 (MQ=255)
                                                                                         cGGGCAACAACACCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTCGCGGTaa                            <  2:681763/140‑1 (MQ=255)
                                                                                               acaacaCCTGACTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCGTTCGACATGTCAGATAGGTGCCGTGGGCGTGGCTCGCGGTAACGGGCt                      >  2:235686/1‑140 (MQ=255)
                                                                                                      cTGACTGCGGGCATGAATCGTGAAGGCAATATATTTTTCACTGGTGGCTGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTCGCGGTAACGGGCTACCcgcg               <  2:282850/140‑1 (MQ=255)
                                                                                                         aCTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTCGCGGTAACGGGCTACCcgcg               <  2:188890/137‑1 (MQ=255)
                                                                                                         aCTGCGGGCATGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTCGCGGTAACGGGCTACCcgcg               >  1:188890/1‑137 (MQ=255)
                                                                                                                   tGAAGCGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTCGCGGTAACGGGCTACCCGCGCCTGATAacgcac  <  1:839595/140‑1 (MQ=255)
                                                                                                                       |                                                                                                                                       
GCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGAGGCAATAAATTGTTCGCTGGTGGCAGCGATATCCTGCGGACTGCGGCCCACCGCCAGAGCTTTCGACATTTCATATAAGTGCCGTGTGCGTTGCTCGCGGTAACGGGCTACCCGCGCCTGATAACGCAC  >  NC_000913/722669‑722923

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