Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,818,534 A→G 12.6% N212D (AAC→GAC)  PP_1622 → M23B subfamily metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,818,5340AG12.6% 132.7 / 8.6 54N212D (AAC→GAC) PP_1622M23B subfamily metallopeptidase
Reads supporting (aligned to +/- strand):  ref base A (25/21);  new base G (3/3);  total (29/25)
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.28e-01

ATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCACG  >  NC_002947/1818399‑1818682
                                                                                                                                       |                                                                                                                                                    
atCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCGTAGGcactga                                                                                                                                        <  1:998836/150‑6 (MQ=255)
      ggTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCgg                                                                                                                                  <  2:990336/150‑1 (MQ=255)
         gATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAg                                                                                                                               <  2:757891/150‑1 (MQ=255)
              gggACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGGGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTCCGGGACGGgcg                                                                                                                          >  1:92347/1‑149 (MQ=255)
                   aGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAgg                                                                                                                     >  2:426700/1‑150 (MQ=255)
                      cTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCa                                                                                                                  >  2:31663/1‑150 (MQ=255)
                            ttGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcag                                                                                                            >  1:1395824/1‑150 (MQ=255)
                             tgctgcGTCTGGTGGTGCGGTGGTTGATGCGGTGGGTGTTTTGGGGGGCTACGGCGAGCTGATCATCGTCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCGTAGCCGCCGGCTGTTGGTGCGGGAGGGGGGGGAGGTCAAggcagg                                                                                                           <  2:92347/150‑1 (MQ=255)
                             tgctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagg                                                                                                           <  2:1395824/150‑1 (MQ=255)
                                    tCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGGTATTACCGTAgg                                                                                                                                              >  1:1417528/1‑108 (MQ=255)
                                    tCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGGTATTACCGTAgg                                                                                                                                              <  2:1417528/108‑1 (MQ=255)
                                          ggTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCt                                                                                              <  1:862805/150‑1 (MQ=255)
                                           gTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCt                                                                                                                                            <  2:204941/103‑1 (MQ=255)
                                           gTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCt                                                                                                                                            >  1:204941/1‑103 (MQ=255)
                                               ggtggtTTATGCCGGGAGTGGATTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAAt                                                                                         >  2:854233/1‑150 (MQ=255)
                                                    tttATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACACGGACGCCTACGTCAGGGACAAGTTTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGt                                                                                    >  1:884703/1‑150 (MQ=255)
                                                    tttATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATAACGCATAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGt                                                                                    <  1:87342/150‑1 (MQ=255)
                                                            gggAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTATTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGc                                                                            >  1:333302/1‑150 (MQ=255)
                                                                            gCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGcagcag                                                                                                                      <  1:1037584/92‑1 (MQ=255)
                                                                            gCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGcagcag                                                                                                                      >  2:1037584/1‑92 (MQ=255)
                                                                             cTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAATGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGACGCt                                                           >  2:499774/1‑150 (MQ=255)
                                                                                      gCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGTGTGAAGCTGCATTTCga                                                  <  1:48077/150‑1 (MQ=255)
                                                                                             atcatcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAg                                                             >  1:1217674/1‑132 (MQ=255)
                                                                                             atcatcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAg                                                             <  2:1217674/132‑1 (MQ=255)
                                                                                              tcatcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGACATGGAGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcgcc                                          <  1:31663/150‑1 (MQ=255)
                                                                                                 tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATgg                                                                                       <  1:24424/102‑1 (MQ=255)
                                                                                                 tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATgg                                                                                       >  2:24424/1‑102 (MQ=255)
                                                                                                 tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCCGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                             >  2:577199/1‑144 (MQ=255)
                                                                                                 tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCCGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                             <  1:577199/144‑1 (MQ=255)
                                                                                                   aaGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAg                                     <  1:854233/150‑1 (MQ=255)
                                                                                                     gCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATTGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                                   >  2:1140215/1‑150 (MQ=255)
                                                                                                       acaGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                             <  2:526889/138‑1 (MQ=255)
                                                                                                       acaGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                             >  1:526889/1‑138 (MQ=255)
                                                                                                           cGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAcc                             >  1:215422/1‑150 (MQ=255)
                                                                                                              cacCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGt                          <  1:426700/150‑1 (MQ=255)
                                                                                                                 caaCGTCACTGTCTAACTTCAGAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGa                       >  2:59905/3‑150 (MQ=255)
                                                                                                                  tACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGAt                      >  1:36156/1‑150 (MQ=255)
                                                                                                                  tACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGAt                      >  2:442630/1‑150 (MQ=255)
                                                                                                                  tACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTAGCAGTAAGGGAAAAACCGTCGAt                      >  2:1149243/1‑150 (MQ=255)
                                                                                                                      tCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAggg                                                                                                                             >  2:1365613/1‑43 (MQ=255)
                                                                                                                      tCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAggg                                                                                                                             <  1:1365613/43‑1 (MQ=255)
                                                                                                                      caacagCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCAc                  >  2:737631/6‑150 (MQ=255)
                                                                                                                        acagCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAccc                            <  2:1363746/135‑1 (MQ=255)
                                                                                                                        acagCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAccc                            >  1:1363746/4‑138 (MQ=255)
                                                                                                                           gCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCa              <  1:628913/149‑1 (MQ=255)
                                                                                                                           gCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCa              >  2:628913/1‑149 (MQ=255)
                                                                                                                                cGGTCATAACCGCAGGCTGTTGGTTAGGGAGGGGCAGCAGGGAAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTATGGGCACAGATTGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCt        <  2:215422/150‑1 (MQ=255)
                                                                                                                                   tCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggca                                                                                                             >  1:137132/1‑46 (MQ=255)
                                                                                                                                   tCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggca                                                                                                             <  2:137132/46‑1 (MQ=255)
                                                                                                                                   tCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCgg                                                                   >  2:370040/1‑88 (MQ=255)
                                                                                                                                   tCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCgg                                                                   <  1:370040/88‑1 (MQ=255)
                                                                                                                                     aTAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTAc                                                                                >  2:462833/1‑73 (MQ=255)
                                                                                                                                     aTAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTAc                                                                                <  1:462833/73‑1 (MQ=255)
                                                                                                                                      tAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCAcg  >  1:641819/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                    
ATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCACG  >  NC_002947/1818399‑1818682

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