Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,113 T→A 6.9% L59M (TTG→ATG) ‡ PP_2445 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,793,1130TA 6.9% 163.4 / 2.5 58L59M (TTG→ATG) ‡PP_2445membrane protein
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base A (2/2);  total (28/30)
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.60e-01

CCCTACCGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGC  >  NC_002947/2792968‑2793252
                                                                                                                                                 |                                                                                                                                           
cccTACCGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTgg                                                                                                                                          >  2:244781/1‑149 (MQ=255)
  cTACCGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGt                                                                                                                                         >  7:351469/1‑148 (MQ=255)
     ccGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCCCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATg                                                                                                                                     <  5:49167/149‑1 (MQ=255)
      cgacTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGc                                                                                                                                    <  5:204451/146‑1 (MQ=255)
      cGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGc                                                                                                                                    <  6:66500/149‑1 (MQ=255)
       gCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGc                                                                                                                                    >  2:97465/1‑148 (MQ=255)
         cTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGAc                                                                                                                                 >  7:42870/1‑149 (MQ=255)
            ggTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACACCCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGGCATGCGACTg                                                                                                                               <  2:29553/148‑1 (MQ=255)
            ggTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACt                                                                                                                                <  5:171085/147‑1 (MQ=255)
                      gCTACGCTGGTGGCAATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTc                                                                                                                      <  7:97209/147‑1 (MQ=255)
                               gtggCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGc                                                                                                           <  7:334460/149‑1 (MQ=255)
                                   ccATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTg                                                                                                        >  6:379181/1‑148 (MQ=255)
                                    cATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGa                                                                                                      <  5:341209/149‑1 (MQ=255)
                                    cATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACAAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGa                                                                                                      >  4:40672/1‑149 (MQ=255)
                                     aTGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCAt                                                                                                                                      >  6:165210/1‑116 (MQ=255)
                                     aTGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCAt                                                                                                                                      <  5:165210/116‑1 (MQ=255)
                                     aTGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGAt                                                                                                     >  4:259245/1‑149 (MQ=255)
                                       ggAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGGCATGAGACTGGCCTTGCTCAGCACCGTGGCGTGGATCg                                                                                                   >  8:356928/1‑149 (MQ=255)
                                         aaGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTc                                                                                                 <  5:174661/149‑1 (MQ=255)
                                                 cTGCGCATCGACAACCTCATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgac                                                                                         <  4:76380/149‑1 (MQ=255)
                                                  tGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTgg                                                                                                                              <  4:260649/111‑1 (MQ=255)
                                                  tGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTgg                                                                                                                              >  3:260649/1‑111 (MQ=255)
                                                  tGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgacc                                                                                        >  2:93499/1‑149 (MQ=255)
                                                    ggaaTCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgacc                                                                                        <  2:115617/144‑1 (MQ=255)
                                                    ggCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGACCAAGGCCATGctga                                                                                                                                 >  1:38270/1‑103 (MQ=255)
                                                    ggCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGACCAAGGCCATGctga                                                                                                                                 <  2:32530/106‑4 (MQ=255)
                                                    ggCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGACCAAGGCCATGctga                                                                                                                                 >  1:32530/1‑103 (MQ=255)
                                                    ggCATAGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGACCAAGGCCATGctga                                                                                                                                 <  2:38270/106‑4 (MQ=255)
                                                          gACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGgt                                                                                >  1:357595/1‑149 (MQ=255)
                                                            cAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGgtgt                                                                              <  2:108887/149‑1 (MQ=255)
                                                                 tGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGAATGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCg                                                                           >  5:378089/1‑147 (MQ=255)
                                                                   aTCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGaa                                                                         <  2:357595/147‑1 (MQ=255)
                                                                        cATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTc                                                                   <  5:151453/148‑1 (MQ=255)
                                                                        cATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCg                                                                  >  7:77090/1‑149 (MQ=255)
                                                                          tctcGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCgg                                                                 >  1:71073/1‑148 (MQ=255)
                                                                            tcGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTGTCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCaa                                                              >  2:30708/1‑149 (MQ=255)
                                                                                    gACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCCCCGTGGCCTGGATCGTCCCGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCGTtct                                                      >  6:214826/1‑149 (MQ=255)
                                                                                     aCCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTtct                                                      <  3:21668/148‑1 (MQ=255)
                                                                                      ccAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCg                                                    >  5:140714/1‑149 (MQ=255)
                                                                                       caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTc                                                                   <  3:306675/133‑1 (MQ=255)
                                                                                       caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTc                                                                   >  4:306675/1‑133 (MQ=255)
                                                                                       caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                   >  7:22251/1‑149 (MQ=255)
                                                                                       caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                   <  1:97465/149‑1 (MQ=255)
                                                                                       caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGAATCGGTATCAGCATGGACTTCGTCATGAGACTGGCCTTGCTCAGCACCGTGGCCTGGAAAGTCCACATTACCGACCCGGTGTACCAACTGATCGGCAACGCCTGCTCGt                                                   >  4:233193/1‑149 (MQ=255)
                                                                                        aacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACAGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTg                                                  >  4:127512/1‑149 (MQ=255)
                                                                                          caaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGa                                                <  5:379181/149‑1 (MQ=255)
                                                                                           aaaTTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTgg                                                 >  1:220572/1‑147 (MQ=255)
                                                                                                gCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGAt                                           <  3:259245/148‑1 (MQ=255)
                                                                                                 ccGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGAtg                                          <  1:244781/148‑1 (MQ=255)
                                                                                                        aGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCt                                   <  6:378089/148‑1 (MQ=255)
                                                                                                         gTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGa                                 <  1:251005/149‑1 (MQ=255)
                                                                                                             cGCTCCAAGGCGCGGCGTCTCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGc                             <  1:30708/149‑1 (MQ=255)
                                                                                                             cGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGc                             >  3:137348/1‑149 (MQ=255)
                                                                                                               cTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCg                           <  8:77090/149‑1 (MQ=255)
                                                                                                                     ggCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTc                     >  7:129937/1‑149 (MQ=255)
                                                                                                                     ggCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTc                     <  8:42870/149‑1 (MQ=255)
                                                                                                                              tATCGGTATTAGCATGTCCTTGGTGATGCGACTGGCGTTGCTCAGCACCGTGGCCTGGATCGTACAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGTTCGTGTTTGCCGGTGGTCTGTTCctgc            <  3:233193/149‑1 (MQ=255)
                                                                                                                                        aGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGc  >  5:157802/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CCCTACCGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGC  >  NC_002947/2792968‑2793252

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