Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,197,287 G→A 7.5% R421R (CGG→CGA degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,197,2870GA 7.5% 385.7 / 15.4 120R421R (CGG→CGAdegQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (74/37);  new base A (3/6);  total (77/43)
Fisher's exact test for biased strand distribution p-value = 6.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 124 total aligned reads displayed.

CAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAA  >  Exported/2197218‑2197351
                                                                     |                                                                
cAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCggg                                                                 >  1:4351188/1‑71 (MQ=255)
 aaaGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGt                                                                <  1:2544436/71‑1 (MQ=255)
 aaaGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGt                                                                <  1:1309202/71‑1 (MQ=255)
   aGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTg                                                               >  1:2262163/1‑70 (MQ=255)
   aGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGtt                                                               >  1:1723931/1‑69 (MQ=255)
     gAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAAc                                                            >  1:2417863/1‑71 (MQ=255)
       aGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGaa                                                             <  1:1214431/68‑1 (MQ=255)
       aGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          >  1:1633437/1‑71 (MQ=255)
       aGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          >  1:4008119/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:3139710/1‑68 (MQ=255)
         cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:873860/1‑68 (MQ=255)
         cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:3575007/1‑68 (MQ=255)
         cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:949032/1‑68 (MQ=255)
         cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACt                                                           >  1:4555545/1‑68 (MQ=255)
         cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:130084/69‑1 (MQ=255)
         cccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:836400/69‑1 (MQ=255)
          ccAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                                        <  1:731395/70‑1 (MQ=255)
             gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:2199082/1‑71 (MQ=255)
             gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:399904/1‑71 (MQ=255)
             gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:4708494/1‑71 (MQ=255)
             gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:285772/1‑71 (MQ=255)
             gctgctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGc                                                    >  1:1007753/1‑71 (MQ=255)
                gctCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTg                                                  >  1:106464/1‑70 (MQ=255)
                     ggctggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                                        >  1:522266/1‑59 (MQ=255)
                     ggctggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGa                                                        >  1:4097010/1‑59 (MQ=255)
                        tggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGt                                          >  1:1827609/1‑70 (MQ=255)
                        tggctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGt                                          >  1:146236/1‑70 (MQ=255)
                          gctTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         >  1:3025641/1‑69 (MQ=255)
                               cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:2622244/41‑1 (MQ=255)
                               cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:3074498/41‑1 (MQ=255)
                               cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:2053287/41‑1 (MQ=255)
                               cAAAAAGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:4293419/41‑1 (MQ=255)
                                aaaaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTg                                 <  1:4758394/71‑1 (MQ=255)
                                aaaaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATCGCTGAAATGCGTa                                         <  1:3922349/63‑1 (MQ=255)
                                aaaaaGACGATGTGATCATTGGCGTCAACCGCAATCGAGt                                                                <  1:3884757/40‑1 (MQ=255)
                                  aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1621444/1‑63 (MQ=255)
                                  aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3738393/1‑63 (MQ=255)
                                  aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1844940/1‑63 (MQ=255)
                                  aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc                               <  1:4614624/71‑1 (MQ=255)
                                  aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc                               >  1:4688562/1‑71 (MQ=255)
                                  aaaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc                               >  1:1528502/1‑71 (MQ=255)
                                   aaGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggcg                              >  1:4180075/1‑71 (MQ=255)
                                     gACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGAtt                                                      >  1:3643351/1‑45 (MQ=255)
                                     gACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggcgg                             >  1:3594395/1‑70 (MQ=255)
                                       cGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa                          <  1:2500264/71‑1 (MQ=255)
                                       cGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa                          <  1:2878843/71‑1 (MQ=255)
                                       cGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCa                           >  1:3474255/1‑70 (MQ=255)
                                       cGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCa                           >  1:82442/1‑70 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:4149739/38‑1 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTc                                                          <  1:1073768/38‑1 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3821256/1‑57 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3364358/1‑57 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:953301/1‑57 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3020210/1‑57 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1385221/1‑57 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:1880845/1‑57 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaaa                                       >  1:3025574/1‑57 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaaa                         <  1:1237370/71‑1 (MQ=255)
                                        gATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa                          >  1:1983881/1‑70 (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)
                                              aTCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         <  1:4159054/49‑1 (MQ=255)
                                              aTCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGc                   >  1:1944778/1‑71 (MQ=255)
                                                  ttGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa                                        <  1:365199/46‑1 (MQ=255)
                                                  ttGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa                                        <  1:1790934/46‑1 (MQ=255)
                                                  ttGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa                                        <  1:1776611/46‑1 (MQ=255)
                                                    ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:1899771/1‑46 (MQ=255)
                                                    ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:1220470/1‑46 (MQ=255)
                                                    ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:1008171/1‑46 (MQ=255)
                                                    ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAg                                      >  1:2724760/1‑46 (MQ=255)
                                                    ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatcat             >  1:549494/1‑71 (MQ=255)
                                                    ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatca              <  1:3753970/70‑1 (MQ=255)
                                                    ggCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCcatca              <  1:1654515/70‑1 (MQ=255)
                                                     gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         <  1:2170593/42‑1 (MQ=255)
                                                     gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTa                                         <  1:4247965/42‑1 (MQ=255)
                                                     gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt                                  >  1:440553/1‑49 (MQ=255)
                                                     gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt                                  >  1:2821123/1‑49 (MQ=255)
                                                     gCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCt                                  >  1:1426099/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)
                                                        tCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGcc         <  1:224296/71‑1 (MQ=255)
                                                          aaCCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCt       >  1:3191171/1‑71 (MQ=255)
                                                          aaCCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCt       >  1:936946/1‑71 (MQ=255)
                                                          aaCCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCt       >  1:2931701/1‑71 (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)
                                                               cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  <  1:1417274/71‑1 (MQ=255)
                                                               cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  <  1:68304/71‑1 (MQ=255)
                                                               cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  >  1:2806958/1‑71 (MQ=255)
                                                               cgATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCaaa  <  1:1024477/71‑1 (MQ=255)
                                                                    gggTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCa    <  1:540601/64‑1 (MQ=255)
                                                                    gggTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCa    <  1:3199587/64‑1 (MQ=255)
                                                                     |                                                                
CAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAACCGGCCATCATCGCCCTGCAAA  >  Exported/2197218‑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: