Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,966 G→A 10.4% G354G (GGC→GGT)  PP_4224 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,796,9660GA10.4% 103.3 / 3.0 40G354G (GGC→GGT) PP_4224sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base A (2/2);  total (21/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 = 9.77e-01

ACCGCAACCACCACTCTGCCCATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTCGCCTCCGCCGGCATC  >  NC_002947/4796822‑4797103
                                                                                                                                                |                                                                                                                                         
aCCGCAACCACCACTCTGCCCATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCgccgcc                                                                                                                                      >  2:294902/1‑150 (MQ=255)
          ccacTCTGCCCATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAg                                                                                                                            <  2:716371/150‑1 (MQ=255)
               ctGCCCATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCAtt                                                                                                                       <  1:564026/150‑1 (MQ=255)
                tGCCCATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCAGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTg                                                                                                                      >  2:488212/1‑150 (MQ=255)
                    cATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATcgcc                                                                                                                                         >  1:138736/1‑127 (MQ=255)
                    cATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATcgcc                                                                                                                                         <  2:138736/127‑1 (MQ=255)
                            cAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTg                                                                                                          >  2:241242/1‑150 (MQ=255)
                                       ggTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGc                                                                                               <  1:241242/150‑1 (MQ=255)
                                          cAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCCTTCGGTCGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCAc                                                                                            <  1:572342/150‑1 (MQ=255)
                                           aGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCAcc                                                                                           >  2:429796/1‑150 (MQ=255)
                                            ggTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATATGCACCTCCAACTCCACGGCTTCATTGACCAGATCAACAAAGGCACGTGCCAGCAGGTTGCGCAGCGCTCCCATCg                                                                                          >  1:354935/1‑150 (MQ=255)
                                                 tCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAAc                                                                                     >  2:113023/1‑150 (MQ=255)
                                                  ccAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGc                                                                                         <  1:287272/145‑1 (MQ=255)
                                                  ccAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGc                                                                                         >  2:287272/1‑145 (MQ=255)
                                                  ccAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACg                                                                                    <  1:429796/150‑1 (MQ=255)
                                                      gTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGc                                                                                >  1:132851/1‑150 (MQ=255)
                                                                gCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATcccc                                                                      >  2:547680/1‑150 (MQ=255)
                                                                            ctgtGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTCGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGc                                                          >  2:187394/4‑150 (MQ=255)
                                                                            ccATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGt                                                                                                            >  2:447359/1‑100 (MQ=255)
                                                                            ccATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGt                                                                                                            <  1:447359/100‑1 (MQ=255)
                                                                            ccATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGc                                                          <  1:294902/150‑1 (MQ=255)
                                                                                        aCCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCctct                                              >  2:466364/1‑150 (MQ=255)
                                                                                        aCCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGCCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCctct                                              >  1:710848/1‑150 (MQ=255)
                                                                                           aTGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCACCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGc                                           >  1:681591/1‑150 (MQ=255)
                                                                                                   cACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCaa                                   <  1:488212/150‑1 (MQ=255)
                                                                                                   cACCTGGCCGTCCACCCCCTGCACCTCCCACTCCACGGCTTCATGCCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGGCGGCATCCCCAGGTAAATATGCGTCGCCTCCTCTGGCAAACGCaa                                   >  2:528998/1‑150 (MQ=255)
                                                                                                               cGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCg                       >  2:207244/1‑150 (MQ=255)
                                                                                                                 ccATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGAt                     >  2:223845/1‑150 (MQ=255)
                                                                                                                  cATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGAtt                    <  2:132851/150‑1 (MQ=255)
                                                                                                                  cATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGAtt                    <  1:113023/150‑1 (MQ=255)
                                                                                                                  cATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGAtt                    <  1:547680/150‑1 (MQ=255)
                                                                                                                    tCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGCCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGt                  <  2:710848/150‑1 (MQ=255)
                                                                                                                     cTGCACCTCTCACCCTCCGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCACCGGTGCAGGCATCCCCAGGTAAATCTTCGTCGCCTCCTCCGGCAAACGCAACACGATATGCCGATTGTc                 <  2:245206/150‑1 (MQ=255)
                                                                                                                     cTGCACCTCCAACTCCACGGCTTCATCACCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTc                 <  2:681591/150‑1 (MQ=255)
                                                                                                                         aCCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTc                 <  2:402550/146‑1 (MQ=255)
                                                                                                                         aCCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTc                 >  1:402550/1‑146 (MQ=255)
                                                                                                                           cTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTCGCCTcc           >  1:648277/1‑150 (MQ=255)
                                                                                                                            tCCAACTCCACGGCTTCATCACCGCCATGGCGCAAGGCAtt                                                                                                                       >  1:283207/1‑41 (MQ=255)
                                                                                                                            tCCAACTCCACGGCTTCATCACCGCCATGGCGCAAGGCAtt                                                                                                                       <  2:283207/41‑1 (MQ=255)
                                                                                                                            tCCAACTCAACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTGGCCTccg          >  1:157089/1‑150 (MQ=255)
                                                                                                                              cAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGTCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTCGCCTccgcc        <  2:120717/150‑1 (MQ=255)
                                                                                                                                  tCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTCGCCTCCGCCGGCa    >  1:204956/1‑150 (MQ=255)
                                                                                                                                    cACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTCGCCTCCGCCGGCATc  <  2:157089/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                         
ACCGCAACCACCACTCTGCCCATTGCGCCAGAAGCGCTCGGTCAGGTGCTCCAGGTCCGCCTCAGCAATGCCTGGCCCATGGTCACGCACCATGAAGCCCACCTGGCCGTCCGCCATCTGCACCTCCAACTCCACGGCTTCATCGCCGCCATGGCGCAAGGCATTGTCCAGCAGGTTGCGCAGCGCTGCCACCGCCAACGGTGCCGGCATCCCCAGGTAAATCTGCGTCGCCTCCTCTGGCAAACGCAACAGGATACGCCGATTGTCGCCTCCGCCGGCATC  >  NC_002947/4796822‑4797103

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