Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,491,089 T→A 7.7% H101L (CAC→CTC)  clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,491,0890TA 7.7% 148.9 / 2.8 52H101L (CAC→CTC) clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base A (2/2);  total (25/27)
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

CCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAAC  >  NC_002947/3490946‑3491229
                                                                                                                                               |                                                                                                                                            
ccAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTgcgg                                                                                                                                          <  3:8715/148‑1 (MQ=255)
          tGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCg                                                                                                                               <  7:199274/149‑1 (MQ=255)
            cccaGGCGATCACCACCCGCCGCACGTCAGCCTCGTCAAGCCGAGCCTGCAGCAGGTGAGCGCCGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCCGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCggg                                                                                                                             <  4:201825/146‑1 (MQ=255)
              gAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGt                                                                                                                           >  6:300461/1‑149 (MQ=255)
                   gATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCggg                                                                                                                      >  7:42148/1‑149 (MQ=255)
                       accacGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTAt                                                                                                                   >  8:317859/1‑148 (MQ=255)
                        ccacGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCggg                                                                                                                             <  4:234848/137‑1 (MQ=255)
                        ccacGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCggg                                                                                                                             >  3:234848/1‑137 (MQ=255)
                        ccacGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTa                                                                                                                 >  7:47578/1‑149 (MQ=255)
                         cacGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTAc                                                                                                                <  7:98664/149‑1 (MQ=255)
                           cgcgACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCt                                                                                                              >  3:285470/1‑149 (MQ=255)
                              gACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTc                                                                                                           <  6:100542/149‑1 (MQ=255)
                               aCGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCt                                                                                                          <  5:300461/149‑1 (MQ=255)
                                 gCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCtt                                                                                                         >  7:304957/1‑148 (MQ=255)
                                   aGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTg                                                                                                        >  6:91492/1‑147 (MQ=255)
                                   aGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTg                                                                                                        >  1:212972/1‑147 (MQ=255)
                                   aGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGaa                                                                                                      <  7:87204/149‑1 (MQ=255)
                                    gTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAg                                                                                                     <  8:47578/149‑1 (MQ=255)
                                        aGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGc                                                                                                 <  4:285470/149‑1 (MQ=255)
                                              gTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAggg                                                                                           <  5:91492/149‑1 (MQ=255)
                                              gTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAgg                                                                                            <  7:317859/148‑1 (MQ=255)
                                               tCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGc                                                                                          >  3:132/1‑149 (MQ=255)
                                                 caGCCGAGCCTGCAGCAGCTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCg                                                                                         <  4:35861/147‑1 (MQ=255)
                                                 aaGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCg                                                                                         <  3:315685/148‑1 (MQ=255)
                                                   gCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGt                                                                                       >  6:103191/1‑148 (MQ=255)
                                                    cAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGa                                                                                     <  7:276232/149‑1 (MQ=255)
                                                       agCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGAc                                                                                    >  1:25770/1‑147 (MQ=255)
                                                       agCCTGCAGCAGGAGAGCGATGCGCACCTGGCCTAAGCCGTATTCGATCGAGGCCAGCAACCAGGCTTCCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTAATACCtat                                                                                                            >  5:284750/1‑121 (MQ=255)
                                                       agCCTGCAGCACGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTATTGTGCGGCGACTGCGCCGGGGTGCGGGTATTACGttt                                                                                                            <  6:284750/123‑1 (MQ=255)
                                                         ccTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGc                                                                                >  3:256280/1‑149 (MQ=255)
                                                                  ggTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTg                                                                        >  3:309287/1‑148 (MQ=255)
                                                                      aGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAgg                                                                    <  4:309287/148‑1 (MQ=255)
                                                                         gcTGCACACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGTCTTGCTCGATCATGCGGGTAATGAGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGt                                                                                       <  1:86492/126‑1 (MQ=255)
                                                                         gcTGCACACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGTCTTGCTCGATCATGCGGGTAATGAGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGt                                                                                       <  1:91076/126‑1 (MQ=255)
                                                                         gcTGCACACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGAGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGt                                                                                       >  2:91076/1‑126 (MQ=255)
                                                                         gcTGCACACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGAGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGt                                                                                       >  2:86492/1‑126 (MQ=255)
                                                                          cTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTc                                                               <  2:25770/149‑1 (MQ=255)
                                                                             cgcACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGc                                                            <  8:42148/149‑1 (MQ=255)
                                                                                  cTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGcgacga                                                       >  4:259407/1‑149 (MQ=255)
                                                                                            ccGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTg                                             <  6:320541/149‑1 (MQ=255)
                                                                                             cGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTgg                                            >  7:18198/1‑149 (MQ=255)
                                                                                                 ttCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGgcgc                                         <  4:132/148‑1 (MQ=255)
                                                                                                         aGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCTCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCa                                        >  6:309739/1‑141 (MQ=255)
                                                                                                         aGGCCAGCAACCAGGCTTGCTCGATCATGCGGGAAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCTCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCa                                        <  5:309739/141‑1 (MQ=255)
                                                                                                           gCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAg                                                                     >  6:369538/1‑110 (MQ=255)
                                                                                                           gCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAg                                                                     <  5:369538/110‑1 (MQ=255)
                                                                                                                   ccAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGc                      >  1:156780/1‑149 (MQ=255)
                                                                                                                    cAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGgcgc                                         <  3:3418/129‑1 (MQ=255)
                                                                                                                    cAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGAGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGgcgc                                         >  4:3418/1‑129 (MQ=255)
                                                                                                                              cGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGttgttg           >  7:217288/1‑149 (MQ=255)
                                                                                                                               gATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCgttgtt            <  8:18198/147‑1 (MQ=255)
                                                                                                                                       gggTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCAAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAAc  >  3:173484/1‑149 (MQ=255)
                                                                                                                                               tgCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTc                  >  2:34184/1‑125 (MQ=255)
                                                                                                                                               tgCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTc                  <  1:34184/125‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
CCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGTGAGCGCTGCGCACCTGGCCTACGCCGTATTCGATCGAGGCCAGCAACCAGGCTTGCTCGATCATGCGGGTAATGTGCGGCGACAGCGCCGGGGTGCGGGTATTACCTTTCTTGAAGGTGCCCAGGGCGGTGACCAGCTGCGCTTGCAGGCGTTCGGCGACGACCTGGTAGTGGCGCAGGATCGGCGCCAGGTCGCTGTCGTTGTTGTCGAGCAAC  >  NC_002947/3490946‑3491229

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