Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,908,653 G→A 7.5% S14S (AGC→AGT)  hasE ← secretion metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,908,6530GA 7.5% 152.3 / 2.5 53S14S (AGC→AGT) hasEsecretion metalloprotease
Reads supporting (aligned to +/- strand):  ref base G (27/22);  new base A (2/2);  total (29/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 = 7.35e-01

GTGACCTTGACCGTACCGGTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCGCGGCGAACACTTCGTCGCGGCTG  >  NC_002947/2908510‑2908794
                                                                                                                                               |                                                                                                                                             
gTGACCTTGACCGTACCGGTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCtctc                                                                                                                                          >  5:360160/1‑149 (MQ=255)
     cTTGACCGTACCGGTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTcggcgg                                                                                                                                     >  3:349985/1‑149 (MQ=255)
      ttGACCGTACCGGTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCGAGATCGATCAGGTGGCTGCTGTCGGCggg                                                                                                                                    >  8:97609/1‑149 (MQ=255)
      ttGACCGTACCGGTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCggg                                                                                                                                    >  8:147552/1‑149 (MQ=255)
               ccGGTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTc                                                                                                                           <  4:349985/149‑1 (MQ=255)
                 ggTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCgg                                                                                                                         <  3:271206/149‑1 (MQ=255)
                     gCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATg                                                                                                                     >  6:357016/1‑149 (MQ=255)
                       gACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCgggg                                                                                                                                   <  5:331546/133‑1 (MQ=255)
                       gACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCAGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCgggg                                                                                                                                   >  6:331546/1‑133 (MQ=255)
                          cACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAg                                                                                                                <  7:147552/149‑1 (MQ=255)
                          cACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAg                                                                                                                >  2:49847/1‑149 (MQ=255)
                           aCGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGt                                                                                                               <  8:211101/149‑1 (MQ=255)
                               cGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATc                                                                                                           <  6:157536/149‑1 (MQ=255)
                                 gCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTg                                                                                                         >  7:186960/1‑149 (MQ=255)
                                   gTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCa                                                                                                       >  6:127001/1‑149 (MQ=255)
                                       aGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGt                                                                                                   >  1:305156/1‑149 (MQ=255)
                                        gTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTc                                                                                                  >  4:379364/1‑149 (MQ=255)
                                        gTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTc                                                                                                  >  4:252071/1‑149 (MQ=255)
                                           gcgcCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGAGCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGt                                                                                               <  8:97881/149‑1 (MQ=255)
                                              ccAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTAc                                                                                            <  3:252071/149‑1 (MQ=255)
                                                aggcaggACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTc                                                                                          <  1:49847/149‑1 (MQ=255)
                                                   caggACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCg                                                                                       <  7:85587/149‑1 (MQ=255)
                                                      gACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGtcc                                                                                    >  3:28092/1‑147 (MQ=255)
                                                        caacaGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACtt                                                                                  >  2:147983/3‑149 (MQ=255)
                                                         cacaGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACttt                                                                                 <  2:305156/149‑1 (MQ=255)
                                                                ggAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGaa                                                                          <  8:186960/149‑1 (MQ=255)
                                                                 gAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAg                                                                         >  7:126805/1‑149 (MQ=255)
                                                                   aaCCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGc                                                                        >  6:94249/1‑148 (MQ=255)
                                                                                       aCCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGc                                                   >  3:176555/1‑149 (MQ=255)
                                                                                        ccAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAg                                                          <  6:154470/141‑1 (MQ=255)
                                                                                        ccAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAg                                                          >  5:154470/1‑141 (MQ=255)
                                                                                        ccAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGcc                                                  <  4:220686/149‑1 (MQ=255)
                                                                                          aGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCgg                                                                                                <  4:290966/101‑1 (MQ=255)
                                                                                          aGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCgg                                                                                                >  3:290966/1‑101 (MQ=255)
                                                                                           gCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCgg                                                >  1:337965/1‑148 (MQ=255)
                                                                                           gCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTACTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCgg                                                >  3:284832/1‑148 (MQ=255)
                                                                                           gCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTACTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCgg                                                >  3:289878/1‑148 (MQ=255)
                                                                                                     ccAGCGTACAGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGGTCTCGGCGGCGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCCGTGCCGCAGCCGTTGCCGGAATCACCGCgg                                     >  8:14831/1‑149 (MQ=255)
                                                                                                      cAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGc                                    >  6:372730/1‑149 (MQ=255)
                                                                                                        gCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCg                                                             >  1:334521/1‑122 (MQ=255)
                                                                                                        gCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCg                                                             <  2:334521/122‑1 (MQ=255)
                                                                                                         cGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCtt                                 <  1:154769/149‑1 (MQ=255)
                                                                                                           tACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTgc                               <  5:127001/149‑1 (MQ=255)
                                                                                                              cGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTACTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCa                            <  4:289878/149‑1 (MQ=255)
                                                                                                              cGGGCGGCCCGCGTCCAGATCGATCAGGGGGCTACTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCa                            <  4:284832/149‑1 (MQ=255)
                                                                                                                    gcCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCgcgg                      <  6:360160/149‑1 (MQ=255)
                                                                                                                      ccGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCgcggcg                    <  4:28092/149‑1 (MQ=255)
                                                                                                                       cgagTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGCCCGAAGCCAGCGGTGCCGCAGCCCTTGCCGGAATCCCCGCGGCCTTCCGCAGCgctgcg                    >  8:141440/4‑148 (MQ=255)
                                                                                                                              aGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCGCGGCGAACACTtc            >  6:143350/1‑149 (MQ=255)
                                                                                                                                   gatcAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCACAGCGCGGCGAACACTTCGTcgc       >  4:214795/1‑149 (MQ=255)
                                                                                                                                     tcAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCGCGGCGAACACTTCGTCGCgg     >  6:151661/1‑149 (MQ=255)
                                                                                                                                       aGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCGCGGCGAACACTTCGTCGCGGCt   <  2:337965/149‑1 (MQ=255)
                                                                                                                                        ggTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCGCGGCGAACACTTCGTCGCGGCTg  <  1:147983/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GTGACCTTGACCGTACCGGTGGCGACCACGCCGGCGTCCAGTGGCGCCAGGCAGGACCACAGCAGGAAACCGCCGAAGCCGGCGATGACCAGCCATACGCCCCAGCGTACCGGGCGGCCCGCGTCCAGATCGATCAGGTGGCTGCTCTCGGCGGGGATCAGTTCGGCATGTTGAGTCATCTGCATGGTCGGTTACTCCCGTACTTTGACCGAAGCCAGCGGTGCCGCAGCCGTTGCCGGAATCACCGCGGCCTTGCGCAGCGCGGCGAACACTTCGTCGCGGCTG  >  NC_002947/2908510‑2908794

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