Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 160,882 T→A 24.1% S392T (TCC→ACC)  degP → serine endoprotease (protease Do), membrane‑associated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE160,8820TA24.1% 188.0 / 56.0 116S392T (TCC→ACC) degPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base T (42/46);  new base A (21/7);  total (63/53)
Fisher's exact test for biased strand distribution p-value = 1.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

ACGGTAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGAAC  >  minE/160817‑160950
                                                                 |                                                                    
aCGGTAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCa                                                                   <  1:2501991/69‑1 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:1036479/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:3429515/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              <  1:1726796/69‑1 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:374687/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:3407888/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:1776039/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:2514281/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:2583136/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:525720/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:799904/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:1382842/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:1300703/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:1273979/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:3249031/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              >  1:3121424/1‑69 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              <  1:2356259/69‑1 (MQ=255)
     aaGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTc                                                              <  1:213581/69‑1 (MQ=255)
             tAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACgg                                                         <  1:1323993/66‑1 (MQ=255)
             tAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCa                                                       >  1:1495623/1‑68 (MQ=255)
             tAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCa                                                       >  1:268341/1‑68 (MQ=255)
                cGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTg                                                    <  1:1351210/68‑1 (MQ=255)
                cGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTg                                                    <  1:1273237/68‑1 (MQ=255)
                cGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTg                                                    <  1:1224268/68‑1 (MQ=255)
                cGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTg                                                    <  1:1598804/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCATCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:2604728/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:1553795/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:355802/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:1602469/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:422543/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:1292418/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:2760972/68‑1 (MQ=255)
                 gTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:2765122/68‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCTGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:854702/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:2859896/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:292468/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:3290584/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:2926489/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:2702702/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:2265567/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:2200653/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:218877/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:3158380/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:904902/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:857387/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:4824/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:3294450/69‑1 (MQ=255)
                         ggAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:1989808/69‑1 (MQ=255)
                         ggAACTGCAGCAAAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTga                                          <  1:2164716/69‑1 (MQ=255)
                            aCTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTg                                                    <  1:200517/56‑1 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   >  1:2827686/1‑54 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:32493/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:1937849/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:2454325/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:2932878/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:702642/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:1209077/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:201015/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:903784/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:2335080/1‑68 (MQ=255)
                               gcagcaGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGa                                     >  1:3137850/1‑68 (MQ=255)
                                   cagagcagCCAGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaac                                <  1:917214/69‑1 (MQ=255)
                                   cagagcagCCAGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaac                                <  1:2851299/69‑1 (MQ=255)
                                   cagagcagCCAGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaac                                <  1:2767419/69‑1 (MQ=255)
                                   cagagcagCCAGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaac                                <  1:1456396/69‑1 (MQ=255)
                                   cagagcagCCAGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaac                                <  1:1901028/69‑1 (MQ=255)
                                   cagagcagCCAGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaac                                <  1:2436764/69‑1 (MQ=255)
                                       gcagcCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                   <  1:1044867/46‑1 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCATCGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:1954852/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:3290966/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:588618/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:522312/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:514492/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:1402638/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:1706736/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:1844989/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:1054854/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:1939565/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:322917/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:1160290/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggta                         >  1:3141918/1‑64 (MQ=255)
                                             aGAATCAGGTTGATTCCAGCACCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAggt                          >  1:173603/1‑64 (MQ=255)
                                                aTCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATc                   >  1:610273/1‑69 (MQ=255)
                                                            ccAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGt         >  1:2204008/1‑67 (MQ=255)
                                                            ccAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTa        >  1:2500109/1‑68 (MQ=255)
                                                            ccAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTa        >  1:2054457/1‑68 (MQ=255)
                                                            ccAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTa        >  1:1416576/1‑68 (MQ=255)
                                                            ccAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTa        >  1:1124251/1‑68 (MQ=255)
                                                               gCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTg     <  1:1555520/68‑1 (MQ=255)
                                                               gCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTg     <  1:2578669/68‑1 (MQ=255)
                                                                 tCCGTCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:3110737/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:1035110/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:1468435/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:1604731/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:1672086/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:3189449/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:2393834/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:3057858/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:2528418/1‑69 (MQ=255)
                                                                 tCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGaac  >  1:2722977/1‑69 (MQ=255)
                                                                 |                                                                    
ACGGTAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTGAAC  >  minE/160817‑160950

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