Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,998,928 T→C 7.4% R20R (CGA→CGG)  PP_2622 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,998,9280TC 7.4% 151.7 / 2.3 54R20R (CGA→CGG) PP_2622hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (2/2);  total (27/27)
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.99e-01

CAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACGGCGAGCGAAACCTGGTAC  >  NC_002947/2998783‑2999074
                                                                                                                                                 |                                                                                                                                                  
caagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCgg                                                                                                                                                 <  7:269984/149‑1 (MQ=255)
caagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGCCTGCTCTGGGAGCACCCCTCGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTAGAGGGAAGCGCTCTCCCGTCgg                                                                                                                                                 >  7:230945/1‑149 (MQ=255)
   gcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAg                                                                                                                                              >  1:374071/1‑149 (MQ=255)
   gcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGCCGagag                                                                                                                                              >  4:360330/1‑149 (MQ=255)
    caagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTATCCCGTCGGGAgg                                                                                                                                             >  6:230096/1‑149 (MQ=255)
    aaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCATGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGTGAgg                                                                                                                                             >  6:224689/2‑149 (MQ=255)
     aagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGc                                                                                                                                            <  4:187728/149‑1 (MQ=255)
       gTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCct                                                                                                                                          >  7:85789/1‑149 (MQ=255)
         cGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCctct                                                                                                                                        >  2:301610/1‑149 (MQ=255)
                 cTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGCCGAGAGGCCTCCcgacgg                                                                                                                                  >  8:178739/1‑147 (MQ=255)
                    aCGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGCCGAGAGGCCTCCCGACgggaga                                                                                                                              <  7:178739/148‑4 (MQ=18)
                    aCGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGCCGAGAGGCCTCCCGACgggag                                                                                                                               <  3:360330/147‑3 (MQ=255)
                       tCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTg                                                                                                                          <  8:303166/149‑1 (MQ=255)
                         aTGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGt                                                                                                                        <  6:75647/149‑1 (MQ=255)
                              cATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTccacc                                                                                                                   <  1:301610/149‑1 (MQ=255)
                                ttGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCAt                                                                                                                 >  3:84871/1‑149 (MQ=255)
                                   aaGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAg                                                                                                              >  1:15409/1‑149 (MQ=255)
                                      tCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCg                                                                                                           >  3:39541/1‑149 (MQ=255)
                                        gcgcAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGAt                                                                                                         <  8:85789/149‑1 (MQ=255)
                                         cgcAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATc                                                                                                        <  8:69097/149‑1 (MQ=255)
                                              gCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTTCACTGGTCCACCATCAGGCGATCaaaaa                                                                                                   <  5:230096/149‑1 (MQ=255)
                                              gCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGGGGCGATTGATGGCTGCGAACTTGTGGGAAGCGCTCTCCCGTCGTGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaaa                                                                                                   <  5:224689/149‑1 (MQ=255)
                                                ccAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGa                                                                                                 >  5:368044/1‑149 (MQ=255)
                                                ccAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGa                                                                                                 >  3:162088/1‑149 (MQ=255)
                                                  aGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGAcc                                                                                               <  4:84871/149‑1 (MQ=255)
                                                      cAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATc                                                                                                                >  6:284417/1‑128 (MQ=255)
                                                      cAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATc                                                                                                                <  5:284417/128‑1 (MQ=255)
                                                          gctgctCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTg                                                                                       >  1:215160/1‑149 (MQ=255)
                                                          gctgctCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTg                                                                                       <  4:148561/149‑1 (MQ=255)
                                                                   ggAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAg                                                                              >  2:104604/1‑149 (MQ=255)
                                                                         ccccTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGccc                                                                        <  1:104604/149‑1 (MQ=255)
                                                                          cccTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGccc                                                                        >  1:11171/1‑148 (MQ=255)
                                                                          cccTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGccc                                                                        >  5:141008/1‑148 (MQ=255)
                                                                           ccTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGa                                                                      <  4:39541/149‑1 (MQ=255)
                                                                            cTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGAc                                                                     >  2:231279/1‑149 (MQ=255)
                                                                               gCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAg                                                                  <  6:71008/149‑1 (MQ=255)
                                                                                ggggAATTCAGCAGCGTTTCGTGGGGGGGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTGTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAgg                                                                 <  2:374071/148‑1 (MQ=255)
                                                                                     aTTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATAGAGCCGCCCGACCAGGGAAAt                                                            >  1:44981/1‑149 (MQ=255)
                                                                                           cagcGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGaa                                                      <  2:15409/149‑1 (MQ=255)
                                                                                                        tGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATg                                         <  6:141008/149‑1 (MQ=255)
                                                                                                                 aTGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGaa                                 <  7:30345/148‑1 (MQ=255)
                                                                                                                  tGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGaa                                 <  3:180523/147‑1 (MQ=255)
                                                                                                                  tGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGaa                                 >  4:180523/1‑147 (MQ=255)
                                                                                                                    gcggcgAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGc                             <  4:162088/149‑1 (MQ=255)
                                                                                                                             ttGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACg                    >  4:380275/1‑149 (MQ=255)
                                                                                                                              tgtgGGATGCGCTCTCCCGTCGTGAGGCATCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACg                    <  6:49438/148‑1 (MQ=255)
                                                                                                                                  ggAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATg                                         <  3:282910/123‑1 (MQ=255)
                                                                                                                                  ggAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATg                                         >  4:282910/1‑123 (MQ=255)
                                                                                                                                              ccGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGCCGACTCATTGCGTATCGAGCCGCCCGCCCAGGGAAATCTCGAAGTTCGCGCCCATGTAGGTGAAGTGCGGGCGTTCGGCGAGCGAAACCTggga   >  4:95627/1‑147 (MQ=255)
                                                                                                                                              ccGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACGGCGAGCGAAACCTgg     <  1:231279/147‑1 (MQ=255)
                                                                                                                                              ccGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACGGCGAGCGAAACCTGGTa   >  2:125400/1‑149 (MQ=255)
                                                                                                                                                gTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCTCGCCCATGTACTTGAAGTGCGGGCGTACGgcgagcga            <  3:52793/138‑1 (MQ=255)
                                                                                                                                                gTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACGgcgagcga            >  4:52793/1‑138 (MQ=255)
                                                                                                                                                gTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACGGCGAGCGAAACCTGGTAc  >  3:213474/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTGGCGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCCGACCAGGGAAATCTCGAAGTTCGCGCCCATGTACTTGAAGTGCGGGCGTACGGCGAGCGAAACCTGGTAC  >  NC_002947/2998783‑2999074

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