Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,630,233 C→T 7.2% P216S (CCC→TCC)  PP_1430 → DegP‑like serine endoprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,630,2330CT 7.2% 159.2 / 2.1 55P216S (CCC→TCC) PP_1430DegP‑like serine endoprotease
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base T (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

AAGGTCGATGGCAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGTTTCATGGGCCTGT  >  NC_002947/1630086‑1630377
                                                                                                                                                   |                                                                                                                                                
aaGGTCGATGGCAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGcc                                                                                                                                                 >  1:186808/1‑149 (MQ=255)
    tCGATGGCAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAAt                                                                                                                                             <  6:58891/149‑1 (MQ=255)
    tCGATGGCAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAAt                                                                                                                                             >  8:195313/1‑149 (MQ=255)
    tCGATGGCAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAAt                                                                                                                                             <  2:186808/149‑1 (MQ=255)
     cGATGGCAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTCCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAAt                                                                                                                                             >  2:264587/1‑148 (MQ=255)
           cAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGCCAcc                                                                                                                                       >  5:88060/1‑148 (MQ=255)
            aaGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGac                                                                                                                                          <  6:122052/144‑1 (MQ=255)
            aaGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGac                                                                                                                                          >  5:122052/1‑144 (MQ=255)
                         cGGTGAAACTGGGTGACTCCGAGAAGCTGAAGGGGGGGGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTGCGACCATTGGGTGTCCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATc                                                                                                                        <  6:230093/149‑1 (MQ=255)
                         cGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTGCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATc                                                                                                                        <  2:82642/149‑1 (MQ=255)
                                  tGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGAc                                                                                                               >  8:2500/1‑149 (MQ=255)
                                          tCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGTCCAATGACACCTACGTACCGTTCATc                                                                                                                        <  3:203451/132‑1 (MQ=255)
                                          tCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGTCCAATGACACCTACGTACCGTTCATc                                                                                                                        >  4:203451/1‑132 (MQ=255)
                                            cGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGTCCATc                                                                                                      <  7:195313/148‑1 (MQ=255)
                                             gagaAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAACGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCaa                                                                                                    >  8:160806/1‑149 (MQ=255)
                                                 aGCTCAAGGTGGGCTAGTGGGTACTGGCCATCGGTTCGCGGTTCGGCTTCGACCATTCGGTGACCAAAGGTAGCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAAccc                                                                                                 <  7:240927/148‑1 (MQ=255)
                                                  gCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGa                                                                                                                    >  8:227889/1‑128 (MQ=255)
                                                  gCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGa                                                                                                                    <  7:227889/128‑1 (MQ=255)
                                                  gCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTCCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCgg                                                                                               <  1:264587/149‑1 (MQ=255)
                                                    tCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACg                                                                                                              <  7:37654/132‑1 (MQ=255)
                                                    tCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACg                                                                                                              >  8:37654/1‑132 (MQ=255)
                                                      aaGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGt                                                                                             <  7:2500/147‑1 (MQ=255)
                                                       aGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTAc                                                                                                                                    >  2:208609/1‑107 (MQ=255)
                                                       aGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTAc                                                                                                                                    <  1:208609/107‑1 (MQ=255)
                                                                aGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCgg                                                                                               >  7:217550/1‑135 (MQ=255)
                                                                aGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCgg                                                                                               <  8:217550/135‑1 (MQ=255)
                                                                 gTGGTTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCcggcg                                                                                   <  3:191153/146‑1 (MQ=255)
                                                                        cTGGCCATCGGATCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGtt                                                                         <  5:163645/149‑1 (MQ=255)
                                                                          ggCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGCCCATCAACCCGGGAACCACCGGCGGCCCGCTGtta                                                                        >  6:97029/1‑147 (MQ=255)
                                                                              aTCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGACGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATAAACCCGGGTAACTCCGGCGGCCCGCTGTTCAAc                                                                     >  1:153037/1‑147 (MQ=255)
                                                                                      gCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCa                                                                       >  1:2352/1‑137 (MQ=255)
                                                                                      gCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCa                                                                       <  2:2352/137‑1 (MQ=255)
                                                                                            gggCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAgtg                                                      <  5:97029/147‑1 (MQ=255)
                                                                                                   gACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGa                                               >  8:27975/1‑148 (MQ=255)
                                                                                                      cATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGTCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCaa                                           >  8:13443/1‑149 (MQ=255)
                                                                                                       aTTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAAc                                          >  3:230276/1‑149 (MQ=255)
                                                                                                       aTTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAAc                                          >  4:13641/1‑149 (MQ=255)
                                                                                                          cGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTc                                                                        >  5:93954/1‑116 (MQ=255)
                                                                                                          cGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTc                                                                        <  6:93954/116‑1 (MQ=255)
                                                                                                           ggTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGc                                      <  7:27975/149‑1 (MQ=255)
                                                                                                                 cAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCt                                >  8:30886/1‑149 (MQ=255)
                                                                                                                           gTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGtt                       <  3:13641/148‑1 (MQ=255)
                                                                                                                            tCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGtggtgg                                                  <  7:196096/120‑1 (MQ=255)
                                                                                                                            tCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGtggtgg                                                  >  8:196096/1‑120 (MQ=255)
                                                                                                                            tCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCa                                     >  7:160540/1‑133 (MQ=255)
                                                                                                                            tCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCa                                     <  8:160540/133‑1 (MQ=255)
                                                                                                                              aGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTgg                   <  2:235025/149‑1 (MQ=255)
                                                                                                                                cgcCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCg                 <  7:160806/149‑1 (MQ=255)
                                                                                                                                   cAAGGGCCGGACCCTGTCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGtt              <  7:13443/149‑1 (MQ=255)
                                                                                                                                   cAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCACCCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGtt              >  7:290593/1‑149 (MQ=255)
                                                                                                                                     aGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTATGGCGGttt             >  2:54861/1‑148 (MQ=255)
                                                                                                                                       ggCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTAAACATGAAAGGCGAAGTGGTGGGGATTAACTCGCTGATCTTCACCAGTTCTGGCGGTTTAAt          >  8:116999/1‑149 (MQ=255)
                                                                                                                                         ccGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGTTTCATgg        >  6:98321/1‑149 (MQ=255)
                                                                                                                                         ccGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGTTTCATgg        <  6:88060/149‑1 (MQ=255)
                                                                                                                                         ccGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGTTTCATgg        <  5:98321/149‑1 (MQ=255)
                                                                                                                                               ccTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGTTTCATGGGCCTGt  >  4:294437/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AAGGTCGATGGCAAGAATTTGCCAACGGTGAAACTGGGTGACTCCGAGAAGCTCAAGGTGGGCGAGTGGGTACTGGCCATCGGTTCGCCGTTCGGCTTCGACCATTCGGTGACCAAAGGTATCGTCAGCGCCAAGGGCCGGACCCTGCCCAATGACACCTACGTACCGTTCATCCAGACCGACGTGGCCATCAACCCGGGTAACTCCGGCGGCCCGCTGTTCAACATGAAAGGCGAAGTGGTGGGGATCAACTCGCAGATCTTCACCCGTTCTGGCGGTTTCATGGGCCTGT  >  NC_002947/1630086‑1630377

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: ‑