Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,546,490 C→T 8.6% P531S (CCG→TCG)  PP_2238 → metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,546,4900CT 8.6% 157.3 / 4.7 58P531S (CCG→TCG) PP_2238metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (23/30);  new base T (2/3);  total (25/33)
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.62e-01

TCGCCGTTGCTCGACGCCATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCT  >  NC_002947/2546349‑2546637
                                                                                                                                              |                                                                                                                                                   
tCGCCGTTGCTCGACGCCATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAATTCATCGCCTACCAGGCGGCACATGCACCGAAGAc                                                                                                                                              >  4:87454/1‑149 (MQ=255)
  gCCGTTGCTCGACGCCATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACgg                                                                                                                                            <  3:187663/149‑1 (MQ=255)
  gCCGTTGCTCGACGCCATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACgg                                                                                                                                            >  4:187663/1‑149 (MQ=255)
   ccGTTGCTCGACGCCATTGCCCGGGCGCACGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACgg                                                                                                                                            >  1:196600/1‑149 (MQ=255)
         cTCGACGCCATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCATCGCCGTCTGCGGTGCa                                                                                                                                     <  3:274182/149‑1 (MQ=255)
             aCGCCATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTg                                                                                                                                 <  7:282052/149‑1 (MQ=255)
                ccATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGc                                                                                                                              <  2:143867/149‑1 (MQ=255)
                   ttGCCCGGGCGC‑CGATTTGGGATGAACAGGTCGACTATGGCTATGGTACGGGCCACGGCGTGGGATACTTAATGAACGTGCATGAGGGGACTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGc                                                                                                                                      <  4:144875/138‑1 (MQ=255)
                   ttGCCCGGGCGC‑CGATCTGGGCTGACCAGGGCGACTATGGCCATGGAAGCGGCCACTGAGTGGGCTACATCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCAACCAGGCGGCACCTGCACCGCAGACGGCGATGc                                                                                                                                      >  3:144875/1‑138 (MQ=255)
                       ccGGGCGC‑CGTTCTTGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCa                                                                                                                       <  3:87454/149‑1 (MQ=255)
                         gggCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGGGGGCTACTTCATGAACGTGCATGAGGGGGCTCAAGTGATCGCCTACCAGGCGGCACCTGCCCCGCAGACGGCGATGCAAGTGGGCTTGCTCCGc                                                                                                                     >  4:260412/1‑149 (MQ=255)
                                cGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCtcgatc                                                                                                               >  5:327061/1‑149 (MQ=255)
                                      gggCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCtcga                                                                                                                 <  2:282914/141‑1 (MQ=255)
                                      gggCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGCCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCtcga                                                                                                                 >  1:282914/1‑141 (MQ=255)
                                       ggCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCtcgatc                                                                                                               <  6:327061/142‑1 (MQ=255)
                                       ggCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCtcgatc                                                                                                               <  8:185598/142‑1 (MQ=255)
                                       ggCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGcc                                                                                                          <  1:104159/147‑1 (MQ=255)
                                       ggCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCgg                                                                                                        <  8:181576/149‑1 (MQ=255)
                                       ggCGGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCgg                                                                                                        <  2:160610/149‑1 (MQ=255)
                                         cTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGgt                                                                                                      <  1:146094/149‑1 (MQ=255)
                                              gaGGTCGCCTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGgtacgta                                                                                                 <  4:264771/148‑1 (MQ=255)
                                                 gTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGCCGGCGAGGCTAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCg                                                                                              >  2:54284/1‑149 (MQ=255)
                                                 gTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCg                                                                                              >  1:140111/1‑149 (MQ=255)
                                                     aCTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCg                                                                                          >  3:36001/1‑149 (MQ=255)
                                                       tATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCggg                                                                                        >  3:40518/1‑149 (MQ=255)
                                                         tGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCgggg                                                                                       <  4:36001/148‑1 (MQ=255)
                                                           gCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCtg                                                                                    >  6:181501/1‑149 (MQ=255)
                                                            ccATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCtgt                                                                                   <  2:180916/149‑1 (MQ=255)
                                                               tGGTACCGGCCACGGCGTGGGCTACTTCATGACCGTGCATGAGGGGCCTCAAGTGCTCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATTATGAGCTCGATCGAGCCGGGTACTTACCGCCCGGGGCTGTggg                                                                                <  7:307462/149‑1 (MQ=255)
                                                                   aCCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCATCGCCGTCTGCGGTGCAGGTgc                                                                                                                                <  8:222644/97‑2 (MQ=38)
                                                                   aCCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCATCGCCGTCTGCGGTGCAGGTgc                                                                                                                                >  7:222644/1‑96 (MQ=38)
                                                                                tGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAAc                                                               <  5:181501/149‑1 (MQ=255)
                                                                                            tGAACGTGCATGACGGGCCTCAAGGGATCGCCTACCCGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAAc                                                   <  3:260412/149‑1 (MQ=255)
                                                                                             gAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCg                                                                    >  5:226364/1‑131 (MQ=255)
                                                                                             gAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCg                                                                    <  6:226364/131‑1 (MQ=255)
                                                                                              aaCGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCtgt                                                                                   <  8:313241/115‑1 (MQ=255)
                                                                                              aaCGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCtgt                                                                                   >  7:313241/1‑115 (MQ=255)
                                                                                              aaCGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAAcc                                                  >  3:222373/1‑148 (MQ=255)
                                                                                                cGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCgag                                               <  4:14775/149‑1 (MQ=255)
                                                                                                      tGAGGGGCCTCAAGTGATCGCCTACCCGGCGGCACCTGCACCGCAGACGGCGATGCAAGTCGGCATGATCAGCTCGATCGAGCCGGGTCCGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCgg                                         <  1:54284/149‑1 (MQ=255)
                                                                                                      tGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCgg                                         >  1:224763/1‑149 (MQ=255)
                                                                                                        aGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCggg                                        <  3:18037/148‑1 (MQ=255)
                                                                                                             ccTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAagag                                  >  4:35515/1‑149 (MQ=255)
                                                                                                                       tCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCtt                            <  1:320881/145‑1 (MQ=255)
                                                                                                                       tCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCtt                            >  2:320881/1‑145 (MQ=255)
                                                                                                                       tCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGt                        <  3:35515/149‑1 (MQ=255)
                                                                                                                       gcGCCTACCACGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATCATCAGCTCGAGCGCGCCGGGTACGGACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGt                        <  4:222373/148‑1 (MQ=255)
                                                                                                                        cGCCTACCAGGCGGCACCTGCACCGCAGACGTCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGTCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGt                        <  8:136983/148‑1 (MQ=255)
                                                                                                                         gCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGTGa                      >  1:192517/1‑149 (MQ=255)
                                                                                                                          ccTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTg                          >  5:41337/1‑144 (MQ=255)
                                                                                                                          ccTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTg                          <  6:41337/144‑1 (MQ=255)
                                                                                                                            tACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCtggtg                                                         <  4:191262/111‑1 (MQ=255)
                                                                                                                            tACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCtggtg                                                         >  3:191262/1‑111 (MQ=255)
                                                                                                                               cAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTg                          >  8:315715/1‑139 (MQ=255)
                                                                                                                               cAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTg                          <  7:315715/139‑1 (MQ=255)
                                                                                                                                      cACTTGCATCGCCGTCTGCGGTGCAGGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGt                        >  2:122044/1‑134 (MQ=255)
                                                                                                                                      cACTTGCATCGCCGTCTGCGGTGCAGGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGt                        <  1:122044/134‑1 (MQ=255)
                                                                                                                                             aCCGCAGACGGCGATGCAAGTGGGCATGATCAGATCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCt  >  2:145183/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCGCCGTTGCTCGACGCCATTGCCCGGGCGC‑CGATCTGGGCTGACCAGGTCGACTATGGCCATGGTACCGGCCACGGCGTGGGCTACTTCATGAACGTGCATGAGGGGCCTCAAGTGATCGCCTACCAGGCGGCACCTGCACCGCAGACGGCGATGCAAGTGGGCATGATCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAGAGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCT  >  NC_002947/2546349‑2546637

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