Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,535,315 G→A 23.3% P185L (CCG→CTG)  envZ ← sensory histidine kinase EnvZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,535,3150GA23.3% 126.7 / 41.5 60P185L (CCG→CTG) envZsensory histidine kinase EnvZ
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base A (7/7);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GGTAACGGAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGA  >  NC_000913/3535218‑3535410
                                                                                                 |                                                                                               
cggcgtgcgcaagtcgtgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAAcgg                                                                                                <  2:1377999‑M2/58‑1 (MQ=255)
 ggcgtgcgcaagtcgtgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAAcggt                                                                                               <  2:626538‑M2/59‑1 (MQ=255)
 gTAACGGAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAAcagt                                                                                               <  2:1345760/99‑1 (MQ=255)
   aaCGGAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAAcagtcg                                                                                             >  1:618307/1‑99 (MQ=255)
    aCGGAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAAcggtcgg                                                                                            >  2:25172/1‑99 (MQ=255)
    aCGGAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAAcagtcgg                                                                                            >  2:1084864/1‑99 (MQ=255)
     tgcgcaagtcgtgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAAcggtcggt                                                                                           >  1:2188521‑M2/37‑99 (MQ=255)
       gAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTc                                                                                         >  1:1903697/1‑99 (MQ=255)
       gAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTc                                                                                         >  2:1470257/1‑99 (MQ=255)
         caagtcgtgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTg                                                                                       >  2:1058668‑M2/33‑99 (MQ=255)
         aCGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTg                                                                                       >  2:952478/1‑99 (MQ=255)
          aagtcgtgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTgg                                                                                      >  1:1094054‑M2/32‑99 (MQ=255)
             aCCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATa                                                                                   >  1:2263633/1‑99 (MQ=255)
               gtgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACg                                                                                 >  1:1648369‑M2/27‑99 (MQ=255)
                tgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGa                                                                                >  1:1671600‑M2/26‑99 (MQ=255)
                tgacttacccccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGa                                                                                >  1:766336‑M2/26‑99 (MQ=255)
                 ccGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGaa                                                                               <  1:1323015/99‑1 (MQ=255)
                 ccGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGaa                                                                               <  1:1640638/99‑1 (MQ=255)
                     aGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATaaa                                                                           <  2:2222931/99‑1 (MQ=255)
                        gcCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAg                                                                        <  1:712712/99‑1 (MQ=255)
                         ccATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGc                                                                       <  1:863318/99‑1 (MQ=255)
                         ccATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGc                                                                       <  2:1580210/99‑1 (MQ=255)
                          ccgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGcc                                                                      >  1:1514266‑M2/16‑99 (MQ=255)
                          cATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGcc                                                                      >  2:392341/1‑99 (MQ=255)
                           cgccatcagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCa                                                                     <  2:1016207‑M2/85‑1 (MQ=255)
                           aTACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCa                                                                     <  1:1472021/99‑1 (MQ=255)
                                 cagcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGcccc                                                               <  1:527546‑M2/91‑1 (MQ=255)
                                 aCGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGcccc                                                               <  1:832541/99‑1 (MQ=255)
                                   gcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGc                                                             >  1:804058‑M2/7‑99 (MQ=255)
                                   gcagcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGc                                                             >  1:1459642‑M2/7‑99 (MQ=255)
                                     agcgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCt                                                           <  2:1973261‑M2/95‑1 (MQ=255)
                                       cgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTAt                                                         <  1:309505‑M2/97‑1 (MQ=255)
                                       cgGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTAt                                                         <  1:1567700‑M2/97‑1 (MQ=255)
                                        gGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATc                                                        <  1:2306243‑M2/98‑1 (MQ=255)
                                          cggcggAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGc                                                      <  2:714794/99‑1 (MQ=255)
                                           ggcggAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGcc                                                     >  1:1888417/1‑99 (MQ=255)
                                            gcggAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCa                                                    >  1:1009026/1‑99 (MQ=255)
                                              ggAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAAt                                                  >  2:1692890/1‑99 (MQ=255)
                                              ggAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAAt                                                  >  1:1994168/1‑99 (MQ=255)
                                              ggAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAAt                                                  >  1:1881012/1‑99 (MQ=255)
                                               gAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATa                                                 >  2:273305/1‑99 (MQ=255)
                                                aataatCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAg                                                <  1:1864520/99‑1 (MQ=255)
                                                aataatCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAg                                                <  1:1472310/99‑1 (MQ=255)
                                                 ataatCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGc                                               >  1:1172035/1‑99 (MQ=255)
                                                      cccTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAAt                                          <  1:1066642/99‑1 (MQ=255)
                                                      cccTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAAt                                          <  1:1827442/99‑1 (MQ=255)
                                                       ccTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATc                                         >  2:2112436/1‑99 (MQ=255)
                                                           tACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCa                                     <  2:1308732/99‑1 (MQ=255)
                                                                aCCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGta                                >  1:1213903/1‑99 (MQ=255)
                                                                aCCTGCAAGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGta                                <  2:1768875/99‑1 (MQ=255)
                                                                      aaGGCTGCGTGTTCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCgg                          >  2:1016936/1‑99 (MQ=255)
                                                                       aGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGa                         >  2:90826/1‑99 (MQ=255)
                                                                       aGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGa                         >  1:1506392/1‑99 (MQ=255)
                                                                        ggCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGaa                        >  2:1039425/1‑99 (MQ=255)
                                                                            gCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGc                    <  2:604849/99‑1 (MQ=255)
                                                                            gCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGc                    >  2:1835933/1‑99 (MQ=255)
                                                                             cGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCg                   <  2:1236648/99‑1 (MQ=255)
                                                                             cGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCg                   <  2:1023901/99‑1 (MQ=255)
                                                                               tgtTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGga                 >  1:270162/1‑99 (MQ=255)
                                                                                 ttCGAGATCGACCAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGgaga               >  1:476861/1‑99 (MQ=255)
                                                                                     agaTCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGaaa           >  2:866869/1‑99 (MQ=255)
                                                                                     agaTCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGaaa           >  2:1448153/1‑99 (MQ=255)
                                                                                        tCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCg        <  2:1961325/99‑1 (MQ=255)
                                                                                           aCCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGccc     >  1:905822/1‑99 (MQ=255)
                                                                                            ccAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCt    <  1:2083796/99‑1 (MQ=255)
                                                                                              aaCGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCtga  <  1:947369/99‑1 (MQ=255)
                                                                                              aaCGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCtga  >  1:187639/1‑99 (MQ=255)
                                                                                                 |                                                                                               
GGTAACGGAACGCACCTCCGAAGCGCCATACTCACGCAGCGGCGGCGGAATAATCCCTTTACCAACCTGCAAGGCTGCGTGTTCGAGATCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGA  >  NC_000913/3535218‑3535410

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