Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,488,049 T→A 12.8% N92Y (AAC→TAC)  degS ← serine endoprotease DegS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,488,0490TA12.8% 118.9 / 8.9 47N92Y (AAC→TAC) degSserine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base A (3/3);  total (23/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.05e-01

GACATGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCG  >  NC_002947/1487903‑1488195
                                                                                                                                                  |                                                                                                                                                  
gACATGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTca                                                                                                                                                  <  2:228670/149‑1 (MQ=255)
   aTGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTcacca                                                                                                                                               >  5:205677/1‑149 (MQ=255)
    tGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGtt                                                                                                                                                    <  8:72051/143‑1 (MQ=255)
    tGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGtt                                                                                                                                                    >  7:72051/1‑143 (MQ=255)
    tGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTcaccac                                                                                                                                              >  8:64846/1‑149 (MQ=255)
    tGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTcaccac                                                                                                                                              <  4:28374/149‑1 (MQ=255)
          ttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtgg                                                                                                                                         <  1:59663/148‑1 (MQ=255)
           tgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggtg                                                                                                                                       >  3:138791/1‑149 (MQ=255)
              tCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTcac                                                                                                                                                 >  5:290194/1‑136 (MQ=255)
              tCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTcac                                                                                                                                                 <  6:290194/136‑1 (MQ=255)
              tCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGgtggtg                                                                                                                                    <  5:145031/149‑1 (MQ=255)
                cgcagGTAGCCCTCTGGGCTCATGATTACCGCCGACCCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTACACCACCAAGgtggtgaa                                                                                                                                  <  8:277934/148‑3 (MQ=255)
                agcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGgtgt                                                                                                                                   <  4:138791/148‑1 (MQ=255)
                 gcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTACACCACCAAGGTGGTGAAc                                                                                                                                 <  6:3072/149‑1 (MQ=255)
                                ggCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggt                                                                                                                                        >  2:83235/1‑127 (MQ=255)
                                ggCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggt                                                                                                                                        <  1:83235/127‑1 (MQ=255)
                                     aTGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCggg                                                                                                               <  6:205677/147‑1 (MQ=255)
                                          tACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggt                                                                                                                                        <  7:279126/117‑1 (MQ=255)
                                          tACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggt                                                                                                                                        >  8:279126/1‑117 (MQ=255)
                                             cgccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAggg                                                                                                     >  3:11277/1‑149 (MQ=255)
                                               ccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTAGCGCCGAAGCAGCGGGGGAACGGCGGGGGGTCGAGCCGCGGGGGGGGGGTTTTGTTCACCCCCGGGGGGGGGTACGGGTTGGCCACCGCGGGGGGGGCCAGGGTc                                                                                                   >  8:73590/1‑149 (MQ=255)
                                                    cccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTggtggt                                                                                                                                     <  2:203763/110‑1 (MQ=255)
                                                       aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtc                                                                                           >  6:94755/1‑149 (MQ=255)
                                                         gCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcgg                                                                                         <  8:2939/149‑1 (MQ=255)
                                                                   cAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACAGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGacac                                                                               >  3:193315/1‑149 (MQ=255)
                                                                       ggcgCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCggg                                                                           <  4:193315/149‑1 (MQ=255)
                                                                                 gggCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGa                                                                 >  3:236063/1‑149 (MQ=255)
                                                                                            cGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTc                                                      >  8:1052/1‑149 (MQ=255)
                                                                                             gCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCt                                                     >  8:200623/1‑149 (MQ=255)
                                                                                             gCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTCGCATGATCCGGGTAGTc                                                      <  8:248269/148‑1 (MQ=255)
                                                                                              ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACa                                                                                                                                >  5:298188/1‑73 (MQ=255)
                                                                                              ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACa                                                                                                                                <  6:298188/73‑1 (MQ=255)
                                                                                              ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCtg                                                    <  7:200623/149‑1 (MQ=255)
                                                                                                        ggcggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtgc                                          <  4:236063/149‑1 (MQ=255)
                                                                                                                  cGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGtt                                  <  7:1052/147‑1 (MQ=255)
                                                                                                                   gggTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGAc                               <  7:64846/149‑1 (MQ=255)
                                                                                                                     gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGt                                                       <  4:85423/123‑1 (MQ=255)
                                                                                                                     gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACAGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGt                                                       >  3:85423/1‑123 (MQ=255)
                                                                                                                        ttCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCct                          >  6:235200/1‑149 (MQ=255)
                                                                                                                            aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAgg                                    <  3:78978/135‑1 (MQ=255)
                                                                                                                            aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAgg                                    >  4:78978/1‑135 (MQ=255)
                                                                                                                            aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggct                      >  6:32105/1‑149 (MQ=255)
                                                                                                                             aCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGAGCCGGGTAGTCTGTGGCGCTTGCTGCTGGTTGACGTCctggctg                     >  7:273802/1‑149 (MQ=255)
                                                                                                                              cAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggctgg                    <  3:188230/149‑1 (MQ=255)
                                                                                                                                         ccaacctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                   >  6:584/7‑139 (MQ=255)
                                                                                                                                         ccaacctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                   <  5:584/133‑1 (MQ=255)
                                                                                                                                            ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGTGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAt      >  3:121483/1‑149 (MQ=255)
                                                                                                                                             cctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGTGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt     >  3:70992/3‑149 (MQ=255)
                                                                                                                                              ctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTCGCGCTTGCTCCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTg    >  8:33538/2‑149 (MQ=255)
                                                                                                                                                 ttCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCg  >  5:64481/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GACATGGTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCG  >  NC_002947/1487903‑1488195

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