Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,197,282 G→A 7.3% D420N (GAT→AAT)  degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,197,2820GA 7.3% 398.0 / 13.9 123D420N (GAT→AAT) degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (77/37);  new base A (3/6);  total (80/43)
Fisher's exact test for biased strand distribution p-value = 6.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

AGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAA  >  Exported/2197212‑2197351
                                                                      |                                                                     
aGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgcg                                                                       <  1:1359472/71‑1 (MQ=255)
aGTTGTCAAAGGAAGAACAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgcg                                                                       <  1:118179/71‑1 (MQ=255)
   tGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATc                                                                    >  1:2455446/1‑71 (MQ=255)
   tGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATc                                                                    >  1:2916368/1‑71 (MQ=255)
    gTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                                                   >  1:1442458/1‑71 (MQ=255)
     tCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgcg                                                                       <  1:2534021/66‑1 (MQ=255)
     tCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                                                   <  1:2772911/70‑1 (MQ=255)
     tCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                                                   <  1:4545552/70‑1 (MQ=255)
     tCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                                                   <  1:4186494/70‑1 (MQ=255)
     tCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                                                   <  1:1356740/70‑1 (MQ=255)
      cAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                                                 >  1:4351188/1‑71 (MQ=255)
      cAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                                                 >  1:1175341/1‑71 (MQ=255)
       aaaGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGt                                                                <  1:1309202/71‑1 (MQ=255)
         aGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGa                                                              <  1:3583421/71‑1 (MQ=255)
         aGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGtt                                                               >  1:3701930/1‑69 (MQ=255)
         aGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGtt                                                               >  1:1723931/1‑69 (MQ=255)
             aGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          >  1:1633437/1‑71 (MQ=255)
             aGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          >  1:2677058/1‑71 (MQ=255)
              gCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCg                                                         <  1:1281757/71‑1 (MQ=255)
               cccAGCTGCTCAGGCTGTCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCTCGATCGGGTGAACTc                                                          <  1:2621896/69‑1 (MQ=255)
               cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:4555545/1‑68 (MQ=255)
               cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:3139710/1‑68 (MQ=255)
               cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:949032/1‑68 (MQ=255)
               cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:873860/1‑68 (MQ=255)
               cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:3575007/1‑68 (MQ=255)
               cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:836400/69‑1 (MQ=255)
               cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:130084/69‑1 (MQ=255)
                ccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                                        <  1:731395/70‑1 (MQ=255)
                   gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                                                     >  1:66832/1‑70 (MQ=255)
                   gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:1007753/1‑71 (MQ=255)
                   gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:399904/1‑71 (MQ=255)
                   gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:2199082/1‑71 (MQ=255)
                   gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:4708494/1‑71 (MQ=255)
                           ggctggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                                        >  1:522266/1‑59 (MQ=255)
                           ggctggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                                        >  1:4097010/1‑59 (MQ=255)
                           ggctggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                                        >  1:3263125/1‑59 (MQ=255)
                              tggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGt                                          >  1:1827609/1‑70 (MQ=255)
                                gctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         >  1:4565017/1‑69 (MQ=255)
                                gctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         >  1:3025641/1‑69 (MQ=255)
                                     cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:2053287/41‑1 (MQ=255)
                                     cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:2622244/41‑1 (MQ=255)
                                     cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:1513480/41‑1 (MQ=255)
                                     cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:4293419/41‑1 (MQ=255)
                                     cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:3074498/41‑1 (MQ=255)
                                      aaaaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATCGCTGAAATGCGTa                                         <  1:3922349/63‑1 (MQ=255)
                                      aaaaaGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:3884757/40‑1 (MQ=255)
                                        aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1844940/1‑63 (MQ=255)
                                        aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1621444/1‑63 (MQ=255)
                                        aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc                               <  1:4614624/71‑1 (MQ=255)
                                        aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc                               >  1:4688562/1‑71 (MQ=255)
                                         aaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggcg                              >  1:4180075/1‑71 (MQ=255)
                                           gACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggcgg                             >  1:3594395/1‑70 (MQ=255)
                                             cGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa                          <  1:2500264/71‑1 (MQ=255)
                                             cGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa                          <  1:2878843/71‑1 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:1826366/38‑1 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:1073768/38‑1 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:4149739/38‑1 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:953301/1‑57 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3025574/1‑57 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3821256/1‑57 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1385221/1‑57 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3364358/1‑57 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3020210/1‑57 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaa                         <  1:1237370/71‑1 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaa                         <  1:489396/71‑1 (MQ=255)
                                              gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa                          >  1:1233728/1‑70 (MQ=255)
                                                 gtgATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTg                                 <  1:3792604/60‑1 (MQ=255)
                                                 gtgATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa                          <  1:2376989/67‑1 (MQ=255)
                                                 gtgATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAAcc                      >  1:3372401/1‑71 (MQ=255)
                                                    aTCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGc                   >  1:1944778/1‑71 (MQ=255)
                                                        ttGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa                                        <  1:365199/46‑1 (MQ=255)
                                                        ttGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa                                        <  1:1776611/46‑1 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:836468/1‑46 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:1008171/1‑46 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:1899771/1‑46 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:1220470/1‑46 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:2724760/1‑46 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatcat             >  1:2878785/1‑71 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatca              >  1:2315871/1‑70 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatca              <  1:1654515/70‑1 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatca              <  1:1976073/70‑1 (MQ=255)
                                                          ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatca              <  1:3753970/70‑1 (MQ=255)
                                                           gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         >  1:2596919/1‑42 (MQ=255)
                                                           gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         <  1:2170593/42‑1 (MQ=255)
                                                           gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt                                  >  1:1426099/1‑49 (MQ=255)
                                                           gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt                                  >  1:2821123/1‑49 (MQ=255)
                                                             gTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:2852398/1‑42 (MQ=255)
                                                             gTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1757832/1‑42 (MQ=255)
                                                             gTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGc          >  1:528209/1‑71 (MQ=255)
                                                              tCAACCGCGATCGGGTGAACTCTATTGCTGAAATGCGt                                          <  1:4088575/38‑1 (MQ=255)
                                                              tCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGcc         <  1:224296/71‑1 (MQ=255)
                                                                aaCCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCt       >  1:3697363/1‑71 (MQ=255)
                                                                  ccGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTg      >  1:1988591/1‑70 (MQ=255)
                                                                   cgcgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaaa                        >  1:212557/1‑51 (MQ=255)
                                                                   cgcgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCg           >  1:1484014/1‑64 (MQ=255)
                                                                    gcgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCa    <  1:428344/70‑1 (MQ=255)
                                                                     cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  <  1:68304/71‑1 (MQ=255)
                                                                     cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  >  1:2806958/1‑71 (MQ=255)
                                                                     cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  <  1:1417274/71‑1 (MQ=255)
                                                                     cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  <  1:1024477/71‑1 (MQ=255)
                                                                      |                                                                     
AGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAA  >  Exported/2197212‑2197351

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