Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,922,427 G→C 11.1% T117S (ACC→AGC) ‡ PP_4331 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,922,4270GC11.1% 153.5 / 8.8 54T117S (ACC→AGC) ‡PP_4331hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/21);  new base C (3/3);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CAGCGGGTACTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGCCACCGCCGCACGCAGCTCAT  >  NC_002947/4922280‑4922574
                                                                                                                                                    |                                                                                                                                                   
cAGCGGGTACTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGTGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGt                                                                                                                                                    >  5:261886/1‑149 (MQ=255)
 aGCGGGTACTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGt                                                                                                                                                    >  1:82508/1‑148 (MQ=255)
 aGCGGGTACTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTc                                                                                                                                                   <  6:125186/149‑1 (MQ=255)
         cTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCg                                                                                                                                             >  6:70382/1‑147 (MQ=255)
         cTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGACCCAGGCCgag                                                                                                                                           <  5:186007/149‑3 (MQ=255)
         cTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGACCCAGGCCgag                                                                                                                                           <  5:106377/149‑3 (MQ=255)
                aGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTgg                                                                                                                                    <  1:153551/149‑1 (MQ=255)
                       cGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACt                                                                                                                             <  8:13382/149‑1 (MQ=255)
                        gACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACt                                                                                                                             >  5:176694/1‑148 (MQ=255)
                          cacCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCAt                                                                                                                          <  5:87338/149‑1 (MQ=255)
                             cTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGa                                                                                                                       <  3:119771/149‑1 (MQ=255)
                              tttCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGACCCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGAACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGa                                                                                                                       >  5:84436/1‑148 (MQ=255)
                              tttCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGACCCAGGCCgag                                                                                                                                           >  5:52035/1‑126 (MQ=255)
                              tttCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGACCCAGGCCgag                                                                                                                                           <  6:52035/128‑3 (MQ=255)
                                  aGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTc                                                                                                                  <  2:88785/149‑1 (MQ=255)
                                   gATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCATGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCg                                                                                                                 >  6:157561/1‑149 (MQ=255)
                                        aGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGcccc                                                                                                             >  2:144621/1‑148 (MQ=255)
                                        aGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGcccc                                                                                                             >  4:39779/1‑148 (MQ=255)
                                                      cgccgcTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCtt                                                                                              >  6:75743/1‑149 (MQ=255)
                                                       gccgcTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTg                                                                                             <  6:176694/149‑1 (MQ=255)
                                                          gcTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGgcc                                                                                          >  8:219040/1‑149 (MQ=255)
                                                                tttCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTgg                                                                                    <  7:219040/149‑1 (MQ=255)
                                                                 ttCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTgg                                                                                    <  7:295375/148‑1 (MQ=255)
                                                                  tCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCt                                                                                  <  1:117057/149‑1 (MQ=255)
                                                                    ccAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCAGTGATCTTCAACTCCCCCACGCCCTGGGTCAGGCAGCAGAGGTGGGTCACCTCATCCATAGACGCCCCCAGCCCCCCCCCCTTGCCCGCCTgcgg                                                                                  >  1:85464/1‑145 (MQ=255)
                                                                                ccTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGc                                                                    >  2:305048/1‑149 (MQ=255)
                                                                                  tgctgcCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTg                                                                  >  3:235796/1‑149 (MQ=255)
                                                                                     tgcCGAGGGTACCCCGTCAAGGTGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTgggg                                                               >  7:49777/1‑149 (MQ=255)
                                                                                       ccGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTc                                                             >  4:215321/1‑149 (MQ=255)
                                                                                            ggTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCgct                                                        >  7:108922/1‑149 (MQ=255)
                                                                                                ccccGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTc                                                    >  4:300116/1‑149 (MQ=255)
                                                                                                 cccGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTc                                                    >  5:298485/1‑148 (MQ=255)
                                                                                                    gTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAgg                                                <  8:49777/149‑1 (MQ=255)
                                                                                                    gTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAgg                                                <  1:144621/149‑1 (MQ=255)
                                                                                                        aGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGtt                                            <  3:215321/149‑1 (MQ=255)
                                                                                                         ggAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTc                                           <  1:305048/149‑1 (MQ=255)
                                                                                                           aGCCCGGTACACCTTGGTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTTGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCa                                         <  2:82508/149‑1 (MQ=255)
                                                                                                             cctggcccacCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGc                                       >  2:9122/8‑149 (MQ=255)
                                                                                                             cccGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGc                                       <  5:75743/149‑1 (MQ=255)
                                                                                                             cccGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGc                                       <  2:90899/149‑1 (MQ=255)
                                                                                                                  tACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTc                                  <  3:300116/149‑1 (MQ=255)
                                                                                                                     acCTTGCTCACCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCgg                               >  4:40875/1‑149 (MQ=255)
                                                                                                                      ccTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCgg                               >  2:225568/1‑148 (MQ=255)
                                                                                                                      ccTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGc                              <  4:235796/149‑1 (MQ=255)
                                                                                                                               atcaACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGTCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTg                        <  6:261886/147‑1 (MQ=255)
                                                                                                                               atcaACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGc                      >  7:158616/1‑149 (MQ=255)
                                                                                                                                tcaACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGc                      >  4:286289/1‑148 (MQ=255)
                                                                                                                                     tCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCgg                                    >  5:173337/1‑129 (MQ=255)
                                                                                                                                     tCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCgg                                    <  6:173337/129‑1 (MQ=255)
                                                                                                                                         ctgcGGCCTGACCCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGCCACCGCcgca            >  4:128576/4‑149 (MQ=255)
                                                                                                                                         ctgcGGCCTGACCCAGGCAGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCGGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGAGTGGCCACCGCcgca            >  8:269540/4‑149 (MQ=255)
                                                                                                                                         ccTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCGTCGATACTAGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTAGGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGCCCCCGCcgca            >  2:266451/1‑149 (MQ=255)
                                                                                                                                               ccTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCCTCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGCCACCGCCGCACGCAGc      >  2:263304/1‑149 (MQ=255)
                                                                                                                                               ccTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGCCACCGCCGCACGCAg       >  4:214891/1‑148 (MQ=255)
                                                                                                                                                   ggTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGCCACCGCCGCACGCAGCTCat  <  5:157561/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCGGGTACTGAACGAGCCTGGCGACACCTTTCAGATCGAGCGGAGCATGGGGCGCCGCTCGATTTCCCAGGCGACGAACCTGCTGCCGAGGGTACCCCGTCAAGGAGCCCGGTACACCTTGCTCATCAACTCCGCCTCGGCCTGGGTCAGGCCGCAGCTCTGGGTCAACTCATCGATACTCGCCCCCATCCCCACCAGCTTGGCCGCCTGGCTGAAAGTTACGTTGCTGGGGTCGCGCTGCTCGAGGCGTTCCAGCCGGTCCGGCAGGGTGGCCACCGCCGCACGCAGCTCAT  >  NC_002947/4922280‑4922574

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: