Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 587,076 C→G 7.9% P43R (CCG→CGG)  PP_0497 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947587,0760CG 7.9% 144.4 / 2.6 50P43R (CCG→CGG) PP_0497Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base G (2/2);  total (25/25)
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.08e-01

GCGTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTG  >  NC_002947/586929‑587221
                                                                                                                                                    |                                                                                                                                                 
gcgTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTcc                                                                                                                                                   >  1:285953/1‑148 (MQ=255)
     aCCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACGTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGc                                                                                                                                             >  1:72360/1‑149 (MQ=255)
       cTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCtt                                                                                                                                           <  2:189354/149‑1 (MQ=255)
        tGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCttt                                                                                                                                          >  5:222831/1‑149 (MQ=255)
         ggATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCttt                                                                                                                                          >  7:123285/1‑148 (MQ=255)
          gATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCtttgc                                                                                                                                        <  1:290530/149‑3 (MQ=255)
          gATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCttt                                                                                                                                          <  5:298150/147‑1 (MQ=255)
          gATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCa                                                                                                                                        <  4:160174/149‑1 (MQ=255)
            ttCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                      >  6:128051/1‑149 (MQ=255)
            ttCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAc                                                                                                                                      >  6:141004/1‑149 (MQ=255)
              ccGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAAcc                                                                                                                                     <  2:72360/148‑1 (MQ=255)
                  atatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCAGACGGCTTTCAACCTGAAg                                                                                                                                >  1:77589/1‑149 (MQ=255)
                   tatGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGc                                                                                                                               <  2:285953/149‑1 (MQ=255)
                           aTCACTTCGCTACAACCCTCAAACTTGCCTGCATGCACTACCGCTCTATTTCCGAAGTTTCCCGCCAGCTGTAGATCAACCGCGCGCAGTTCAACAAGTACCTGAGGGGTCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGc                                                                                                                       <  5:171629/149‑1 (MQ=255)
                                ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCtttgc                                                                                                                                        >  5:314032/1‑125 (MQ=255)
                                ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCGGACGGCtttgc                                                                                                                                        <  6:314032/127‑3 (MQ=255)
                                ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGa                                                                                                                     >  5:307433/1‑146 (MQ=255)
                                ttCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGa                                                                                                                     <  6:307433/146‑1 (MQ=255)
                                    cTACAAACCTCAAACTTGCCTGCAGTCACTNCCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACttcttc                                                                                                              >  1:27916/1‑149 (MQ=255)
                                    cTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACttcttc                                                                                                              >  1:32897/1‑149 (MQ=255)
                                    cTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACttcttc                                                                                                              <  3:47615/149‑1 (MQ=255)
                                       cAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCgg                                                                                                            <  7:64941/148‑1 (MQ=255)
                                         aaCCTGAAACTTGCCTTCAGTCAGTACCGCTCTATTTCGGTAGTTTGGCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTTCCTGTGCGGGGAAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGgtg                                                                                                          <  1:240275/149‑1 (MQ=255)
                                          aCCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTGCCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTc                                                                                                        >  1:297932/1‑149 (MQ=255)
                                             tCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                             >  6:222084/1‑141 (MQ=255)
                                             tCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCg                                                                                                             <  5:222084/141‑1 (MQ=255)
                                               aaaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAg                                                                                                     <  8:123285/147‑1 (MQ=255)
                                                aaCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGAt                                                                                                  <  5:128051/149‑1 (MQ=255)
                                                 aCTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCATCAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGAtt                                                                                                 <  4:106618/149‑1 (MQ=255)
                                                   ttGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTAc                                                                                               >  3:130214/1‑149 (MQ=255)
                                                   ttGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGAAGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTAc                                                                                               >  7:59154/1‑149 (MQ=255)
                                                                           tttCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTGCCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCa                                                                        <  2:297932/148‑1 (MQ=255)
                                                                               cGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCg                                                                   >  5:198119/1‑149 (MQ=255)
                                                                                           ccAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCg                                                       <  6:198119/149‑1 (MQ=255)
                                                                                                  tCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGcgcg                                                <  2:32897/149‑1 (MQ=255)
                                                                                                  tCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGcgcg                                                <  2:27916/149‑1 (MQ=255)
                                                                                                            gcgcgcAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGccc                                      >  3:227158/1‑149 (MQ=255)
                                                                                                              gcgcAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAAc                                                                                      >  4:33091/1‑99 (MQ=255)
                                                                                                              gcgcAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAAc                                                                                      <  3:33091/99‑1 (MQ=255)
                                                                                                               ctcAGTTCAACAAGTACCTGATCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAg                                   <  2:77589/147‑1 (MQ=255)
                                                                                                                    ttCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATAGCCCCAGGcc                                >  8:68906/1‑147 (MQ=255)
                                                                                                                    ccgatctaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGCAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                              <  5:223/142‑1 (MQ=255)
                                                                                                                     tCAACAAGTAGCTGATCTGGCAAAGCCGTCCGTCTGATTTCAACCTGAAGCGTATCGGCGAATTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCGGCCTAGAAGTTCGCCCGTCTGATCTGTGCGCGGGTATCGCCCCAGGCCCa                              <  8:59154/148‑1 (MQ=255)
                                                                                                                         caaGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAAc                         >  5:169291/1‑149 (MQ=255)
                                                                                                                            gTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCa                              >  6:223/1‑141 (MQ=255)
                                                                                                                               ccTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGt                    >  6:53944/1‑148 (MQ=255)
                                                                                                                               ccTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGt                    >  4:121435/1‑148 (MQ=255)
                                                                                                                                  gAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGAcc                >  8:1162/1‑149 (MQ=255)
                                                                                                                                       gttgaaaGCCGTCGGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGAt            >  2:184395/5‑149 (MQ=255)
                                                                                                                                               ccGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAgt    <  6:222831/149‑1 (MQ=255)
                                                                                                                                                 gTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAgttg  >  5:182596/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGTCACCTGGATTCCGCATATGAGCGATCACTTCGCTACAAACCTCAAACTTGCCTGCAGTCACTACCGCTCTATTTCCGAAGTTTGCCGCCAGCTGTCGATCAACCGCGCGCAGTTCAACAAGTACCTGAGCGGGCAAAGCCGTCCGACGGCTTTCAACCTGAAGCGTATCGGCGACTTCTTCGGTGTCGAGGATTACGAACTGAACCTGCCCGCCGAGCAGTTCGCCCGTCTGATCGGTGCGCGGGTATCGCCCCAGGCCCAGCAACCCAGTGACCCGATCAGCGAGTTG  >  NC_002947/586929‑587221

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: