Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,538,465 T→G 9.4% K192Q (AAG→CAG)  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,4650TG 9.4% 160.1 / 6.2 53K192Q (AAG→CAG) mrpABK(+)/H(+) antiporter subunit A/B
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base G (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 = 1.00e+00

GCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATCAGCAG  >  NC_002947/2538330‑2538606
                                                                                                                                       |                                                                                                                                             
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                                                                                                                          <  2:244067/133‑3 (MQ=38)
                        aGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCTATGACCTGGACAAGGTCCTGGCTGCCggg                                                                                                                          >  1:244067/1‑131 (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                                                                                                              <  6:87884/126‑1 (MQ=255)
                                           tgCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGag                                                                                                              >  5:87884/1‑126 (MQ=255)
                                             cgcTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGca                                                                                      <  6:336841/148‑1 (MQ=255)
                                             cgcTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGca                                                                                      >  6:312510/1‑148 (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                                                                                         <  1:176409/103‑1 (MQ=255)
                                                                                       ccGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAg                                                                                         >  2:176409/1‑103 (MQ=255)
                                                                                          ggTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGagc                                                                                                             >  1:204715/1‑80 (MQ=255)
                                                                                          ggTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGagc                                                                                                             <  2:204715/80‑1 (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)
                                                                                                                                       |                                                                                                                                             
GCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATCAGCAG  >  NC_002947/2538330‑2538606

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