Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,336,884 C→T 8.3% H30Y (CAC→TAC)  PP_2934 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,336,8840CT 8.3% 133.2 / 2.6 48H30Y (CAC→TAC) PP_2934hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/24);  new base T (2/2);  total (22/26)
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.36e-01

CACAACGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGCCGAAGAGGCGCCCGTCGCCATGGGGCGATTCCC  >  NC_002947/3336737‑3337032
                                                                                                                                                   |                                                                                                                                                    
cacaACGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCa                                                                                                                                                     <  1:324984/149‑1 (MQ=255)
  caACGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCAcc                                                                                                                                                   <  1:69364/149‑1 (MQ=255)
  caACGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCCTGGGTTGCCGGCCTGcccc                                                                                                                                                   >  2:332894/1‑149 (MQ=255)
   aaCGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGTACCt                                                                                                                                                  >  2:95682/1‑149 (MQ=255)
   aaCGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGTACCt                                                                                                                                                  >  6:83524/1‑149 (MQ=255)
     cGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTa                                                                                                                                                >  4:286026/1‑149 (MQ=255)
      gccggccgGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTAt                                                                                                                                               >  6:43504/1‑149 (MQ=255)
         ggcgggTTTTATGCTGGGTTATCCCTTCGGCCTGCACGGAACCGCTGGCCTAGTACCCTTGCACGCTGCGGGTTGGCACTGGTCCTGCGCGGCGACCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGAc                                                                                                                                            <  1:90395/149‑1 (MQ=255)
         ggccggTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCtt                                                                                                                                                 <  6:314895/144‑1 (MQ=255)
         ggccggTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCtt                                                                                                                                                 >  5:314895/1‑144 (MQ=255)
              gTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGTACCTTATGACCTTga                                                                                                                                       <  1:95682/149‑1 (MQ=255)
              gTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGTACCTTATGACCTTga                                                                                                                                       <  5:83524/149‑1 (MQ=255)
               ttttATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTgac                                                                                                                                      >  8:320324/1‑149 (MQ=255)
                      cTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTTACCGAccc                                                                                                                                <  4:328234/148‑1 (MQ=255)
                                ccTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCTGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAgcc                                                                                                                     <  1:332894/149‑1 (MQ=255)
                                 cTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTgaccgac                                                                                                                                  <  3:197997/135‑1 (MQ=255)
                                 cTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTgaccgac                                                                                                                                  >  4:197997/1‑135 (MQ=255)
                                      tCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGa                                                                                                               >  5:226687/1‑149 (MQ=255)
                                        cTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATg                                                                                                             >  3:272847/1‑149 (MQ=255)
                                          gCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGc                                                                                                            >  4:303316/1‑148 (MQ=255)
                                            aaGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAg                                                                                                         >  5:106941/1‑149 (MQ=255)
                                               gAACCGCTGCCCGATGACCCTTGCACGCTGCTTGTTGTCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGccc                                                                                                      <  2:197475/149‑1 (MQ=255)
                                                aaCCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCt                                                                                                     <  6:226687/149‑1 (MQ=255)
                                                          cGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTAc                                                                                           >  7:117201/1‑149 (MQ=255)
                                                           gATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTAcc                                                                                          <  7:320324/149‑1 (MQ=255)
                                                               aCCCTTGCAGGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCgg                                                                                       <  6:92461/148‑1 (MQ=255)
                                                                   ttGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGtt                                                                                  >  5:231380/1‑149 (MQ=255)
                                                                      cACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTgg                                                                                >  5:197358/1‑148 (MQ=255)
                                                                        cGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGCTGCCAGCCCTGGCGTTCTACCCGGCCGTTGGCa                                                                              >  7:244212/1‑148 (MQ=255)
                                                                             cGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAAccc                                                                        <  7:116920/149‑1 (MQ=255)
                                                                              gtgtTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCg                                                                       <  3:303316/149‑1 (MQ=255)
                                                                                gtTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTAc                                                                                           <  4:5899/127‑1 (MQ=255)
                                                                                gtTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTAc                                                                                           >  3:5899/1‑127 (MQ=255)
                                                                                      gcTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGcc                                                               >  8:178913/1‑149 (MQ=255)
                                                                                          gTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGta                                                            <  2:9144/148‑1 (MQ=255)
                                                                                                cTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACgg                                                     <  3:286026/149‑1 (MQ=255)
                                                                                                atCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACgg                                                     <  4:100111/148‑1 (MQ=255)
                                                                                                 tCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACggg                                                    >  3:40780/1‑149 (MQ=255)
                                                                                                   gggggcCTGTTGGCCCGCGCGGAGGCCGGGCGGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTGGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAG‑CGTATAGACGGGt                                                   <  4:272847/144‑1 (MQ=255)
                                                                                                   ggcggcCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGt                                                   <  6:201686/148‑1 (MQ=255)
                                                                                                     cggcCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATc                                                <  8:117201/149‑1 (MQ=255)
                                                                                                       gcCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATccc                                              <  7:178913/149‑1 (MQ=255)
                                                                                                            ttGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGa                                         <  4:40780/149‑1 (MQ=255)
                                                                                                            ttGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGa                                         >  8:93510/1‑149 (MQ=255)
                                                                                                                 ccGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTgg                                    <  6:106941/149‑1 (MQ=255)
                                                                                                                  cgcgCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGc                                   >  6:159463/1‑149 (MQ=255)
                                                                                                                           gCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGCCGAAGAGgc                          >  1:55028/1‑149 (MQ=255)
                                                                                                                             cgcgCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGCCGAAGAGgcgc                        <  1:162534/149‑1 (MQ=255)
                                                                                                                              gcgcCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGCCGAAGAGGCGcc                       >  5:300735/1‑149 (MQ=255)
                                                                                                                                                   cACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGCCGAAGAGGCGCCCGTCGCCATGGGGCGATTccc  >  6:141125/1‑149 (MQ=255)
                                                                                                                                                   cACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGCCGAAGAGGCGCCCGTCGCCATGGGGCGATTccc  >  5:60905/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACAACGCCGGCCGGTTTTATGCTGGGTTATCCCTTCGTCCTGCAAGGAACCGCTGGCCGATGACCCTTGCACGCTGCGTGTTGGCGCTGGTCCTGCTCGGCGGCCTGTTGGCCCGCGCCGAGGCCGCGCCGTGGGTTGCCGGCCTGCACCTTATGACCTTGACCGACCCGGTCGATAGCCGCCCGATGCAAGCCCTGGCGTTCTACCCGGCCGTTGGCAAAACCCGTATCAGCCGTATAGACGGGTATCCCGTCGACGTGGCCGAAGAGGCGCCCGTCGCCATGGGGCGATTCCC  >  NC_002947/3336737‑3337032

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