Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,295,019 Δ1 bp 13.6% intergenic (+29/+41) slyB → / ← nhaA‑I outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,295,0190A.13.6% 124.4 / 15.6 44intergenic (+29/+41)slyB/nhaA‑Iouter membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base . (3/3);  total (24/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 = 6.32e-01

TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGT  >  NC_002947/1294882‑1295163
                                                                                                                                         |                                                                                                                                                 
tgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCa                                                                                                                                       >  2:919791/1‑150 (MQ=255)
    agcaTGCGCGCCTATGTGCAGGCCGTGGAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCatgtcg                                                                                                                                  <  2:961349/150‑6 (MQ=255)
    agcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCAACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCg                                                                                                                                  >  1:783601/1‑150 (MQ=255)
     gcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCgg                                                                                                                                 <  1:475828/150‑1 (MQ=255)
               cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                       <  1:64784/150‑1 (MQ=255)
                tATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTCCTGCCAGTCCGGGTTTTTTCAt                                                                                                                      <  2:955022/150‑1 (MQ=255)
                             tGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGt                                                                                                         >  2:1015770/1‑150 (MQ=255)
                                    gaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCaa                                                                                                  <  1:541326/150‑1 (MQ=255)
                                     aaCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCa                                                                                                                <  1:900964/135‑1 (MQ=255)
                                     aaCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCa                                                                                                                >  2:900964/1‑135 (MQ=255)
                                     aaCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAAt                                                                                                 >  1:206099/1‑150 (MQ=255)
                                     aaCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCggg                                                                                                                                <  2:56018/119‑1 (MQ=38)
                                     aaCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCggg                                                                                                                                >  1:56018/1‑119 (MQ=38)
                                        gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCt                                                                                              >  1:442699/1‑150 (MQ=255)
                                        gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCt                                                                                              >  1:716357/1‑150 (MQ=255)
                                          aaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTgg                                                                                            >  1:1210111/1‑150 (MQ=255)
                                                                gTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGt                                                                      <  1:1015770/150‑1 (MQ=255)
                                                                         atGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTACCGGACCAATGCTGCACCCCCCAAGTCGGATAATGCTGGTTTTGCGGCTGGCCATTGCGGGGACGGCGTa                                                             >  1:472933/1‑150 (MQ=255)
                                                                          tGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAg                                                            >  1:1247254/1‑150 (MQ=255)
                                                                           gACTGTGGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGcc                                                                                                                   >  2:479340/1‑94 (MQ=255)
                                                                           gACTGTGGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGcc                                                                                                                   <  1:479340/94‑1 (MQ=255)
                                                                           gACTGTAGACGGTACCAGTCGCATTTCGCACTCACCGTAGTGCAGTAATAAAATACCCCGACATGGCTGCTGCCAGTCGCGGTTTTTTCATGCCTCAGTCAAGTCAGGCAATGCTTGTTGTGCGGCTGGCCATTACCACGACGGCGTAGc                                                           <  2:633898/150‑1 (MQ=255)
                                                                            aCTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGcc                                                          <  2:783601/150‑1 (MQ=255)
                                                                                  gACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAc                                                    >  1:758974/1‑150 (MQ=255)
                                                                                  gACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAc                                                    <  1:854658/150‑1 (MQ=255)
                                                                                        aCCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcagc                                              <  2:716357/150‑1 (MQ=255)
                                                                                            gTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAg                                          >  1:328201/1‑150 (MQ=255)
                                                                                            gTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTCGCCATTGCGGTGACGCCCTAGCCGATCACCgcacccacg                                          >  2:205863/1‑144 (MQ=255)
                                                                                               gcgTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAGgaac                                       <  2:758974/150‑1 (MQ=255)
                                                                                                  tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaa                                    >  1:528522/1‑150 (MQ=255)
                                                                                                       cACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGt                               >  2:1184633/1‑150 (MQ=255)
                                                                                                         cTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTaa                             <  2:1247254/150‑1 (MQ=255)
                                                                                                            aCCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGc                                                                >  2:751094/1‑112 (MQ=255)
                                                                                                            aCCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGc                                                                <  1:751094/112‑1 (MQ=255)
                                                                                                                    gCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCGGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATg                  >  1:1122521/1‑150 (MQ=255)
                                                                                                                       gTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCgg               <  2:206099/150‑1 (MQ=255)
                                                                                                                        taataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                   <  1:90905/129‑1 (MQ=255)
                                                                                                                        taataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                   >  2:90905/1‑129 (MQ=255)
                                                                                                                               aaaCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGcc       >  1:799592/1‑150 (MQ=39)
                                                                                                                                 aCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt              >  2:27698/1‑141 (MQ=39)
                                                                                                                                 aCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt              <  1:27698/141‑1 (MQ=39)
                                                                                                                                  ccccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacg                                      <  2:872591/116‑1 (MQ=255)
                                                                                                                                  ccccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacg                                      >  1:872591/1‑116 (MQ=255)
                                                                                                                                    ccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGt  >  1:46488/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGT  >  NC_002947/1294882‑1295163

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: