Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,490,195 G→A 7.7% A399V (GCG→GTG)  clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,490,1950GA 7.7% 151.2 / 2.2 52A399V (GCG→GTG) clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base G (19/29);  new base A (2/2);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

TTCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTACTTCTTGTATTCGGCCC  >  NC_002947/3490047‑3490339
                                                                                                                                                    |                                                                                                                                                  
ttCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCg                                                                                                                                                    >  3:108703/1‑149 (MQ=255)
ttCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCg                                                                                                                                                    >  7:179463/1‑149 (MQ=255)
ttCCAGTGCGCCGGCCAGCGAAGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCg                                                                                                                                                    >  1:238969/1‑149 (MQ=255)
           cGGGCAGCGAAGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTgtgg                                                                                                                                         <  2:57591/149‑1 (MQ=255)
              gCAGCGAAGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGc                                                                                                                                      <  4:108703/149‑1 (MQ=255)
                    aaGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCCCCGCCACCTTGTGGTGCTcgcg                                                                                                                                 >  2:65030/1‑148 (MQ=255)
                        ttGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCt                                                                                                                            >  5:211208/1‑149 (MQ=255)
                           cgcCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGc                                                                                                                         <  5:229406/149‑1 (MQ=255)
                              cAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGa                                                                                                                      <  2:193675/149‑1 (MQ=255)
                               aGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAg                                                                                                                     >  2:43148/1‑149 (MQ=255)
                                    aaTGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGc                                                                                                                >  1:45746/1‑149 (MQ=255)
                                        cgcgcgcAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCcgcg                                                                                                            <  3:297681/149‑1 (MQ=255)
                                         gcgcgcAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCgcgc                                                                                                           <  2:45746/149‑1 (MQ=255)
                                          cgcgcAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCgca                                                                                                          <  1:284821/149‑1 (MQ=255)
                                             gcAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCgcagca                                                                                                      >  6:347323/1‑149 (MQ=255)
                                              cAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCAt                                                                                                     <  3:169548/149‑1 (MQ=255)
                                                ggCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCAtgt                                                                                                   <  4:41712/149‑1 (MQ=255)
                                                  ccGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCAtgtg                                                                                                  >  1:301213/1‑148 (MQ=255)
                                                        cGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGcc                                                                                           <  7:109834/149‑1 (MQ=255)
                                                         gAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGcc                                                                                           <  2:252633/148‑1 (MQ=255)
                                                             aCGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTc                                                                                      <  8:40821/149‑1 (MQ=255)
                                                             aCGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCACCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTc                                                                                      >  7:127060/1‑149 (MQ=255)
                                                             aCGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCACCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTc                                                                                      >  7:131495/1‑149 (MQ=255)
                                                             aCGCTGACCGCCGTGTCGGGCAGCTGGCGGCCGGTTATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTc                                                                                      <  7:139643/149‑1 (MQ=255)
                                                                 tGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCAGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTtcg                                                                                  >  2:333511/1‑149 (MQ=255)
                                                                    ccGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTtcgtcg                                                                               <  5:212087/149‑1 (MQ=255)
                                                                    ccGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTtcgtc                                                                                <  5:140428/148‑1 (MQ=255)
                                                                     cGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCgg                                                                              >  5:219565/1‑149 (MQ=255)
                                                                      gCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCggg                                                                             <  2:301213/149‑1 (MQ=255)
                                                                                 cAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCACCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACtt                                                                  <  8:131495/149‑1 (MQ=255)
                                                                                 cAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCACCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACtt                                                                  <  8:127060/149‑1 (MQ=255)
                                                                                  aGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACttt                                                                 <  1:43148/149‑1 (MQ=255)
                                                                                    cTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGa                                                               <  6:169263/149‑1 (MQ=255)
                                                                                       gcggcCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGac                                                              <  6:211208/147‑1 (MQ=255)
                                                                                         ggcCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCAACTTTGaccacc                                                          >  8:239676/1‑149 (MQ=255)
                                                                                             ggTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTgg                                                      >  8:218462/1‑148 (MQ=255)
                                                                                               tGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAg                                                   <  7:218462/149‑1 (MQ=255)
                                                                                                 aTGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAgcg                                                 >  6:264026/1‑149 (MQ=255)
                                                                                                        ggTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGctc                                                                          >  5:343682/1‑118 (MQ=255)
                                                                                                        ggTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGctc                                                                          <  6:343682/118‑1 (MQ=255)
                                                                                                        ggTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCggcgggc                                           >  3:141065/1‑148 (MQ=255)
                                                                                                          ttGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAg                                        <  5:347323/149‑1 (MQ=255)
                                                                                                                aGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGa                                                               >  6:111124/1‑121 (MQ=255)
                                                                                                                aGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGa                                                               <  5:111124/121‑1 (MQ=255)
                                                                                                                 ggCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGgcggcg                                 <  5:264026/149‑1 (MQ=255)
                                                                                                                   cGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTc                               >  1:40439/1‑149 (MQ=255)
                                                                                                                         gCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGa                                                               <  1:166145/112‑1 (MQ=255)
                                                                                                                         gCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGa                                                               >  2:166145/1‑112 (MQ=255)
                                                                                                                                 cAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTAct                 <  2:238969/149‑1 (MQ=255)
                                                                                                                                     gcTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTActtctt             <  8:179463/149‑1 (MQ=255)
                                                                                                                                      cTTCGTCCATCCCCGCCACCTTGTGGTGCTCGCGCATCTTGCCGACGATGCCGCGGAGCATGTGAATGGGCTTGGCTTCCTCCGGCTCTTCCACTTTGCCCACCACGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTActtctt             <  1:65030/149‑1 (MQ=255)
                                                                                                                                          gTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCTGTGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTACTTCTTGTAt         >  1:93436/1‑149 (MQ=255)
                                                                                                                                             cATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGTGCGAGGGCGGCGTCCTTCTCGAAGTACTTCTTGTATTCgg     >  3:201550/1‑149 (MQ=255)
                                                                                                                                               tCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTACTTCTTGTATTCGGcc   >  3:51766/1‑149 (MQ=255)
                                                                                                                                                 aCCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTACTTCTTGTATTCGGccc  <  7:239676/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAATGCGCGCGCAGGCCGTGTCGAGCACGCTGACCGCCTTGTCGGGCAGCTGGCGGCCGGTGATGTAGCGGTTGGACAGGCGTACCGCCTGCACCAGCGCTTCGTCCATCACCGCCACCTTGTGGTGCTCGCGCATCTTGCCGAGCAGGCCGCGCAGCATGTGAATGGCCTTGTCTTCGTCGGGCTCTTCCACTTTGACCACCTGGAAGCGGCGGGCGAGGGCGGCGTCCTTCTCGAAGTACTTCTTGTATTCGGCCC  >  NC_002947/3490047‑3490339

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: