Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,488,064 T→A 16.5% K89M (AAG→ATG) ‡ degS ← serine endoprotease DegS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,488,0640TA16.5% 120.1 / 17.1 50K89M (AAG→ATG) ‡degSserine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base A (4/4);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

TTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTG  >  NC_002947/1487920‑1488212
                                                                                                                                                  |                                                                                                                                                    
ttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGTTTATCGCCTAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtcgtg                                                                                                                                                 <  2:185540/150‑1 (MQ=255)
  gttCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTggtggt                                                                                                                                               >  2:601199/3‑150 (MQ=255)
      agcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTAc                                                                                                                                           <  2:90831/150‑1 (MQ=255)
           gTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTTTAGCAGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCAAGGTGGTGAACaaaag                                                                                                                                      <  1:146595/150‑5 (MQ=255)
             aGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGAAGGTGAACGGCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTgg                                                                                                                                    <  1:677388/150‑1 (MQ=255)
                ccTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTACACCACCAAGgtggtgaac                                                                                                                                           >  2:309294/1‑137 (MQ=255)
                ccTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTACACCACCAAGGTGGTGAAc                                                                                                                                           <  1:309294/140‑1 (MQ=255)
                 ctctGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTACCGTCGTCGGCATGGTGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCAc                                                                                                                                >  2:635092/1‑150 (MQ=255)
                               ttACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCa                                                                                                                  >  1:606933/1‑150 (MQ=255)
                                tACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAg                                                                                                                 <  2:162375/150‑1 (MQ=255)
                                 aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                <  1:897031/150‑1 (MQ=255)
                                 aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                <  1:1133546/150‑1 (MQ=255)
                                     ccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCa                                                                                                            >  2:716609/1‑150 (MQ=255)
                                          cccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcg                                                                                                       >  1:747783/1‑150 (MQ=255)
                                          cccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcg                                                                                                       >  1:654834/1‑150 (MQ=255)
                                             aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACa                                                                                                                                          <  1:208363/112‑1 (MQ=255)
                                             aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACa                                                                                                                                          >  2:208363/1‑112 (MQ=255)
                                             aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcg                                                                                                    <  2:1187266/150‑1 (MQ=255)
                                               gCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTAGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcggc                                                                                                  <  1:489184/150‑1 (MQ=255)
                                                   gATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTAc                                                                                              >  1:825567/1‑150 (MQ=255)
                                                        ccAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                <  1:492706/127‑1 (MQ=255)
                                                        ccAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                >  2:492706/1‑127 (MQ=255)
                                                                   gTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                >  1:373343/1‑116 (MQ=255)
                                                                   gTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                <  2:373343/116‑1 (MQ=255)
                                                                                aTCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGt                                                                 <  2:271394/150‑1 (MQ=255)
                                                                                       aagaagCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGgc                                                          >  1:187491/1‑150 (MQ=255)
                                                                                                  gAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACGGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                               <  1:566518/150‑1 (MQ=255)
                                                                                                     ctccGGGT‑‑CTTCCAAGAGCGGGCGGGCGCTCCTCTCCACCACCCTGGCGGTCTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCGTGCATGATCCGGGTAGTCTGTGGCCCCTGCCGCAggct                                            <  2:39858/147‑3 (MQ=255)
                                                                                                     cTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGtt                                            <  2:747783/150‑1 (MQ=255)
                                                                                                      tGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTg                                           >  2:679445/1‑150 (MQ=255)
                                                                                                                ttCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctg                                   <  1:601199/150‑1 (MQ=255)
                                                                                                                ttCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctg                                   <  1:716609/150‑1 (MQ=255)
                                                                                                                    aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggctg                               <  2:606933/150‑1 (MQ=255)
                                                                                                                      cAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                             <  2:207368/150‑1 (MQ=255)
                                                                                                                                gcgcTTTTGTTCACCACCTTGGTGGTGTACAAGTTGGCCACCGCCGGTGCGTCCAGGGTCACGGCGTCGTCTTACGACCTTCATCTTCGGATGATGCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAcc                   >  1:390793/1‑150 (MQ=255)
                                                                                                                                    ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt               >  2:493810/1‑150 (MQ=255)
                                                                                                                                    ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt               >  1:763938/1‑150 (MQ=255)
                                                                                                                                       tgtacaccacCAAGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggct                                <  2:1092306/126‑1 (MQ=255)
                                                                                                                                       tGTACACCACCAAGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggct                                >  1:1092306/1‑130 (MQ=255)
                                                                                                                                       tGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCg            >  1:320490/1‑150 (MQ=255)
                                                                                                                                         ttCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGa          >  2:196687/1‑150 (MQ=255)
                                                                                                                                           caccacCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcgg                                                                                                   <  2:43408/59‑1 (MQ=37)
                                                                                                                                           caccacCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcgg                                                                                                   >  1:43408/1‑59 (MQ=38)
                                                                                                                                            accaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                               >  2:293055/1‑110 (MQ=255)
                                                                                                                                            accaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                               <  1:293055/110‑1 (MQ=255)
                                                                                                                                            accaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaa       <  1:196687/150‑1 (MQ=255)
                                                                                                                                             ccaccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCCGCCGGGGTCACGGCGTCGGCGTACGAAACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGaaaag      >  2:20439/1‑149 (MQ=255)
                                                                                                                                              caccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACg     >  1:1164587/1‑150 (MQ=255)
                                                                                                                                              caccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACg     >  1:745921/1‑150 (MQ=255)
                                                                                                                                                 cTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTg  <  2:1164587/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGT‑‑CTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTG  >  NC_002947/1487920‑1488212

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