Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,614,535 G→T 8.5% E302D (GAG→GAT)  PP_4085 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,614,5350GT 8.5% 167.6 / 4.5 55E302D (GAG→GAT) PP_4085rhs family protein
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base T (2/3);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TCCGGGACATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCCGCG  >  NC_002947/4614390‑4614674
                                                                                                                                                 |                                                                                                                                           
tCCGGGACATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCAACGGGAGCAAGCTGAAGGAGGGt                                                                                                                                          >  4:305693/1‑149 (MQ=255)
      aCATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCggag                                                                                                                                    <  1:233907/149‑1 (MQ=255)
           cGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACAt                                                                                                                               >  8:202453/1‑149 (MQ=255)
                 ggCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACtt                                                                                                                         <  8:299856/149‑1 (MQ=255)
                 ggCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACtt                                                                                                                         <  2:114817/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                                      >  2:123563/1‑149 (MQ=255)
                                                                                                    aTTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATgg                                      >  5:150669/1‑149 (MQ=255)
                                                                                                    aTTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATgg                                      >  3:163091/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:17600/1‑106 (MQ=255)
                                                                                                                    gCTGTTGCCACCGGGAGCAAGCTGCAGG‑TGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACAt                                                                <  3:12535/106‑1 (MQ=255)
                                                                                                                    gCTGTTGCCACCGGGAGCAAGCTGCAGG‑TGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACAt                                                                <  3:17600/106‑1 (MQ=255)
                                                                                                                    gCTGTTGCCACCGGGAGCAAGCTGCAGG‑TGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACAt                                                                >  4:12535/1‑106 (MQ=255)
                                                                                                                     cTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTAt                     >  6:214919/1‑149 (MQ=255)
                                                                                                                        ttGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGATACGACAGCAGGTACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGaa                  <  3:51518/149‑1 (MQ=255)
                                                                                                                            cACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTg                                             >  5:318599/1‑118 (MQ=255)
                                                                                                                            cACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTg                                             <  6:318599/118‑1 (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)
                                                                                                                                                 |                                                                                                                                           
TCCGGGACATTCGCAACGGCAAGAGCAAGATTGCGATGTTTACCGGGATTATTGCGGGGATGCTTATCTCCCGAAGGTTTGGCCGCATCAAAGGGTGCACATTAGGCAACCCGGTTGCTGTTGCCACCGGGAGCAAGCTGCAGGAGGGTCCGGAGGACATCGACTTCAATCTCCCAGGCTTGCTGGCCATCCAATGGGCACGGCGCTACGACAGCAGGGACATGCGCAGCGACGGTCTACTGGGGATGGGTTGGAGTGTCCCCTATGAAGTGGAGCTGGCCCGCG  >  NC_002947/4614390‑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: ‑