Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,672,777 C→A 7.9% E52* (GAG→TAG)  PP_1467 ← NhaP‑type Na+(K+)/H+ antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,672,7770CA 7.9% 148.8 / 2.7 51E52* (GAG→TAG) PP_1467NhaP‑type Na+(K+)/H+ antiporter
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base A (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GACCAGTGCGCCCACCGTGACCAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATCTGCTGTTC  >  NC_002947/1672631‑1672924
                                                                                                                                                  |                                                                                                                                                   
gACCAGTGCGCCCACCGGGACCAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCggg                                                                                                                                                  <  2:377669/150‑1 (MQ=255)
     gTGCGCCCACCGTGACCAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAACGCCTCTGTGTCGa                                                                                                                                             >  2:932576/1‑150 (MQ=255)
          cccACCGTGACCAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCAt                                                                                                                                        <  1:257619/150‑1 (MQ=255)
            cACCGTGACCAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATcc                                                                                                                                      >  1:323385/1‑150 (MQ=255)
                  gACCAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAgcgc                                                                                                                                >  1:1010052/1‑150 (MQ=255)
                    ccAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTg                                                                                                                              >  1:505360/1‑150 (MQ=255)
                       ggCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCAACGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGc                                                                                                                           <  1:514830/150‑1 (MQ=255)
                         cGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCg                                                                                                                         <  2:1010052/150‑1 (MQ=255)
                           aTGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGc                                                                                                                       >  2:803181/1‑150 (MQ=255)
                             gCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCaa                                                                                                                     <  2:939060/150‑1 (MQ=255)
                                ccacGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGa                                                                                                                  >  2:387460/1‑150 (MQ=255)
                                 cacGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCaa                                                                                                                     <  2:1035234/146‑1 (MQ=255)
                                 cacGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCaa                                                                                                                     >  1:1035234/1‑146 (MQ=255)
                                 cacGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTAGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGAt                                                                                                                 >  1:1152817/1‑150 (MQ=255)
                                               cgcgCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAacaaca                                                                                                   <  1:387460/150‑1 (MQ=255)
                                               cgcgCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTAGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAacaaca                                                                                                   <  2:1152817/150‑1 (MQ=255)
                                                gcgcCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACaacag                                                                                                  >  1:444504/1‑150 (MQ=255)
                                                  gcCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACaacagaa                                                                                                <  2:734106/150‑1 (MQ=255)
                                                  gcCACTGCGACAGGTGCAGGGACAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACaacagaa                                                                                                >  2:583271/1‑150 (MQ=255)
                                                        gCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAacagaaacagaa                                                                                          <  1:803181/150‑1 (MQ=255)
                                                           aCAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTAGACGATTCCGAGCGCTGGGCCGGCAAGTATCCCGGCCAACAACAGAAACAGAATTg                                                                                       <  1:932576/150‑1 (MQ=255)
                                                            cAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGc                                                                                      <  2:760675/150‑1 (MQ=255)
                                                            cAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGc                                                                                      <  2:1205004/150‑1 (MQ=255)
                                                              ggTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCg                                                                                    <  1:583271/150‑1 (MQ=255)
                                                                        cAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAg                                                                          <  1:1027558/150‑1 (MQ=255)
                                                                          ggCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCcaa                                                                                                        >  2:316581/1‑118 (MQ=255)
                                                                          ggCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCcaa                                                                                                        <  1:316581/118‑1 (MQ=255)
                                                                               ccTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCgccag                                                                   <  2:545442/150‑1 (MQ=255)
                                                                                 ttCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCACGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCgccaggc                                                                 >  2:208840/1‑150 (MQ=255)
                                                                                   cGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCgccaggcca                                                               <  2:323385/150‑1 (MQ=255)
                                                                                     aaCAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGgccagc                                                             <  2:700205/150‑1 (MQ=255)
                                                                                      aCAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGgccagcc                                                            >  1:905576/1‑150 (MQ=255)
                                                                                           aTCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAGCAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTgg                                                       >  1:567582/1‑150 (MQ=255)
                                                                                           aTCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGCCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGCAACAGAATTGCCGGCAACTTCAGACGCCAGGCCAGCCACTgg                                                       >  2:616366/1‑150 (MQ=255)
                                                                                               gcgcCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAgg                                                   >  1:1074898/1‑150 (MQ=255)
                                                                                                      cGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAgcgc                                            >  1:1240981/1‑150 (MQ=255)
                                                                                                               cacCAGCGGCATCAGCAATGGGCCGAACAGCGCCTAGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGcc                                   >  2:133172/1‑150 (MQ=255)
                                                                                                                  cAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCACCTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGAt                                <  1:616366/150‑1 (MQ=255)
                                                                                                                          tCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCca                        >  1:1211512/1‑150 (MQ=255)
                                                                                                                            aGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCcaca                      <  2:505360/150‑1 (MQ=255)
                                                                                                                                  gggCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGa                <  1:1030240/150‑1 (MQ=255)
                                                                                                                                  gggCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAAGTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGa                <  2:305653/150‑1 (MQ=255)
                                                                                                                                  gggCCGAACAGCGCCTAGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGa                <  1:133172/150‑1 (MQ=255)
                                                                                                                                         aCAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATct         >  1:812405/1‑150 (MQ=255)
                                                                                                                                           aGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATctgc       <  2:905576/150‑1 (MQ=255)
                                                                                                                                              gcCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATctg        >  2:91867/1‑146 (MQ=255)
                                                                                                                                              gcCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATctg        <  1:91867/146‑1 (MQ=255)
                                                                                                                                               ccTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATct         <  1:169612/144‑1 (MQ=255)
                                                                                                                                               ccTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATct         >  2:169612/1‑144 (MQ=255)
                                                                                                                                                cTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAGCAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATCTGCTGTTc  <  2:567582/150‑1 (MQ=255)
                                                                                                                                                cTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATCTGCTGTTc  >  1:690654/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCAGTGCGCCCACCGTGACCAGGCGATGCACCACGCTGCCGATCTCGCGCCACTGCGACAGGTGCAGGGTCAGGCTGCCTTCGAACAGAATCAGCGCCACCGACAGCGACACCAGCGGCATCAGCAATGGGCCGAACAGCGCCTCGGGGTCGAGCCATCCCAGCGCTGGGCCGGCAAGGATCCCGGCCAACAACAGAAACAGAATTGCCGGCAACTTCAGGCGCCAGGCCAGCCACTGGCAGGCGAGCGCCGCGGCGCCGATGCCGCCCACACTCAGAAGGATCTGCTGTTC  >  NC_002947/1672631‑1672924

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: ‑