Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,535,042 C→T 7.8% A262T (GCC→ACC)  mrpD ← K(+)/H(+) antiporter subunit D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,535,0420CT 7.8% 143.7 / 3.2 51A262T (GCC→ACC) mrpDK(+)/H(+) antiporter subunit D
Reads supporting (aligned to +/- strand):  ref base C (20/27);  new base T (2/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAGGCCATGGCCAGGATTGCCGCGCCGGCATG  >  NC_002947/2534901‑2535186
                                                                                                                                             |                                                                                                                                                
gCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGc                                                                                                                                          <  2:163860/150‑1 (MQ=255)
  tCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCg                                                                                                                                        >  1:1094045/1‑150 (MQ=255)
   ccAGTCGCTGAGCTGCCAGGACCGACACCGCTGCCCCACCCAGGGCGACCAGGCCGCCGTCAACACGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGa                                                                                                                                       <  2:1082255/150‑1 (MQ=255)
   ccAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGa                                                                                                                                       >  1:719795/1‑150 (MQ=255)
        ggcTGAGCTGCCAGGAGCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTg                                                                                                                                  <  2:245457/149‑1 (MQ=255)
        cgcTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGGCCGGAGAACAGCAGGGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTg                                                                                                                                  <  2:194946/150‑1 (MQ=255)
        cgcTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTg                                                                                                                                  <  2:434529/150‑1 (MQ=255)
                        tCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAACCCAGCCATTCGCCGCCTAAATGCGCCGAGGCCAAGGAATCCCACGAGTACAGCAAATTGCAGAGTGGCAGCTTGGCTTACAGGCCGACTTTGGTCATGATGGCGAACa                                                                                                                  <  1:570452/150‑1 (MQ=255)
                        tCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACa                                                                                                                  <  1:752253/150‑1 (MQ=255)
                             accacTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCa                                                                                                             <  1:765269/146‑1 (MQ=255)
                               cGCTGCCCCCCCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGc                                                                                                           <  1:432664/150‑1 (MQ=255)
                               cGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGc                                                                                                           >  1:257858/1‑150 (MQ=255)
                                            aGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCgctg                                                                                              >  2:1017008/1‑150 (MQ=255)
                                             gggTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCgctgg                                                                                             <  1:709697/150‑1 (MQ=255)
                                              ggTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCgctggc                                                                                            <  2:1094045/150‑1 (MQ=255)
                                                     aGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGAcgcg                                                                                     >  1:755092/1‑150 (MQ=255)
                                                                aaCCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCAc                                                                          <  1:970166/150‑1 (MQ=255)
                                                                 aagaagaaTTCGCCGCCGAAATGCGCCGAGGCCCCGGACTGCACGCAGAGGCGCAGACCCCCGAAGCGCAGCACGGCGTACCGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGAAGCGTAAGCTGGAAcc                                                                         >  2:38588/8‑150 (MQ=255)
                                                                  ccAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCAcca                                                                        >  1:691557/1‑150 (MQ=255)
                                                                  ccAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCAcca                                                                        >  2:80372/1‑150 (MQ=255)
                                                                     gCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCAccagcc                                                                     <  1:1017008/150‑1 (MQ=255)
                                                                     gCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTATGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAGGCTGGCAccagcc                                                                     >  1:783260/1‑150 (MQ=255)
                                                                      ccATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCatgatg                                                                                                                         >  1:797536/1‑97 (MQ=255)
                                                                      ccATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCatgatg                                                                                                                         <  2:797536/97‑1 (MQ=255)
                                                                      ccATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCAccagcca                                                                    >  2:174912/1‑150 (MQ=255)
                                                                            gccgccGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCAc                                                                                                        <  1:449655/108‑1 (MQ=255)
                                                                            gccgccGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCAc                                                                                                        >  2:449655/1‑108 (MQ=255)
                                                                                               gCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGcc                                           <  2:691557/150‑1 (MQ=255)
                                                                                                ccccGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCt                                          <  2:719795/150‑1 (MQ=255)
                                                                                                 cccGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCtt                                         >  2:112742/1‑150 (MQ=255)
                                                                                                  ccGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTg                                        >  2:1097641/1‑150 (MQ=255)
                                                                                                          cccGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTGCGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGaaa                                <  2:30334/150‑1 (MQ=255)
                                                                                                             ggAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCa                                    >  2:99517/1‑143 (MQ=255)
                                                                                                             ggAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCa                                    <  1:99517/143‑1 (MQ=255)
                                                                                                                    agcagcGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAGGCCAt                          >  1:1025010/1‑146 (MQ=255)
                                                                                                                    agcagcGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAGGCCAt                          <  2:1025010/146‑1 (MQ=255)
                                                                                                                        gcGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCAcc                                                                                                       <  2:959899/65‑1 (MQ=255)
                                                                                                                        gcGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCAcc                                                                                                       >  1:959899/1‑65 (MQ=255)
                                                                                                                        gcGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAGGCCATGGCCAGGa                  <  2:755092/150‑1 (MQ=255)
                                                                                                                          gTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCgct                                                                                               <  2:117039/71‑1 (MQ=255)
                                                                                                                          gTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCgct                                                                                               >  1:117039/1‑71 (MQ=255)
                                                                                                                             cacaGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTg                                        <  2:1071973/123‑1 (MQ=255)
                                                                                                                             cacaGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTg                                        >  1:1071973/1‑123 (MQ=255)
                                                                                                                                  gcgcAGCACGGTGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGAcgc                                                                                      <  2:1034456/72‑1 (MQ=255)
                                                                                                                                  gcgcAGCACGGTGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGAcgc                                                                                      >  1:1034597/1‑72 (MQ=255)
                                                                                                                                  gcgcAGCACGGTGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGAcgc                                                                                      >  1:1034456/1‑72 (MQ=255)
                                                                                                                                  gcgcAGCACGGTGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGAcgc                                                                                      <  2:1034597/72‑1 (MQ=255)
                                                                                                                                       gcaCGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAGGCCATGGCCAGGATTGCCGCGCCGGCAt   <  1:967675/150‑1 (MQ=255)
                                                                                                                                        caCGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAgg                              >  1:1125012/1‑122 (MQ=255)
                                                                                                                                        caCGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAgg                              <  2:1125012/122‑1 (MQ=255)
                                                                                                                                        caCGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAACGTCTTGGCCAGGATTGCCGCGCCGGCATg  <  1:579744/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGCGTACAGGCCGACTTTGGTCATGATGGCGAACAGTGCAGCCACCGGCGCGCTGGCTGACGCGTAAGCTGGCACCAGCCAGAAGTTCAGCGGCCAGATACCGGCCTTGACCAGAAAGGCCATGGCCAGGATTGCCGCGCCGGCATG  >  NC_002947/2534901‑2535186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑