Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,489,714 G→C 7.0% R574G (CGC→GGC)  clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,489,7140GC 7.0% 181.1 / 5.3 57R574G (CGC→GGC) clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base G (30/23);  new base C (2/2);  total (32/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 = 8.78e-01

GCCGACGCCGCTGGGGCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGT  >  NC_002947/3489565‑3489861
                                                                                                                                                     |                                                                                                                                                   
gCCGACGCCGCTGGGGCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGgcg                                                                                                                                                     <  1:625803/150‑1 (MQ=255)
       ccgcTGGGGCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACg                                                                                                                                              <  2:248952/150‑1 (MQ=255)
         gcTGGGGCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGcttggt                                                                                                                                            <  1:238699/150‑6 (MQ=255)
         gcTGGGGCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGt                                                                                                                                            >  2:969478/1‑150 (MQ=255)
              ggCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGa                                                                                                                                       <  2:223381/150‑1 (MQ=255)
               gCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATTCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGAt                                                                                                                                      >  2:658861/1‑150 (MQ=255)
                      agcagGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCa                                                                                                                               <  2:492045/150‑1 (MQ=255)
                       gcagGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGACCGCCGACGCTGCACGGTATCGATTTCGTCAc                                                                                                                              >  1:413836/1‑150 (MQ=255)
                          ggAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACgca                                                                                                                           >  1:423639/1‑150 (MQ=255)
                                  ccAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTg                                                                                                                   <  1:1298892/150‑1 (MQ=255)
                                             gTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGAt                                                                                                        >  1:845503/1‑150 (MQ=255)
                                              ttCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCg                                                                                                                                  <  1:46461/123‑1 (MQ=255)
                                              ttCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCg                                                                                                                                  >  2:46461/1‑123 (MQ=255)
                                                 gggTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCg                                                                                                    >  1:1578793/1‑150 (MQ=255)
                                                   gTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGt                                                                                                  >  1:1469902/1‑150 (MQ=255)
                                                      ttCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCa                                                                                               >  1:1314821/1‑150 (MQ=255)
                                                        ccATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGc                                                                                             <  2:1370331/150‑1 (MQ=255)
                                                                gcTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTc                                                                                                                                >  1:1184825/1‑107 (MQ=255)
                                                                gcTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTc                                                                                                                                <  2:1184825/107‑1 (MQ=255)
                                                                  tCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACg                                                                                                                                              >  1:1228058/1‑91 (MQ=255)
                                                                  tCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACg                                                                                                                                              <  2:1228058/91‑1 (MQ=255)
                                                                  tCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGCGCCTGCCGGCCCTGcttg                                                                                                                                              <  1:396983/91‑4 (MQ=255)
                                                                  tCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGCGCCTGCCGGCCCTGcttg                                                                                                                                              >  2:396983/1‑88 (MQ=255)
                                                                  tCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGCGCCTGCCGGCCCTGcttg                                                                                                                                              <  2:1238767/91‑4 (MQ=255)
                                                                  tCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGCGCCTGCCGGCCCTGcttg                                                                                                                                              >  1:1238767/1‑88 (MQ=255)
                                                                    gggAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCt                                                                                 >  2:839392/1‑150 (MQ=255)
                                                                     ggAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTg                                                                                >  1:571600/1‑150 (MQ=255)
                                                                          tCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGa                                                                           <  1:263788/150‑1 (MQ=255)
                                                                            ttGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATg                                                                         >  1:814300/1‑150 (MQ=255)
                                                                                   gcTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGcgccgc                                                                  >  1:354457/1‑150 (MQ=255)
                                                                                     ttGCCGATCTCATGCAGGGCGTGGTCCTGGCGCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATgg                                                                        >  2:844424/1‑142 (MQ=255)
                                                                                     ttGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATgg                                                                        <  1:844424/142‑1 (MQ=255)
                                                                                       gCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTc                                                              >  1:163709/1‑150 (MQ=255)
                                                                                       gCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTc                                                              >  2:1071585/1‑150 (MQ=255)
                                                                                         cGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCGGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCggg                                                                                                          >  2:1061775/1‑104 (MQ=255)
                                                                                         cGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCggg                                                                                                          <  1:1061775/104‑1 (MQ=255)
                                                                                         cGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGa                                                            >  1:491709/1‑150 (MQ=255)
                                                                                           aTCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGAcc                                                          <  1:1071585/150‑1 (MQ=255)
                                                                                                aTGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGAt                                                                                                                                      >  2:38198/1‑69 (MQ=255)
                                                                                                aTGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGAt                                                                                                                                      <  1:38198/69‑1 (MQ=255)
                                                                                                   cAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTg                                                  <  1:338615/150‑1 (MQ=255)
                                                                                                         gTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGt                                                               >  2:1415072/1‑131 (MQ=255)
                                                                                                         gTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGt                                                               <  1:1415072/131‑1 (MQ=255)
                                                                                                            gTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGc                                                                                             <  1:1497012/98‑1 (MQ=255)
                                                                                                            gTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGc                                                                                             >  2:1497012/1‑98 (MQ=255)
                                                                                                            gTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTgg                                         >  2:570579/1‑150 (MQ=255)
                                                                                                             tCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCCGCGACCAGTGGc                                        >  2:1274765/1‑150 (MQ=255)
                                                                                                                          cgcgCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCa                           <  1:839392/150‑1 (MQ=255)
                                                                                                                                 aCCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGAt                                                                          <  2:995685/96‑1 (MQ=255)
                                                                                                                                 aCCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGAt                                                                          >  1:995685/1‑96 (MQ=255)
                                                                                                                                 aCCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACgcgcgc                    >  1:684776/1‑150 (MQ=255)
                                                                                                                                   cAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCg                                                                                                    <  1:1451145/68‑1 (MQ=255)
                                                                                                                                   cAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCg                                                                                                    >  2:1451145/1‑68 (MQ=255)
                                                                                                                                   cAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCa                  >  1:345693/1‑150 (MQ=255)
                                                                                                                                   cAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCa                  <  2:568491/150‑1 (MQ=255)
                                                                                                                                         gggCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCgg            >  2:1274811/1‑150 (MQ=255)
                                                                                                                                               gCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCa      <  2:1578793/150‑1 (MQ=255)
                                                                                                                                                cAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAg     <  2:1469902/150‑1 (MQ=255)
                                                                                                                                                   gcgcTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGgc                                                                       <  1:862769/81‑1 (MQ=255)
                                                                                                                                                   gcgcTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGgc                                                                       >  2:862769/1‑81 (MQ=255)
                                                                                                                                                   gcgcTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGt  <  2:814300/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCGACGCCGCTGGGGCCGAGCAGCAGGAACACGCCAATCGGTTTGTTCGGGTCTTCCATGCGCGCTCGGGAAATCTTGATGCGCTTGCCGATCTCATGCAGGGCGTGGTCCTGGCCCAGCACGCGCTCACCAAGCAGGGCCGGCAGGCGCTGCACGGTATCGATTTCGTCACGCAGCATCTTGCCCAGCGGGATGCCGGTCCAGCCACTGATGACCTGGGCGATGGCGCCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGT  >  NC_002947/3489565‑3489861

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: