Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,538,467 T→A 9.4% D191V (GAC→GTC)  mrpAB ← K(+)/H(+) antiporter subunit A/B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,538,4670TA 9.4% 150.1 / 5.7 53D191V (GAC→GTC) mrpABK(+)/H(+) antiporter subunit A/B
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base A (2/3);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATCAGCAGGAAGGAGAA  >  NC_002947/2538330‑2538615
                                                                                                                                         |                                                                                                                                                    
gcgcTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAg                                                                                                                                           <  3:6848/149‑1 (MQ=255)
 cgcTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGc                                                                                                                                          >  1:70015/1‑149 (MQ=255)
  gcTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCt                                                                                                                                         >  2:325803/1‑149 (MQ=255)
  gcTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCt                                                                                                                                         >  2:218865/1‑149 (MQ=255)
      ccATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCTATGACCTGGACAAGGTCCTGGCTGcc                                                                                                                                      <  2:75202/148‑1 (MQ=255)
        aTGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGAc                                                                                                                                   >  5:87403/1‑149 (MQ=255)
         tGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCTATGACCTGGACAAGGTCCTGGCTGCCgggc                                                                                                                                  <  7:7890/149‑4 (MQ=38)
           gCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACca                                                                                                                                 <  1:218865/148‑1 (MQ=255)
                gggCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGc                                                                                                                           <  8:225526/149‑1 (MQ=255)
                    aGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGag                                                                                                                       <  2:70015/149‑1 (MQ=255)
                      ccAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGag                                                                                                                       >  4:5732/1‑147 (MQ=255)
                        aGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCTATGACCTGGACAAGGTCCTGGCTGCCggg                                                                                                                                   >  1:244067/1‑131 (MQ=38)
                        aGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCTATGACCTGGACAAGGTCCTGGCTGCCggg                                                                                                                                   <  2:244067/133‑3 (MQ=38)
                        aGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGagcagc                                                                                                                   <  2:376681/149‑1 (MQ=255)
                              ggAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGcc                                                                                                             >  7:245177/1‑149 (MQ=255)
                                       aCTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCaa                                                                                                     >  3:191091/1‑148 (MQ=255)
                                        cTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCATCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACa                                                                                                   <  1:325803/149‑1 (MQ=255)
                                           tgCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGag                                                                                                                       >  5:87884/1‑126 (MQ=255)
                                           tgCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGag                                                                                                                       <  6:87884/126‑1 (MQ=255)
                                             cgcTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGca                                                                                               >  6:312510/1‑148 (MQ=255)
                                             cgcTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGca                                                                                               <  6:336841/148‑1 (MQ=255)
                                                     gTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAgccgcc                                                                                      >  6:87245/1‑149 (MQ=255)
                                                         gggCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCAACCAACAGGCACAAGCCGCCTGCa                                                                                  >  3:63795/1‑149 (MQ=255)
                                                              ccGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGt                                                                             <  6:198561/149‑1 (MQ=255)
                                                              cagATGAGTATCAGGGGCAGCATCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCTGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGt                                                                             <  4:50050/147‑1 (MQ=255)
                                                               cGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGt                                                                             >  6:235718/1‑148 (MQ=255)
                                                                 aTGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTAc                                                                          >  4:221025/1‑149 (MQ=255)
                                                                   gAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCa                                                                        >  4:345296/1‑149 (MQ=255)
                                                                          aGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAgcgc                                                                 >  4:28060/1‑149 (MQ=255)
                                                                            gggCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCa                                                               <  2:116257/149‑1 (MQ=255)
                                                                              gcagcagcaCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCata                                                             <  5:235718/149‑1 (MQ=255)
                                                                                agcagcaCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCAtata                                                           <  4:191091/149‑1 (MQ=255)
                                                                                 gcagcaCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAg                                                          <  3:28060/149‑1 (MQ=255)
                                                                                     caCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTAcc                                                                         <  1:13687/130‑1 (MQ=255)
                                                                                     caCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTAcc                                                                         >  2:13687/1‑130 (MQ=255)
                                                                                      aCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCGAGGTCATa                                                                                                                                            >  3:399593/1‑62 (MQ=255)
                                                                                      aCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATa                                                                                                                                            <  4:399593/62‑1 (MQ=255)
                                                                                      aCCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGcgcc                                                     <  3:5732/149‑1 (MQ=255)
                                                                                       ccGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAg                                                                                                  >  2:176409/1‑103 (MQ=255)
                                                                                       ccGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAg                                                                                                  <  1:176409/103‑1 (MQ=255)
                                                                                          ggTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGagc                                                                                                                      <  2:204715/80‑1 (MQ=255)
                                                                                          ggTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGagc                                                                                                                      >  1:204715/1‑80 (MQ=255)
                                                                                                gcgcATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCg                                           >  6:249522/1‑149 (MQ=255)
                                                                                                  gcATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGATTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTc                                         <  8:245177/149‑1 (MQ=255)
                                                                                                    aTGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCgg                                       <  3:345296/149‑1 (MQ=255)
                                                                                                        tGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGcgcg                                   >  4:390990/1‑149 (MQ=255)
                                                                                                          gCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCggt                                 <  5:312510/149‑1 (MQ=255)
                                                                                                                    tGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAg                         >  1:117469/1‑147 (MQ=255)
                                                                                                                     gCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAg                      >  4:217581/1‑149 (MQ=255)
                                                                                                                          gCAGCTATGACCTGGACAAGGTCCTGGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCTCCGCGACTGGCGTCGGCGCGGTGGTGCCAGTAGCCGa                  >  5:326742/1‑148 (MQ=255)
                                                                                                                               cAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATcag              >  2:23427/1‑147 (MQ=255)
                                                                                                                               cAGGACCTTGTCCAGGTCATAGCTGCCGACAACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATcag              >  8:183183/1‑147 (MQ=255)
                                                                                                                                aGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATcagcag           <  7:183183/149‑1 (MQ=255)
                                                                                                                                         tCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGACTGCACCGGTTACCATCAGAGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGAGCAGCAGGAAGGAGaa  >  7:264781/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
GCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATCAGCAGGAAGGAGAA  >  NC_002947/2538330‑2538615

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