Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,035,850 T→C 7.8% A55A (GCA→GCG)  PP_5477 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,035,8500TC 7.8% 142.8 / 2.5 51A55A (GCA→GCG) PP_5477hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base C (2/2);  total (23/28)
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.52e-01

TACGCAATGTGGCACTCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGACTACACATTTGCAAAAAAAGC  >  NC_002947/2035717‑2035997
                                                                                                                                     |                                                                                                                                                    
tACGCAATGTGGCACTCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACt                                                                                                                                       >  7:310201/1‑149 (MQ=255)
     aaTGTGGCACTCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCcac                                                                                                                                  <  5:114994/149‑1 (MQ=255)
            cACTCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCAGTTCTGCTGCGCGGACCGCTACTGCCACCTT‑AGGc                                                                                                                          >  8:307077/1‑149 (MQ=255)
             aCTCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCCGCGCTGAACGCTACTGCGACATT‑ACGCt                                                                                                                         >  2:44656/1‑149 (MQ=255)
             aCTCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCCGCGCTGAACGCTACTGCCACATT‑ACGCt                                                                                                                         >  2:50011/1‑149 (MQ=255)
               tCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCg                                                                                                                       >  7:290996/1‑149 (MQ=255)
                 aaGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGt                                                                                                                      >  7:126789/1‑148 (MQ=255)
                          tctTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTccac                                                                                                            >  7:319053/1‑149 (MQ=255)
                          tctTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATA‑ACGCTCGTAACTCTccac                                                                                                            >  7:283333/1‑149 (MQ=255)
                                  aCGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATcgc                                                                                                    <  8:290996/149‑1 (MQ=255)
                                  aCGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATcgc                                                                                                    <  3:336541/149‑1 (MQ=255)
                                          tatCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGAt                                                                                            >  7:257384/1‑149 (MQ=255)
                                          tatCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCCGCGCTGAACGCTACTGCCACATTCACGCTCGTAACTCTCCACCAGATCGCGAGTCGa                                                                                             <  1:44656/149‑1 (MQ=255)
                                          tatCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCCGCGCTGAACGCTACTGCCACATTAACGCTCGTAACTCTCCACCAGATCGCGAGTCGa                                                                                             <  1:50011/149‑1 (MQ=255)
                                           atCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATcgc                                                                                                    <  8:182998/140‑1 (MQ=255)
                                           atCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATcgc                                                                                                    >  7:182998/1‑140 (MQ=255)
                                              ggTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCg                                                                                        >  6:218776/1‑149 (MQ=255)
                                              ggTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCg                                                                                        <  2:25738/149‑1 (MQ=255)
                                               gTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCgt                                                                                       <  8:257384/149‑1 (MQ=255)
                                                  aaGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTcca                                                                                                             <  7:286049/124‑1 (MQ=255)
                                                  aaGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTcca                                                                                                             >  8:286049/1‑124 (MQ=255)
                                                   aGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTcc                                                                                   <  8:65917/149‑1 (MQ=255)
                                                        ttCTGAAATTGGTGCTCTTTCACAATGATCAGGGCCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGt                                                                              <  3:160174/149‑1 (MQ=255)
                                                              aaTTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGc                                                                        >  2:155096/1‑149 (MQ=255)
                                                                   gTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCa                                                                    <  5:218776/148‑1 (MQ=255)
                                                                    tGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGt                                                                  <  3:648/149‑1 (MQ=255)
                                                                     gCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGt                                                                  >  7:226087/1‑148 (MQ=255)
                                                                      ctctTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGa                                                                <  2:138398/149‑1 (MQ=255)
                                                                         tttCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCa                                                                    >  1:351486/1‑142 (MQ=255)
                                                                         tttCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCa                                                                    <  2:351486/142‑1 (MQ=255)
                                                                            cacaATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCcgc                                                          >  2:348523/1‑149 (MQ=255)
                                                                               aaTGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCgcggc                                                       >  5:353852/1‑149 (MQ=255)
                                                                                           cACGCAGTTTGCACCTTTTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCt                                           <  4:144721/149‑1 (MQ=255)
                                                                                            aCGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCtt                                          <  1:348523/149‑1 (MQ=255)
                                                                                             cGCAGTTTGCACCTTGTGAAAGAAGGTCTTGCATTTCTCCTGCGCTGATCGCTACTGCGACATT‑ACGCGCGTAACTCTCCCCCTGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGGCCGCGGCGATTTGCCATCTTg                                         <  5:185535/149‑1 (MQ=255)
                                                                                              gCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTg                                         >  4:340829/1‑148 (MQ=255)
                                                                                                aGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGAc                                      >  2:353472/1‑149 (MQ=255)
                                                                                                  tttGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTg                                    >  5:33899/1‑149 (MQ=255)
                                                                                                      cACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCCCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGtt                                <  3:269198/149‑1 (MQ=255)
                                                                                                       aCCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTa                               <  8:126789/149‑1 (MQ=255)
                                                                                                             tgAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACAt                         <  8:310201/149‑1 (MQ=255)
                                                                                                              gAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATcgc                                                                                                    >  1:178428/1‑73 (MQ=255)
                                                                                                              gAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATcgc                                                                                                    <  2:178428/73‑1 (MQ=255)
                                                                                                                aagaagGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGAc                      <  6:33899/149‑1 (MQ=255)
                                                                                                                aagaagGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGAc                      <  1:353472/149‑1 (MQ=255)
                                                                                                                aagaagGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGAc                      <  6:353852/149‑1 (MQ=255)
                                                                                                                    agGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGACTaca                  <  5:22072/149‑1 (MQ=255)
                                                                                                                          ttCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGAttt                                                  >  2:80322/1‑111 (MQ=255)
                                                                                                                          ttCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGAttt                                                  <  1:80322/111‑1 (MQ=255)
                                                                                                                           tCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCTAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGACTACACAttt             <  8:283333/147‑1 (MQ=255)
                                                                                                                              tttCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACTAGATCGCGAGTCGATGGGGTGGGCATTTTAAAAGCCCCAGTGATTCCGCGGCGGTTTTCCATATTGGGCTTGGTTATAAAATGACTAAACATTTTCaa         >  1:275578/1‑148 (MQ=255)
                                                                                                                                     tgcGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGACTACACATTTGCAAAAAAAGc  >  4:89460/1‑148 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
TACGCAATGTGGCACTCAAGACGACCTCTTGCTGACGGTTGATATCGGTCAAGAAATTCTGAAATTGGTGCTCTTTCACAATGATCAGGACCACGCAGTTTGCACCTTGTGAAAGAAGGTCTTTCATTTCTGCTGCGCTGAACGCTACTGCCACATT‑ACGCTCGTAACTCTCCACCAGATCGCGAGTCGATGGCGTGTCCATAGTCAAAGCCCCAGTGAGTCCGCGGCGATTTGCCATCTTGGACTGGGTTAGAACATGACTACACATTTGCAAAAAAAGC  >  NC_002947/2035717‑2035997

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