Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,117 T→C 7.3% V60A (GTC→GCC)  PP_2445 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,793,1170TC 7.3% 151.5 / 2.9 56V60A (GTC→GCC) PP_2445membrane protein
Reads supporting (aligned to +/- strand):  ref base T (23/27);  new base C (2/2);  total (26/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 = 6.06e-01

CTACCGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGC  >  NC_002947/2792970‑2793252
                                                                                                                                                   |                                                                                                                                       
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                                                                                                                              >  3:260649/1‑111 (MQ=255)
                                                tGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTgg                                                                                                                              <  4:260649/111‑1 (MQ=255)
                                                tGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgacc                                                                                        >  2:93499/1‑149 (MQ=255)
                                                  ggaaTCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTgacc                                                                                        <  2:115617/144‑1 (MQ=255)
                                                  ggCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGACCAAGGCCATGctga                                                                                                                                 >  1:32530/1‑103 (MQ=255)
                                                  ggCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGACCAAGGCCATGctga                                                                                                                                 >  1:38270/1‑103 (MQ=255)
                                                  ggCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGACCAAGGCCATGctga                                                                                                                                 <  2:32530/106‑4 (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                                                                   >  4:306675/1‑133 (MQ=255)
                                                                                     caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTc                                                                   <  3:306675/133‑1 (MQ=255)
                                                                                     caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                   <  1:97465/149‑1 (MQ=255)
                                                                                     caacaaATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGt                                                   >  7:22251/1‑149 (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)
                                                                                                                                                   |                                                                                                                                       
CTACCGCCTGGGTTGCCCTGGCTACGCTGGTGGCCATGGAAGTTGTACTGGGCATCGACAACCTGATCTTCATCTCGATCCTGACCAACAAATTGCCGGTGGAGTACCGCTCCAAGGCGCGGCGTATCGGTATCAGCATGGCCTTGGTCATGCGACTGGCCTTGCTCAGCACCGTGGCCTGGATCGTCCAGTTGACCGACCCGGTGTTCGAAGTGTTCGGCAACGCCTTCTCGTGGAAGGATGTGATCCTGATTGCCGGTGGTCTGTTCCTGCTGTGGAAGGC  >  NC_002947/2792970‑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: ‑