Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,321,315 C→G 7.4% P85A (CCG→GCG)  PP_2921 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,321,3150CG 7.4% 152.9 / 2.4 51P85A (CCG→GCG) PP_2921hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base G (2/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCTGTTTGCCGAAGTCAGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCATGAAACCCTGATGCGCTTTC  >  NC_002947/3321171‑3321462
                                                                                                                                                |                                                                                                                                                   
cACCTGTTTGCCGAAGTCAGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGa                                                                                                                                                 >  2:343272/1‑149 (MQ=255)
cACCTGTTTGCCGAAGTCAGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGa                                                                                                                                                 >  2:218185/1‑149 (MQ=255)
      tttGCCGAAGTCAGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGt                                                                                                                                             <  5:5599/147‑1 (MQ=255)
          ccGAAGTCAGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGc                                                                                                                                        <  6:360946/148‑1 (MQ=255)
               gTCAGCGAAAGCATGTCGCTGGAGTCCTTTTTCTAGGGTATGGCCGTCACCAGTGCTGTGGAGCGCATGGCCGTGGATTTCGTGCTTGATGTGTGCCTGGCCTATGCGTCGTTCAAGGTCAGCTTCGACGCGGATCGTATCTCCCCGGa                                                                                                                                  <  2:286592/149‑1 (MQ=255)
                tCAGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAg                                                                                                                                 >  6:45673/1‑149 (MQ=255)
                  aGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGGGTCGAAGCGGAcct                                                                                                                                              <  4:313361/134‑2 (MQ=255)
                  aGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGGGTCGAAGCGGAcct                                                                                                                                              >  3:313361/1‑133 (MQ=255)
                  aGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAgg                                                                                                                                >  6:149662/1‑148 (MQ=255)
                   gCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCa                                                                                                                              >  3:79750/1‑149 (MQ=255)
                      aaaGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTg                                                                                                                           <  6:346330/149‑1 (MQ=255)
                           aTGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCg                                                                                                                       >  1:349194/1‑148 (MQ=255)
                             gTCGCTGGAGGCCTTTTTCAAGGGTATGGCAGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTACAGGTCCGCTTCGACCCTTACCGTATCGCCCCCGAGGCACTGCTCGAAg                                                                                                                    >  5:162441/1‑149 (MQ=255)
                                         cTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCgcggcg                                                                                                        >  7:340431/1‑149 (MQ=255)
                                         cTTTTTCAAGGGTATGGCCGTCACCCGTCCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCTTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGCGGCACTGCTCGAAGCCGTTCtcggct                                                                                                        >  6:274174/1‑148 (MQ=255)
                                                   ggTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGgcg                                                                                                >  5:23722/1‑147 (MQ=255)
                                                         gCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGt                                                                                        >  3:41901/1‑149 (MQ=255)
                                                          ccGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGAGTCGTTCCAGGTACGCTTCGACCCGGACCGTATCGCCCCCGAGGAACTGCTCGACGCCGTGCGCGGCGCCGAGGCGCAggacgtg                                                                                       >  4:31199/1‑144 (MQ=255)
                                                             tCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGcc                                                                                    <  2:349194/149‑1 (MQ=255)
                                                                    tGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGcgtt                                                                             <  1:343272/149‑1 (MQ=255)
                                                                        gTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGcgttcgtt                                                                         <  8:180497/149‑1 (MQ=255)
                                                                         tGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTg                                                                        <  5:45673/149‑1 (MQ=255)
                                                                               gcATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGTACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGAcg                                                                  <  6:162441/149‑1 (MQ=255)
                                                                               gcATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGAcg                                                                  <  6:23722/149‑1 (MQ=255)
                                                                                     cccTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGatc                                                            <  8:340431/149‑1 (MQ=255)
                                                                                      ccTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGatc                                                            >  8:86144/1‑148 (MQ=255)
                                                                                        tGGATGGCGTGGTTGATGTGTGCCTGGCAATTGGGTCGTTCCAGGTCCGCTTCGACCCGTAGCGTATCGCCACGGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGatcatc                                                         <  5:274174/149‑1 (MQ=255)
                                                                                         ggATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGAAAGTATCGCCCACGAGGCACTGATCGAAGCCGTGCGCGGCGACGAGGAGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGatcatc                                                         >  6:370477/1‑148 (MQ=255)
                                                                                          gATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGAcgcg                                                                <  1:144072/140‑1 (MQ=255)
                                                                                          gATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGAcgcg                                                                >  8:319330/1‑140 (MQ=255)
                                                                                          gATGGCGGGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGAcgcg                                                                <  7:319330/140‑1 (MQ=255)
                                                                                                      gATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATcc                                                  >  8:2956/1‑142 (MQ=255)
                                                                                                      gATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATcc                                                  <  7:2956/142‑1 (MQ=255)
                                                                                                      gATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCt                                           >  3:332906/1‑149 (MQ=255)
                                                                                                       aTGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTg                                          >  4:277129/1‑149 (MQ=255)
                                                                                                          tgtgCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCgcggc                                                                                                         <  1:238331/83‑1 (MQ=255)
                                                                                                          tgtgCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCgcggc                                                                                                         >  2:238331/1‑83 (MQ=255)
                                                                                                            tgCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACaa                                     >  2:87115/1‑149 (MQ=255)
                                                                                                               cTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGa                                  >  1:18234/1‑149 (MQ=255)
                                                                                                                tGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGAc                                 >  8:319487/1‑149 (MQ=255)
                                                                                                                  gCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGAc                                 <  4:41901/147‑1 (MQ=255)
                                                                                                                    cAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTg                                                                        <  7:30313/106‑1 (MQ=255)
                                                                                                                    cAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTg                                                                        >  8:163758/1‑106 (MQ=255)
                                                                                                                    cAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTg                                                                        <  7:163758/106‑1 (MQ=255)
                                                                                                                    cAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTg                                                                        >  8:30313/1‑106 (MQ=255)
                                                                                                                             gtTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCAt                     <  3:277129/148‑1 (MQ=255)
                                                                                                                             gtTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCAt                     >  8:34059/1‑148 (MQ=255)
                                                                                                                              ttCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCATg                    >  2:50314/1‑148 (MQ=255)
                                                                                                                                 cAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCATGaaa                 <  4:332906/148‑1 (MQ=255)
                                                                                                                                 acggTCCGGGTCGAAGCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCATGAAAc                >  4:287889/3‑149 (MQ=255)
                                                                                                                                       cGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCATGAAACCCTGAt          <  7:86144/149‑1 (MQ=255)
                                                                                                                                               cccGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCATGAAACCCTGATGCGCTTTc  <  1:50314/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CACCTGTTTGCCGAAGTCAGCGAAAGCATGTCGCTGGAGGCCTTTTTCAAGGGTATGGCCGTCACCCGTGCGGTGGAGCGCATGGCCCTGGATGGCGTGCTTGATGTGTGCCTGGCCAATGCGTCGTTCCAGGTCCGCTTCGACCCGGACCGTATCGCCCCCGAGGCACTGCTCGAAGCCGTGCGCGGCGCCGAGGCGCAGGCGGTGGCCGAGCGTTCGTTGCAGACGCGGATCATCGAGATCCCGGTGCTGTACAACGACCCCTGGACCCATGAAACCCTGATGCGCTTTC  >  NC_002947/3321171‑3321462

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