Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,245,722 G→A 7.5% A980V (GCG→GTG)  PP_4623 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,245,7220GA 7.5% 150.3 / 2.0 53A980V (GCG→GTG) PP_4623hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base A (2/2);  total (27/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 = 6.42e-01

CCGCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCTGGACAGGATC  >  NC_002947/5245577‑5245870
                                                                                                                                                 |                                                                                                                                                    
ccGCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGccgccg                                                                                                                                                   >  6:133382/1‑149 (MQ=255)
 cGCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCACGGCGGCGGAACGGTCGCCATGGCTCCACGCGATGGCCCGGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCcgg                                                                                                                                                  <  1:230359/149‑1 (MQ=255)
  gCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCcggc                                                                                                                                                 <  1:141528/149‑1 (MQ=255)
  gCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCcggc                                                                                                                                                 >  6:86123/1‑149 (MQ=255)
  gCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCACcggc                                                                                                                                                 >  3:63062/1‑149 (MQ=255)
  gCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCACcggc                                                                                                                                                 >  3:58771/1‑149 (MQ=255)
     ccAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGa                                                                                                                                              <  8:43007/149‑1 (MQ=255)
       aaGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGaaa                                                                                                                                            <  1:87210/149‑1 (MQ=255)
        aGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAAc                                                                                                                                           <  5:133382/149‑1 (MQ=255)
         gACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACg                                                                                                                                          >  4:95162/1‑149 (MQ=255)
            ccccGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCa                                                                                                                                       >  8:17684/1‑149 (MQ=255)
            ccccGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCa                                                                                                                                       <  4:155740/149‑1 (MQ=255)
            ccccGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCACCGGCGGAAACGGCa                                                                                                                                       <  4:63062/149‑1 (MQ=255)
            ccccGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCACCGGCGGAAACGGCa                                                                                                                                       <  4:58771/149‑1 (MQ=255)
              ccGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACg                                                                                                                                     >  3:143465/1‑149 (MQ=255)
              ccGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACg                                                                                                                                     >  7:41486/1‑149 (MQ=255)
               cGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACgg                                                                                                                                    <  6:93788/149‑1 (MQ=255)
                   gAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGa                                                                                                                                 <  5:86123/148‑1 (MQ=255)
                      gCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGAt                                                                                                                             <  7:288476/149‑1 (MQ=255)
                        cACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAg                                                                                                                                <  1:279053/144‑1 (MQ=255)
                        cACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAg                                                                                                                                >  2:279053/1‑144 (MQ=255)
                        cACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGCCCAAGGCCGCCGGCGGACACCGCACGGTGAGGATGc                                                                                                                           >  4:89751/1‑149 (MQ=255)
                          ccGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCg                                                                                                                         <  8:80092/149‑1 (MQ=255)
                                         ccgccgCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATcccc                                                                                                          >  6:195155/1‑149 (MQ=255)
                                                       aGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCcaccaccac                                                                                            >  6:206619/1‑149 (MQ=255)
                                                                         tCGCCATCGCTCCACGCGATGGCCCCGGTTATGAAGGACCTGGCCTCACAGAACAGCCCTACGGCCATGGCCGCCTGCGGAAACGGAACGGTGAGGATGCCGGCCAGTACGATCCCCCACCCCACCACCACTGACGAAATACGCTCggg                                                                          <  3:184181/149‑2 (MQ=255)
                                                                           gCCAACGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAGACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAg                                                                        >  6:75153/1‑149 (MQ=255)
                                                                            ccATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGc                                                                       >  8:285167/1‑149 (MQ=255)
                                                                             cATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCg                                                                      >  3:265753/1‑149 (MQ=255)
                                                                                      aCGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTg                                                             <  5:62192/149‑1 (MQ=255)
                                                                                                  cGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGCAAACGGCACGGTGAGGATGCCCGCCAGTACGATCCACCACTCCACCACCACGGACGAAATACGATCGGCACCGTCGGCGCAGTCTTCACCgcca                                                 >  3:188108/1‑146 (MQ=255)
                                                                                                       atgaAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAGACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGt                                             <  5:75153/148‑1 (MQ=255)
                                                                                                       atgaAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGtt                                            <  4:143465/149‑1 (MQ=255)
                                                                                                         gaAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCGCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGa                                          >  5:198116/1‑149 (MQ=255)
                                                                                                         gaAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAg                                                      >  3:2818/1‑137 (MQ=255)
                                                                                                         gaAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAg                                                      <  4:2818/137‑1 (MQ=255)
                                                                                                          aaGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGa                                          <  3:95162/148‑1 (MQ=255)
                                                                                                          aaGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGa                                          <  7:285167/148‑1 (MQ=255)
                                                                                                           aGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCAGGCCAGTACGATCCCCCACTCCACCACCTCGGACGAAATACGCTCGTCATCGTCGGCGCTGTCTTCATCGGCTTCGGTGag                                         >  6:266523/1‑147 (MQ=255)
                                                                                                              aCCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGatcatc                                     >  1:192751/1‑149 (MQ=255)
                                                                                                                cTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCCGCGGCTTCGGTGATCATCCg                                   >  4:39742/1‑149 (MQ=255)
                                                                                                                cTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCg                                   >  4:98985/1‑149 (MQ=255)
                                                                                                                   gCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCCCTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCt                                >  6:68428/1‑149 (MQ=255)
                                                                                                                    ccGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCt                                >  2:37543/1‑148 (MQ=255)
                                                                                                                     cGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGa                              <  8:20804/149‑1 (MQ=255)
                                                                                                                        acaGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATAc                           >  1:110605/1‑149 (MQ=255)
                                                                                                                          aGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGt                         <  2:110605/149‑1 (MQ=255)
                                                                                                                            aaCAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGttt                       >  5:143470/1‑149 (MQ=255)
                                                                                                                              cAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTtct                     >  1:198409/1‑149 (MQ=255)
                                                                                                                               aGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTtctc                    <  1:37543/149‑1 (MQ=255)
                                                                                                                                    gACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGt               <  6:143470/149‑1 (MQ=255)
                                                                                                                                       ggCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCt            <  4:265753/149‑1 (MQ=255)
                                                                                                                                          cATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCTGGa         <  6:67979/149‑1 (MQ=255)
                                                                                                                                             ggCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCTGGACAg      >  1:44058/1‑149 (MQ=255)
                                                                                                                                                cgccgGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCTGGACAGGAt   >  7:1312/1‑149 (MQ=255)
                                                                                                                                                 gccgGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCTGGACAGGATc  >  4:182569/1‑149 (MQ=255)
                                                                                                                                                 gccgGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCTGGACAGGATc  >  2:288183/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCGCACCAAGACCCCCGTGGAAGCCACCGTTCCTGCAACCCCCGCCGCGAACGACAGCAAGGCGGCGGAACGGTCGCCATCGCTCCACGCGATGGCCCCGGTGATGAAGGACCTGGCCGCACAGAACAGCCCGACGGCCATGGCCGCCGGCGGAAACGGCACGGTGAGGATGCCGGCCAGTACGATCCCCCACTCCACCACCACGGACGAAATACGCTCGGCAGCGTCGGCGCTGTCTTCAGCGGCTTCGTTGATCATCCGCCTGATACGTTTCTCCAGGTCCTGGACAGGATC  >  NC_002947/5245577‑5245870

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