Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,628,993 T→G 8.0% M168R (ATG→AGG) ‡ 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,9930TG 8.0% 141.2 / 2.4 48M168R (ATG→AGG) ‡mucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base T (18/26);  new base G (2/2);  total (20/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 = 1.00e+00

TGGGTTGGTACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCCGTCGC  >  NC_002947/1628846‑1629137
                                                                                                                                                   |                                                                                                                                                
tGGGTTGGTACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATg                                                                                                                                                 <  2:298786/149‑1 (MQ=255)
        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)
                                                                                                                                                   |                                                                                                                                                
TGGGTTGGTACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCCGTCGC  >  NC_002947/1628846‑1629137

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