Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,614,534 Δ1 bp 7.3% coding (905/4593 nt) PP_4085 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,614,5340A. 7.3% 173.6 / 7.4 55coding (905/4593 nt)PP_4085rhs family protein
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base . (2/2);  total (29/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.94e-01

ACGGTCCGGGACATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCCGCG  >  NC_002947/4614386‑4614674
                                                                                                                                                    |                                                                                                                                            
aCGGTCCGGGACATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCagga                                                                                                                                              <  7:216163/149‑1 (MQ=255)
    tCCGGGACATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCAACGGGAGCAAGCTGAAGGAGGGt                                                                                                                                          >  4:305693/1‑149 (MQ=255)
          aCATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCggag                                                                                                                                    <  1:233907/149‑1 (MQ=255)
               cGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACAt                                                                                                                               >  8:202453/1‑149 (MQ=255)
                     ggCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACtt                                                                                                                         <  2:114817/149‑1 (MQ=255)
                     ggCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACtt                                                                                                                         <  8:299856/149‑1 (MQ=255)
                        aagagcaagaTTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCaa                                                                                                                      <  6:194150/149‑1 (MQ=255)
                                gaTTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAg                                                                                                              <  6:190764/149‑1 (MQ=255)
                                   tGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCt                                                                                                           >  8:211407/1‑149 (MQ=255)
                                            tACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCt                                                                                                       >  4:352613/1‑144 (MQ=255)
                                            tACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCt                                                                                                       <  3:352613/144‑1 (MQ=255)
                                             aCCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCAt                                                                                                 >  1:35928/1‑149 (MQ=255)
                                              ccGTGATTATTGCATGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGATCAAGCTGCATGATGGTCGGGATGACATCGACTTCAATCTCCCAGGCTTGCTTGCCATc                                                                                                <  4:316833/149‑1 (MQ=255)
                                                   attattGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTTCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATg                                                                                           <  1:275572/149‑1 (MQ=255)
                                                    ttattGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATgg                                                                                          >  6:248769/1‑149 (MQ=255)
                                                        tGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCAc                                                                                      >  6:306259/1‑149 (MQ=255)
                                                              gATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCt                                                                                <  2:35928/149‑1 (MQ=255)
                                                                   ttATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGa                                                                            >  2:278877/1‑148 (MQ=255)
                                                                   ttATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGa                                                                            >  2:331844/1‑148 (MQ=255)
                                                                          ccGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAgg                                                                     <  7:202453/148‑1 (MQ=255)
                                                                                gTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATgc                                                              >  2:143950/1‑149 (MQ=255)
                                                                                   tGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATgcg                                                             >  8:98801/1‑147 (MQ=255)
                                                                                      ccGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGc                                                         >  1:255278/1‑148 (MQ=255)
                                                                                          aTCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACgg                                                    <  5:248769/149‑1 (MQ=255)
                                                                                           tCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGt                                                   >  1:335853/1‑149 (MQ=255)
                                                                                            cAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGCACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAAAGGGCACGGCGATACGACAGAAGGGACATTCGCTGCGACGGTc                                                  >  8:162260/1‑149 (MQ=255)
                                                                                               aGGGTGCACATTAGGCAACCCGGTTTCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTAc                                               <  3:305693/149‑1 (MQ=255)
                                                                                                      acaTTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGCCGGTCTACGGGGGag                                        >  4:225466/1‑148 (MQ=255)
                                                                                                        aTTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATgg                                      >  3:163091/1‑149 (MQ=255)
                                                                                                        aTTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATgg                                      >  5:150669/1‑149 (MQ=255)
                                                                                                        aTTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATgg                                      >  2:123563/1‑149 (MQ=255)
                                                                                                          tAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGAAATGCGCAGCGAAGGTCTACTGGGGATGGGt                                    >  4:51518/1‑149 (MQ=255)
                                                                                                            ggCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTg                                  <  1:278877/149‑1 (MQ=255)
                                                                                                              cAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGa                                <  4:163091/149‑1 (MQ=255)
                                                                                                              cAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGa                                <  5:306259/149‑1 (MQ=255)
                                                                                                               aaCCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAAT‑GGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAgt                              >  3:316507/1‑149 (MQ=255)
                                                                                                                 cccGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGt                                                   <  2:146287/127‑1 (MQ=255)
                                                                                                                 cccGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGt                                                   >  1:146287/1‑127 (MQ=255)
                                                                                                                  ccGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAgtgt                            >  4:254979/1‑149 (MQ=255)
                                                                                                                   cGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTc                           <  2:255278/149‑1 (MQ=255)
                                                                                                                        gCTGTTGCCACCGGGAGCAAGCTGCAGG‑TGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACAt                                                                >  4:12535/1‑106 (MQ=255)
                                                                                                                        gCTGTTGCCACCGGGAGCAAGCTGCAGG‑TGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACAt                                                                >  4:17600/1‑106 (MQ=255)
                                                                                                                        gCTGTTGCCACCGGGAGCAAGCTGCAGG‑TGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACAt                                                                <  3:17600/106‑1 (MQ=255)
                                                                                                                        gCTGTTGCCACCGGGAGCAAGCTGCAGG‑TGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACAt                                                                <  3:12535/106‑1 (MQ=255)
                                                                                                                         cTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTAt                     >  6:214919/1‑149 (MQ=255)
                                                                                                                            ttGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGATACGACAGCAGGTACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGaa                  <  3:51518/149‑1 (MQ=255)
                                                                                                                                cACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTg                                             <  6:318599/118‑1 (MQ=255)
                                                                                                                                cACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTg                                             >  5:318599/1‑118 (MQ=255)
                                                                                                                                 aCCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGa             <  7:211407/149‑1 (MQ=255)
                                                                                                                                   cGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGc           <  1:331844/149‑1 (MQ=255)
                                                                                                                                       aGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGc       <  1:123563/149‑1 (MQ=255)
                                                                                                                                         cAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGccc     >  1:146005/1‑149 (MQ=255)
                                                                                                                                          aaGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCcg    >  7:165590/1‑149 (MQ=255)
                                                                                                                                           aGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCcgc   <  6:150669/149‑1 (MQ=255)
                                                                                                                                           aGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCcgc   <  4:99854/149‑1 (MQ=255)
                                                                                                                                            gCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCcgcg  >  8:133843/1‑149 (MQ=255)
                                                                                                                                                 aggaggGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGa                                >  8:227228/1‑114 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
ACGGTCCGGGACATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCCGCG  >  NC_002947/4614386‑4614674

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