Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,818,531 C→T 15.7% H211Y (CAT→TAT)  PP_1622 → M23B subfamily metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,818,5310CT15.7% 116.2 / 13.9 51H211Y (CAT→TAT) PP_1622M23B subfamily metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base T (4/4);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCT  >  NC_002947/1818382‑1818676
                                                                                                                                                     |                                                                                                                                                 
tGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTc                                                                                                                                                   <  1:500405/150‑1 (MQ=255)
                 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                                                       <  2:1217674/132‑1 (MQ=255)
                                                                                                              atcatcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAg                                                       >  1:1217674/1‑132 (MQ=255)
                                                                                                               tcatcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGACATGGAGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcgcc                                    <  1:31663/150‑1 (MQ=255)
                                                                                                                  tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATgg                                                                                 <  1:24424/102‑1 (MQ=255)
                                                                                                                  tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATgg                                                                                 >  2:24424/1‑102 (MQ=255)
                                                                                                                  tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCCGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                       <  1:577199/144‑1 (MQ=255)
                                                                                                                  tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCCGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                       >  2:577199/1‑144 (MQ=255)
                                                                                                                    aaGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAg                               <  1:854233/150‑1 (MQ=255)
                                                                                                                      gCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATTGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                             >  2:1140215/1‑150 (MQ=255)
                                                                                                                        acaGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                       >  1:526889/1‑138 (MQ=255)
                                                                                                                        acaGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                       <  2:526889/138‑1 (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                                                                                                                       <  1:1365613/43‑1 (MQ=255)
                                                                                                                                       tCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAggg                                                                                                                       >  2:1365613/1‑43 (MQ=255)
                                                                                                                                       caacagCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCAc            >  2:737631/6‑150 (MQ=255)
                                                                                                                                         acagCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAccc                      >  1:1363746/4‑138 (MQ=255)
                                                                                                                                         acagCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAccc                      <  2:1363746/135‑1 (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)
                                                                                                                                                     |                                                                                                                                                 
TGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCT  >  NC_002947/1818382‑1818676

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