Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,628,174 T→C 8.1% V96A (GTC→GCC)  mucA → sigma factor AlgU negative regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,628,1740TC 8.1% 134.1 / 3.1 49V96A (GTC→GCC) mucAsigma factor AlgU negative regulator
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

CCATGCACAAGGAACTGCTGCTGCCTAACCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGAAAGCAGTGAGC  >  NC_002947/1628030‑1628320
                                                                                                                                                 |                                                                                                                                                    
ccATGCACAAGGAACTGCTGCTGCCTAACCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTcgcg                                                                                                                                                  <  5:32521/149‑1 (MQ=255)
    gCACAAGGAACTGCTGCTGCCTAACCTGGATATCGCCTCGGCGGTGTCCGCGGCGGGTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑Gt                                                                                                                                              >  1:238482/1‑148 (MQ=255)
          ggAACTGCTGCTGCCTAACCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCg                                                                                                                                      <  5:282738/149‑1 (MQ=255)
              ctgctgctgCCTAACCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGgtgt                                                                                                                                  <  3:203394/149‑1 (MQ=255)
                    ctgcCTAACCTGGATATCGCCTCGGCGGTGTCGGAGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCTCAGCATTGGCCGTCTGGCGGTAGACGCCTCGGTTACCGTCGCT‑‑GTACTGGCCGGTGTGCGTAt                                                                                                                            <  7:56566/149‑1 (MQ=255)
                    ctgcCTAACCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTAt                                                                                                                            >  8:77834/1‑149 (MQ=255)
                       ccTAACCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGt                                                                                                                          <  5:14304/148‑1 (MQ=255)
                          aaCCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGc                                                                                                                                >  6:115856/1‑139 (MQ=255)
                          aaCCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGc                                                                                                                                <  5:115856/139‑1 (MQ=255)
                           aCCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACaa                                                                                                                      <  5:281060/148‑1 (MQ=255)
                           aCCTGGATATCGCCTCGGCGGTGTACGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAAc                                                                                                                     <  6:222893/149‑1 (MQ=255)
                                    tCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGGTCGTATGGACAACCAGGACGa                                                                                                             >  7:120539/1‑148 (MQ=255)
                                    tCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCAGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCAGCCTCGGTTACCGTCGCG‑‑GTACTGGACGGTGTGAGTATGTACAAGAAGGACGa                                                                                                             >  7:124580/1‑148 (MQ=255)
                                      gCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAAt                                                                                                          >  1:4897/1‑149 (MQ=255)
                                           ggcggTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCAc                                                                                                       <  7:77834/147‑1 (MQ=255)
                                           ggcggTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCAc                                                                                                       <  2:238482/147‑1 (MQ=255)
                                                tgtCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTg                                                                                                  >  2:35884/1‑147 (MQ=255)
                                                tgtCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCa                                                                                                <  1:35884/149‑1 (MQ=255)
                                                    cgcgGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGc                                                                                            >  2:70084/1‑149 (MQ=255)
                                                      cggc‑tCTGGCCTACTAAGACGTGCCAGCCAAGCTCAATAAACGCCCTTGTCGAAGCAATGGCCGTCTGTCGGTAGACGCCTCTGTTACCGTCGCG‑‑GTACTGTCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCtt                                                                                          <  8:124580/149‑1 (MQ=255)
                                                      cggc‑gCTGTCCGACGAATCCGTGCCATCCAAGGGCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGTTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCtt                                                                                          <  8:120539/149‑1 (MQ=255)
                                                       ggc‑gCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTtg                                                                                         <  8:89941/149‑1 (MQ=255)
                                                           gcTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTtgctg                                                                                      >  3:56803/1‑149 (MQ=255)
                                                            cTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTtgctgc                                                                                     <  2:21036/149‑1 (MQ=255)
                                                                   acgaAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGc                                                                                                                                >  5:78473/1‑99 (MQ=255)
                                                                   acgaAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGc                                                                                                                                <  6:78473/99‑1 (MQ=255)
                                                                    cgaAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCagcagc                                                                             <  1:59674/149‑1 (MQ=255)
                                                                        gCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGt                                                                         <  1:70084/149‑1 (MQ=255)
                                                                            tGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGccc                                                                                   <  8:218751/135‑1 (MQ=255)
                                                                            tGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGccc                                                                                   >  7:218751/1‑135 (MQ=255)
                                                                                  ccaAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGc                                                               >  2:31326/1‑149 (MQ=255)
                                                                                          aaGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAgtgt                                                       <  1:31326/149‑1 (MQ=255)
                                                                                           aGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCGGAGAACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCAGAgtg                                                        >  1:87729/1‑149 (MQ=255)
                                                                                           aGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGtgtg                                                      >  8:209991/1‑149 (MQ=255)
                                                                                            gAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGCCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGc                                                     >  5:117149/1‑149 (MQ=255)
                                                                                            gAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGCCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGc                                                     >  5:112468/1‑149 (MQ=255)
                                                                                                 gCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCcac                                                  >  7:258978/1‑147 (MQ=255)
                                                                                                  cccTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAg                                               >  7:52655/1‑149 (MQ=255)
                                                                                                     ttGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAgcg                                             >  2:215595/1‑148 (MQ=255)
                                                                                                                     cGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGtttt                            <  8:258978/149‑1 (MQ=255)
                                                                                                                           gCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAgg                      >  4:211801/1‑149 (MQ=255)
                                                                                                                             ggTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCt                    <  7:209991/149‑1 (MQ=255)
                                                                                                                                 gccgccTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCtata                 <  4:56803/148‑1 (MQ=255)
                                                                                                                                     ccTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGAGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGaa            <  1:215595/149‑1 (MQ=255)
                                                                                                                                     ccTCGGTTACCGCCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGaa            <  6:117149/149‑1 (MQ=255)
                                                                                                                                     ccTCGGTTACCGCCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGaa            <  6:112468/149‑1 (MQ=255)
                                                                                                                                         cgTTACGGTCGCG‑‑GTACTGGCAGGTGTGCGTCTGTACAACCTGGACGAAATCCCGGTTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGAAAGc         <  1:125504/147‑1 (MQ=255)
                                                                                                                                               cGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGAAAGCAGTGAg   >  5:136707/1‑148 (MQ=255)
                                                                                                                                                gTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGAAAGCAGTGAGc  >  7:137062/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCATGCACAAGGAACTGCTGCTGCCTAACCTGGATATCGCCTCGGCGGTGTCCGCGGC‑GCTGGCCGACGAAGCCGTGCCAGCCAAGGTCAAGAAAGGCCCTTGGCGCAGCATTGGCCGTCTGGCGGTAGCCGCCTCGGTTACCGTCGCG‑‑GTACTGGCCGGTGTGCGTATGTACAACCAGGACGAAATCACGGGTGCAGAGCTTGCTGCCCAGCAGCCTGTCCAACAGGGCCTGAGTGTGCCACAAGCGCAAGGCCCGGCCGTTTTGGCAGGCTATAGTGAAAGCAGTGAGC  >  NC_002947/1628030‑1628320

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