Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,628,995 C→A 8.5% L169I (CTC→ATC)  mucB → sigma factor AlgU regulator MucB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,628,9950CA 8.5% 144.7 / 3.8 47L169I (CTC→ATC) mucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base A (2/2);  total (20/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 = 7.76e-01

TACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCCGTCGCGGGTT  >  NC_002947/1628854‑1629142
                                                                                                                                             |                                                                                                                                                   
tACGATTTGAGTGTGGCAGGCAAGTCACGGGTCACAGGGCGCGATGACGTGATAATTACGCTCACGCCACGCGACCAGAATCGCTATGCGTTCGAATTGCAACTGGACCGTCACTAGGTTTTGCCGCTGCGTTCGTTGATGCTCAACg                                                                                                                                               >  7:45903/1‑148 (MQ=255)
  cGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGAc                                                                                                                                             <  1:149729/148‑1 (MQ=255)
        gAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAgggg                                                                                                                                       <  6:139041/148‑1 (MQ=255)
          gtgtGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGATGATGCTCAACGACAAGGGGCAg                                                                                                                                    <  2:208020/149‑1 (MQ=255)
             tgGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTg                                                                                                                                 >  1:148717/1‑149 (MQ=255)
             tgGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGATTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTg                                                                                                                                 >  7:71021/1‑149 (MQ=255)
               gcaggcaAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCt                                                                                                                               <  8:246663/149‑1 (MQ=255)
                 aggcaAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTg                                                                                                                              >  7:96039/1‑148 (MQ=255)
                   gcaAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCAGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGaa                                                                                                                           >  4:106332/1‑149 (MQ=255)
                          aCGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTcc                                                                                                                    <  3:106332/149‑1 (MQ=255)
                          aCGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTc                                                                                                                     >  1:146949/1‑148 (MQ=255)
                          aCGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTc                                                                                                                     <  8:63335/148‑1 (MQ=255)
                              gTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGAt                                                                                                                <  2:176783/149‑1 (MQ=255)
                                   cGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGAccc                                                                                                           <  3:142981/149‑1 (MQ=255)
                                       cgcgATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGtct                                                                                                       <  8:71021/149‑1 (MQ=255)
                                       cgagATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGtct                                                                                                       <  8:128576/146‑1 (MQ=255)
                                                 tGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGAt                                                                                             <  4:96392/149‑1 (MQ=255)
                                                               acgccacgcGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGa                                                                               >  6:128843/1‑149 (MQ=255)
                                                                    acgcGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACg                                                                             >  1:135773/1‑146 (MQ=255)
                                                                        gACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCTTAACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGTGCAGTTGCTGGAACGTTTCCAGATGACCCGTCTCGATATAGATGAAATACCCAGCGATGCAGACtt                                                                      <  8:45903/149‑1 (MQ=255)
                                                                            aGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTgcgc                                                                  <  8:250575/149‑1 (MQ=255)
                                                                             gCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAATGCTTCCAGATGACCCGTCTCGATACCGATGAAAt                                                                                        <  3:81165/126‑1 (MQ=255)
                                                                             gCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAATGCTTCCAGATGACCCGTCTCGATACCGATGAAAt                                                                                        >  4:81165/1‑126 (MQ=255)
                                                                                  gCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGt                                                            <  2:148717/149‑1 (MQ=255)
                                                                                          ttCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGctgcctg                                                    >  2:201135/1‑149 (MQ=255)
                                                                                               agtGCACCGGGACCGTCACTCTGGTTTGCCGCTGCGTGCGTTGATGCTCAACGACAAGGGTCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGAGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGc                                               <  8:96039/147‑1 (MQ=255)
                                                                                                ttGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGc                                               >  3:168615/1‑148 (MQ=255)
                                                                                                ttGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGc                                               >  7:243733/1‑148 (MQ=255)
                                                                                                      cTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTg                                          >  3:61670/1‑147 (MQ=255)
                                                                                                       tGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACg                                                                             <  2:172438/111‑1 (MQ=255)
                                                                                                       tGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACg                                                                             >  1:172438/1‑111 (MQ=255)
                                                                                                         gACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTgcag                                       <  2:146949/147‑1 (MQ=255)
                                                                                                                     ggTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCgg                         <  4:61670/149‑1 (MQ=255)
                                                                                                                     ggTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCgg                         >  1:156822/1‑149 (MQ=255)
                                                                                                                       tttGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGgca                       <  8:283729/149‑1 (MQ=255)
                                                                                                                         tGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGgcag                      <  4:168615/148‑1 (MQ=255)
                                                                                                                            cGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAat                  <  1:43014/149‑1 (MQ=255)
                                                                                                                                cttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGgca                       >  4:239265/5‑140 (MQ=255)
                                                                                                                                cttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGgca                       <  3:239265/136‑1 (MQ=255)
                                                                                                                                cttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAatga                >  2:94328/5‑147 (MQ=255)
                                                                                                                                cttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAatga                <  1:94328/143‑1 (MQ=255)
                                                                                                                                 cgttcgttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATAc                                                                                                 >  6:74074/1‑65 (MQ=255)
                                                                                                                                 cgttcgttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATAc                                                                                                 <  5:74074/65‑1 (MQ=255)
                                                                                                                                  gttcgttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGTCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGcc            <  1:201135/149‑1 (MQ=255)
                                                                                                                                   ttcgttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCAATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCCg           >  8:194058/1‑149 (MQ=255)
                                                                                                                                      gttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTgc                                         <  4:179144/116‑1 (MQ=255)
                                                                                                                                       ttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCCGTcgc       >  3:169817/1‑149 (MQ=255)
                                                                                                                                             cTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCCGTCGCGGGtt  <  2:156822/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCCGTCGCGGGTT  >  NC_002947/1628854‑1629142

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