Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,568,693 T→A 7.0% Q285H (CAA→CAT PP_2252 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,568,6930TA 7.0% 160.6 / 3.1 57Q285H (CAA→CATPP_2252hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base A (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

CACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTT  >  NC_002947/2568552‑2568839
                                                                                                                                             |                                                                                                                                                   
cACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCCCGTGCAGGATCACCCGCAGTTGTTCGa                                                                                                                                              >  7:325653/1‑149 (MQ=255)
   cGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGATGGTCGAACa                                                                                                                                           <  2:15434/149‑1 (MQ=255)
          ggcgCACTGGCCAGCATGCGGTCCTCACACTGCCGCGACAGTTAATAGACCCCCGCCAGGCCCCCGCTGTGGCGGTTCTGCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCgg                                                                                                                                    <  1:109320/149‑1 (MQ=255)
               aCTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGggtggt                                                                                                                               >  7:285726/1‑149 (MQ=255)
                   gCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                          <  5:7470/149‑1 (MQ=255)
                        cAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCAt                                                                                                                     >  8:184179/1‑149 (MQ=255)
                         aaGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATg                                                                                                                    <  7:299388/149‑1 (MQ=255)
                          aGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                   <  4:105006/149‑1 (MQ=255)
                           gCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATggg                                                                                                                  <  8:285726/149‑1 (MQ=255)
                              gTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACgg                                                                                                                         <  5:350972/139‑1 (MQ=255)
                              gTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACgg                                                                                                                         >  6:350972/1‑139 (MQ=255)
                                 cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcatc                                                                                                             <  7:184179/148‑1 (MQ=255)
                                  tCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                           >  7:109890/1‑149 (MQ=255)
                                     aCCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcccc                                                                                                        >  6:240809/1‑149 (MQ=255)
                                                 aGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACTATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGc                                                                                            >  2:112741/1‑149 (MQ=255)
                                                    tAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACGATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                          >  7:186304/1‑148 (MQ=255)
                                                         gACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATc                                                                                     <  8:109890/148‑1 (MQ=255)
                                                             cccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATg                                                                                                                    >  4:347873/1‑113 (MQ=255)
                                                             cccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGATCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAat                                                                                <  5:240809/149‑1 (MQ=255)
                                                                       ccccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCt                                                                      <  8:186304/149‑1 (MQ=255)
                                                                         ccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCgg                                                                                                <  3:64264/121‑1 (MQ=255)
                                                                         ccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCgg                                                                                                >  4:64264/1‑121 (MQ=255)
                                                                         ccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGc                                                                    >  3:28975/1‑149 (MQ=255)
                                                                         ccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGc                                                                    >  1:139006/1‑149 (MQ=255)
                                                                                       ctcCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGATGGTCGaac                                                                                                                                            >  5:269343/1‑62 (MQ=39)
                                                                                       ctcCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGATGGTCGaac                                                                                                                                            <  6:269343/64‑3 (MQ=39)
                                                                                        tcCATCAACTCCTGCCATTGCTGGTTGAACGCCACGTTCAGGATCACCCGCAGTTGTTGGCCCATCGGCGGGTGGTGACGGTCATGGGTCATGACCCCGTAACCGGCTGCGGCAATCGCAATAATCGCCCCTGCCCCCGTGCCCACTGc                                                     >  2:140906/1‑149 (MQ=255)
                                                                                         ccATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGc                                                            <  7:317132/141‑1 (MQ=255)
                                                                                         ccATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGc                                                            >  8:317132/1‑141 (MQ=255)
                                                                                          cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                          <  4:213660/110‑1 (MQ=255)
                                                                                          cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                          >  3:213660/1‑110 (MQ=255)
                                                                                           aTCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGa                                                  >  3:219933/1‑149 (MQ=255)
                                                                                            tCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGAt                                                 <  8:325653/149‑1 (MQ=255)
                                                                                                 tCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCgg                                                                                                                                    >  5:158657/1‑62 (MQ=255)
                                                                                                 tCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCgg                                                                                                                                    <  6:158657/62‑1 (MQ=255)
                                                                                                 tCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCgg                                                                                                                                    <  4:47364/62‑1 (MQ=255)
                                                                                                 tCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCgg                                                                                                                                    >  3:47364/1‑62 (MQ=255)
                                                                                                    tGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACgg                                          >  1:190724/1‑148 (MQ=255)
                                                                                                     gCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCa                                        >  3:145627/1‑149 (MQ=255)
                                                                                                      ccATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGc                                                                          >  8:235694/1‑114 (MQ=255)
                                                                                                       cATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACt                                      >  6:13021/1‑149 (MQ=255)
                                                                                                           ccTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGTGGCATCGGCAGTggg                                   >  1:238888/1‑147 (MQ=255)
                                                                                                           ccTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGc                                   <  4:219933/148‑1 (MQ=255)
                                                                                                           ccTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGc                                   <  4:145627/148‑1 (MQ=255)
                                                                                                            cTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGTGGCATCGGCAGTgggc                                  <  2:238888/148‑3 (MQ=255)
                                                                                                            cTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGAt                                                 >  7:127280/1‑133 (MQ=255)
                                                                                                            cTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGAt                                                 <  8:127280/133‑1 (MQ=255)
                                                                                                            cTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCg                                  >  3:61862/1‑148 (MQ=255)
                                                                                                             tGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACAGGCTGCGGAAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGAt                                >  5:171304/1‑149 (MQ=255)
                                                                                                                     cGCCACGTTCATGATCACCCGCAGTTGTTCGACTATCTGCGGGTGGTCACCCGTCATGGGTCATCACCCCGAAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGt                         <  1:112741/149‑1 (MQ=255)
                                                                                                                      gCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTgg                       <  1:140906/149‑1 (MQ=255)
                                                                                                                                 gATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTgg            >  6:29843/1‑149 (MQ=255)
                                                                                                                                     aCCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCAGCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcgc         <  1:24135/148‑1 (MQ=255)
                                                                                                                                       ccGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCaa      <  4:28975/149‑1 (MQ=255)
                                                                                                                                       ccGCAGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGcc        >  7:331043/1‑147 (MQ=255)
                                                                                                                                           aGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCtt  >  3:91535/1‑149 (MQ=255)
                                                                                                                                           aGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCtt  >  3:96406/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTT  >  NC_002947/2568552‑2568839

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: