Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 122,168 C→T 8.7% A705V (GCG→GTG)  katE → hydroperoxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947122,1680CT 8.7% 125.6 / 2.6 46A705V (GCG→GTG) katEhydroperoxidase
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base T (2/2);  total (21/25)
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.16e-01

CTTGAAGGCGCTGGTGCTGACCAAGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGC  >  NC_002947/122022‑122314
                                                                                                                                                  |                                                                                                                                                   
cTTGAAGGCGCTGGTGCTGACCAAGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAggcggcgg                                                                                                                                                   >  4:57378/1‑149 (MQ=255)
 ttGAAGGCGCTGGTGCTGACCAAGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAggcggcgg‑c                                                                                                                                                 >  8:64077/1‑149 (MQ=255)
     aGGCGCTGGTGCTGACCAAGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGa                                                                                                                                              <  7:64077/148‑1 (MQ=255)
      ggCGCTGGTGCTGTCCAAGGATTTGGCGAAAGGGTTGGGGCTAAAGGAGCACAAGGGGCTGTTGCGCGGGGATGATCAGAAGTCAGTGGTTGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAg                                                                                                                                            <  5:181350/149‑1 (MQ=255)
             gTGCTGACCAAGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGAGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTcc                                                                                                                                     >  3:18544/1‑149 (MQ=255)
                     cAAGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGAc                                                                                                                             <  8:196843/149‑1 (MQ=255)
                       aGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGAc                                                                                                                             >  3:161619/1‑147 (MQ=255)
                                cGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTc                                                                                                                  <  3:57378/149‑1 (MQ=255)
                                 gAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTcg                                                                                                                 >  8:12895/1‑149 (MQ=255)
                                   aaGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTcgcg                                                                                                               <  4:18544/149‑1 (MQ=255)
                                    aGCGTTGGGGCTAAAGGAGGACATGGGGCTGTTGCTCGGGGATGTGCAGAAGGCAGGGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTcgcg                                                                                                               <  3:68273/148‑1 (MQ=255)
                                     gCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCggg                                                                                                             <  5:78614/149‑1 (MQ=255)
                                        ttGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGAAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCATGCCTGACAAGGCATTTTCCCGGGTa                                                                                                           >  3:126892/1‑148 (MQ=255)
                                             gCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGa                                                                                                                              >  5:22248/1‑124 (MQ=255)
                                             gCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGa                                                                                                                              <  6:22248/124‑1 (MQ=255)
                                                aaaGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGTCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAAccccc                                                                                                    <  1:74742/148‑1 (MQ=255)
                                                aaaGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGGAGAAGGCCACCGGGTGTGGGAGGGGGAGGCGGCGG‑CGGAAGCTGTTCCTGCCGGGCAAGGCGTTTTCGCGGGTAAGccccg                                                                                                    >  2:74742/1‑146 (MQ=255)
                                                   ggaggaCAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCcctcct                                                                                                >  7:1468/1‑148 (MQ=255)
                                                      ggaCTAGGGGCTGTTGCTCGGGGATGAGCAGAAGTCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCTGTAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCTCGGGTAAACCCCCTCCTTATg                                                                                            <  1:147173/149‑1 (MQ=255)
                                                         cAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAg                                                                                         <  5:292639/149‑1 (MQ=255)
                                                               gCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCgg                                                                                   >  8:147409/1‑149 (MQ=255)
                                                                cTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCggg                                                                                  <  5:84032/149‑1 (MQ=255)
                                                                   ttGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGt                                                                                 <  8:52043/147‑1 (MQ=255)
                                                                   ttGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGt                                                                                 <  1:57817/147‑1 (MQ=255)
                                                                     gCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTa                                                                              <  8:1468/148‑1 (MQ=255)
                                                                                 gCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGtt                                                                 >  6:33352/1‑149 (MQ=255)
                                                                                       ggCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTAccc                                                                           >  4:162630/1‑133 (MQ=255)
                                                                                       ggCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTAccc                                                                           <  3:162630/133‑1 (MQ=255)
                                                                                       ggCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGGGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTActgc                                                           <  7:147409/149‑1 (MQ=255)
                                                                                          aGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGa                                                        >  8:255607/1‑149 (MQ=255)
                                                                                                tGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCggg                                                  <  7:12895/149‑1 (MQ=255)
                                                                                                   cTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTGGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                >  4:25512/1‑148 (MQ=255)
                                                                                                    tttGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                <  8:200723/147‑1 (MQ=255)
                                                                                                    tttGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                >  7:200723/1‑147 (MQ=255)
                                                                                                    tttGTCAAAGCGGTAGAAGGCCACCGGGAGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                >  1:189931/1‑147 (MQ=255)
                                                                                                           aaGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcaccac                                       <  2:189931/149‑1 (MQ=255)
                                                                                                                 tAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGTGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGc                                 >  6:248255/1‑149 (MQ=255)
                                                                                                                    aaGGGCACCGGTTTTGGTAGCGGCAGGCGGCGG‑CTTAAGCTGTTCGTGCCTTACATGTCCTTTTCGCGTTTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCaccacc                                      <  5:33833/141‑1 (MQ=255)
                                                                                                                    aaGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGATCCTGCCTGACAAGGCCTTTTCGCGGGTAAAACCCCTCCTGATTTAGGACCCGGTATACCCGCGAATCGTTTAAAGCTGATGCGGGACCAGCAACACCtagaacg                               >  6:33833/1‑142 (MQ=255)
                                                                                                                       gCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCag                           <  5:33352/149‑1 (MQ=255)
                                                                                                                         cACCGGGTGTGGGAGCGGGAGGCGGTGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGc                                 <  3:161362/141‑1 (MQ=255)
                                                                                                                         cACCGGGTGTGGGAGCGGGAGGCGGTGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGc                                 >  4:161362/1‑141 (MQ=255)
                                                                                                                         cACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCCCCTGTGCCGGCagc                          >  4:192708/1‑148 (MQ=255)
                                                                                                                         cACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAAACCCCTCCTGAAGTAGGAGCGGCTTTACCCGCGAAGCGTTTTCTGCTGATGCGGGACCAGCACCCCCTGTGCCTGCagcg                         >  4:25825/1‑149 (MQ=255)
                                                                                                                         cACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCCTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCagc                          >  7:122190/1‑148 (MQ=255)
                                                                                                                                       gCGGGAGGCGGTGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGa           <  5:248255/149‑1 (MQ=255)
                                                                                                                                         gggAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtg         <  6:331397/149‑1 (MQ=255)
                                                                                                                                         gggAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTACAGATGTAGGAGCGGGTTTACCCGCGCAGCGTTTACTGCGGAAGCGGGACCAGCACCACCTGTTCCGGCAGCGAGCGCTGGATCTGAtg         >  8:128451/1‑149 (MQ=255)
                                                                                                                                                ggcgg‑cTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGGGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgcttctgc  <  4:126892/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTTGAAGGCGCTGGTGCTGACCAAGGAATTGGCGAAAGCGTTGGGGCTAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGG‑CTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGC  >  NC_002947/122022‑122314

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