Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,538,457 C→T 9.6% L194L (CTG→CTA) ‡ 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,4570CT 9.6% 143.6 / 4.7 52L194L (CTG→CTA) ‡mrpABK(+)/H(+) antiporter subunit A/B
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base T (2/3);  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.48e-01

GGCCGACACCGGTGTGGGCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATCAG  >  NC_002947/2538313‑2538603
                                                                                                                                                |                                                                                                                                                  
ggCCGACACCGGTGTGGGCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAgg                                                                                                                                                 <  2:168720/148‑1 (MQ=255)
  ccGACACCGGTGTGGGCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAACCTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGAcc                                                                                                                                              <  7:177739/149‑1 (MQ=255)
                 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: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                                                                                      <  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)
                                                                                                                                                |                                                                                                                                                  
GGCCGACACCGGTGTGGGCGCTGCCATGGCGTGGGGCAGCCAGAACTGGAAGGGAAACTGTGCGCTCTTGGTCAGGGCGCCGATGAGTATCAGGGGCAGCAGCACCGGGTAAAGCGCATGCTGGCGGATGATATGCCCGGCAGCCAGGACCTTGTCCAGGTCATAGCTGCCGACCACATGGCCGAGCAGCAGGGCCCCTACCAACAGGCACAAGCCGCCTGCACCGGTTACCATCAGCGCCATATAGGCGCCGCGACGGGCGTCGGCGCGGTGGTGCCAGTAGCCGATCAG  >  NC_002947/2538313‑2538603

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