Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,061,704 C→T 9.3% A211A (GCC→GCT)  PP_3579 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,061,7040CT 9.3% 118.2 / 2.8 43A211A (GCC→GCT) PP_3579membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/16);  new base T (2/2);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

ACTACTTCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTGCGCCCATGACAAGG  >  NC_002947/4061556‑4061851
                                                                                                                                                    |                                                                                                                                                   
actactTCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCgccgcc                                                                                                                                                     <  7:318949/149‑1 (MQ=255)
actactTCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCgccgcc                                                                                                                                                     >  4:360863/1‑149 (MQ=255)
    ctTCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCggg                                                                                                                                                  <  8:170988/148‑1 (MQ=255)
     ttCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCTGGGtt                                                                                                                                                >  7:7630/1‑149 (MQ=255)
     ttCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCTGGGtt                                                                                                                                                >  7:2241/1‑149 (MQ=255)
     ttCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGtt                                                                                                                                                >  6:160719/1‑149 (MQ=255)
          gTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTgg                                                                                                                                           <  7:332564/149‑1 (MQ=255)
             cACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCTGGGTTGCTGGGt                                                                                                                                         <  8:7630/148‑1 (MQ=255)
             cACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCTGGGTTGCTGGGt                                                                                                                                         <  8:2241/148‑1 (MQ=255)
              aCCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGtt                                                                                                                                        >  1:188390/1‑148 (MQ=255)
                  tCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGc                                                                                                                                              >  3:75880/1‑138 (MQ=255)
                  tCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGc                                                                                                                                              <  4:75880/138‑1 (MQ=255)
                     gCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTaagcaagc                                                                                                                                <  5:160719/149‑1 (MQ=255)
                       gTGGCGGTGCAAACCTCAGACGTGGGCGGCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCg                                                                                                                              <  4:306961/149‑1 (MQ=255)
                             gTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTc                                                                                                                        <  3:360863/149‑1 (MQ=255)
                                        aGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTg                                                                                                             >  1:176609/1‑149 (MQ=255)
                                                 cGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTc                                                                                                       <  6:97699/146‑1 (MQ=255)
                                                 cGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTc                                                                                                       >  5:97699/1‑146 (MQ=255)
                                                     gCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTc                                                                                                >  2:338217/1‑149 (MQ=255)
                                                                 ttGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCAc                                                                                    >  5:151379/1‑149 (MQ=255)
                                                                      gcggGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAAt                                                                               <  5:4740/149‑1 (MQ=255)
                                                                        gggTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTc                                                                             <  4:143749/149‑1 (MQ=255)
                                                                           tggtgCTGGCGCATTCGCTGGGGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGaa                                                                          <  8:272479/149‑1 (MQ=255)
                                                                                   gcgcATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCt                                                                  >  7:244740/1‑149 (MQ=255)
                                                                                   gcgcATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCt                                                                  >  6:11960/1‑149 (MQ=255)
                                                                                     gcATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCTCGAAAGTCAACAAATTCGAAGGTTACCTgg                                                                >  8:266972/1‑149 (MQ=255)
                                                                                                 gggCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAACGTCACCAAATTCGAAGGTTACCTGGCAACTTGCGGcc                                                    >  8:352756/1‑149 (MQ=255)
                                                                                                  ggCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGcc                                                    <  4:207616/148‑1 (MQ=255)
                                                                                                  ggCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCg                                                   <  6:151379/149‑1 (MQ=255)
                                                                                                  ggCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCg                                                   >  4:43160/1‑149 (MQ=255)
                                                                                                   gCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCgg                                                  >  3:160640/1‑149 (MQ=255)
                                                                                                    cTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCgg                                                  >  8:188436/1‑148 (MQ=255)
                                                                                                              aaTACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGGTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCGGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGttt                                        <  7:162121/148‑1 (MQ=255)
                                                                                                                      gATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGcag                               >  6:165782/1‑149 (MQ=255)
                                                                                                                        ttCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGcagcc                             >  2:271765/1‑149 (MQ=255)
                                                                                                                           tGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGcagccagc                          >  6:13540/1‑149 (MQ=255)
                                                                                                                             gggTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCg                        >  5:159554/1‑149 (MQ=255)
                                                                                                                              ggTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGt                       >  5:368049/1‑149 (MQ=255)
                                                                                                                                  cACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTcac                    >  4:130123/1‑148 (MQ=255)
                                                                                                                                   aCTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCacac                  >  6:84639/1‑149 (MQ=255)
                                                                                                                                   aCTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCacac                  >  3:263882/1‑149 (MQ=255)
                                                                                                                                       tCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTgcg              <  3:130123/149‑1 (MQ=255)
                                                                                                                                       tCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTgcg              <  6:159554/149‑1 (MQ=255)
                                                                                                                                          aCATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTGCGcc            >  7:216718/1‑148 (MQ=255)
                                                                                                                                            aTCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTGCGCCCAt         <  8:216718/149‑1 (MQ=255)
                                                                                                                                                 cgccgGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGttt                                        >  7:271242/1‑113 (MQ=255)
                                                                                                                                                 cgccgGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGttt                                        <  8:271242/113‑1 (MQ=255)
                                                                                                                                                 cgccgGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTGCGCCCATGACaa    <  1:338217/149‑1 (MQ=255)
                                                                                                                                                 cgccgGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTGCGCCCATGACaa    <  5:165782/149‑1 (MQ=255)
                                                                                                                                                   ccgGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTGCGCCCATGACAAgg  <  2:176609/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTACTTCTCGTTCACCATCAGCGTGGCGGTGCAAACCTCAGACGTGGGCGTCGCCGGGCGTGGGTTGCGGCGGGTGGTGCTGGCGCATTCGCTGGTGGGCTTCGTGTTCAATACGGCGATTCTGGGGTTCACTATCAACATCGCCGCCGGGTTGCTGGGTTAAGCAAGCCGCTCGTCTGAAAAACCTGCGCGTCGAACGTCCGCGAAAGTCACCAAATTCGAAGGTTACCTGGCACCTTGCGGCCGGGAAAACGTTTGCCATGCAGCCAGCCGTCACACTGCGCCCATGACAAGG  >  NC_002947/4061556‑4061851

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