Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,489,541 G→A S171S (AGC→AGT glpG ← predicted intramembrane serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,489,5410GA98.4% 728.6 / ‑1.0 190S171S (AGC→AGTglpGpredicted intramembrane serine protease
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (99/88);  total (101/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 200 total aligned reads displayed.

ACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGC  >  REL606/3489468‑3489613
                                                                         |                                                                        
aCGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACt                                                                        <  1:4210108/76‑1 (MQ=255)
 cGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTa                                                                       <  2:2850655/76‑1 (MQ=255)
 cGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTa                                                                       <  2:4923180/76‑1 (MQ=255)
 cGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTa                                                                       <  1:2034823/76‑1 (MQ=255)
  ggCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTAc                                                                      >  2:6609327/1‑76 (MQ=255)
  ggCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTAc                                                                      >  2:6095494/1‑76 (MQ=255)
    cccGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCg                                                                    <  2:3088182/76‑1 (MQ=255)
     ccGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGa                                                                   <  2:1523456/76‑1 (MQ=255)
      cGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAg                                                                  >  1:1746884/1‑76 (MQ=255)
       gCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAgg                                                                 <  2:1297161/76‑1 (MQ=255)
          gAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGt                                                              >  2:2905843/1‑76 (MQ=255)
           aaTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGtt                                                             <  2:135814/76‑1 (MQ=255)
           aaTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGtt                                                             <  2:1954780/76‑1 (MQ=255)
            aTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGttt                                                            >  1:5004038/1‑76 (MQ=255)
             tttttGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGtttt                                                           <  2:3922290/76‑1 (MQ=255)
               tttGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTc                                                         >  2:4432538/1‑76 (MQ=255)
                ttGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTcc                                                        <  1:2275342/76‑1 (MQ=255)
                 tgctgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCa                                                       <  1:2442733/76‑1 (MQ=255)
                 tgctgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCa                                                       >  2:2216487/1‑76 (MQ=255)
                 tgctgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCa                                                       <  2:6512510/76‑1 (MQ=255)
                  gctgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCAc                                                      >  1:4837193/1‑76 (MQ=255)
                   ctgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCAcc                                                     >  2:2995716/1‑76 (MQ=255)
                   ctgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCAcc                                                     >  2:2686243/1‑76 (MQ=255)
                    tgtacaTAGACGCTTAACAGGGGGATAATGACAGTAATAAGAATTAGCTTAAAGCTACCGAGGCGTTTTTCAACcg                                                    >  2:1603411/4‑76 (MQ=255)
                    tgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACcg                                                    >  1:4061532/1‑76 (MQ=255)
                    tgcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACcg                                                    <  2:3429237/76‑1 (MQ=255)
                     gcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACcgc                                                   >  1:4768954/1‑76 (MQ=255)
                     gcACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACcgc                                                   >  2:3394755/1‑76 (MQ=255)
                       acaTAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCgcgc                                                 >  2:1912735/1‑76 (MQ=255)
                       acaTAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCgcgc                                                 >  2:562341/1‑76 (MQ=255)
                         aTAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGcgccg                                               >  1:4949101/1‑76 (MQ=255)
                         aTAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGcgccg                                               >  2:6180853/1‑76 (MQ=255)
                          tAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACAGcgccgc                                              >  2:4908798/1‑76 (MQ=255)
                           aGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCgccgcc                                             >  2:1196615/1‑76 (MQ=255)
                           aGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCgccgcc                                             <  2:5353865/76‑1 (MQ=255)
                            gCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCCCCGCGccgccg                                            <  2:694458/76‑1 (MQ=255)
                             ccGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCga                                           <  1:1157943/76‑1 (MQ=255)
                               gCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCgaga                                         <  1:1014872/76‑1 (MQ=255)
                                 ttAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATa                                       <  1:984841/76‑1 (MQ=255)
                                 ttAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATa                                       >  1:6452828/1‑76 (MQ=255)
                                   aaCAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATacc                                     >  1:3950168/1‑76 (MQ=255)
                                   aaCAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACAGAGGCGTTTTTCCACCGCGCCGCCGAGATacc                                     >  2:5414810/1‑76 (MQ=255)
                                    agcgggCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATacca                                    <  2:2262360/73‑1 (MQ=255)
                                    aCAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATacca                                    <  1:2828987/76‑1 (MQ=255)
                                     cAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATaccac                                   >  2:244383/1‑76 (MQ=255)
                                      aGGGCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATaccacc                                  >  2:4505082/1‑76 (MQ=255)
                                       gggCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATaccacca                                 <  1:4588595/76‑1 (MQ=255)
                                       gggCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATaccacca                                 <  2:3946883/76‑1 (MQ=255)
                                        ggCGCTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATaccaccac                                <  1:3689841/76‑1 (MQ=255)
                                          cgcTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATAccaccacca                              >  2:2149381/1‑76 (MQ=255)
                                          cgcTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATAccaccacca                              <  1:1727874/76‑1 (MQ=255)
                                          cgcTAATGAGAGAAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATAccaccacca                              >  1:1495487/1‑76 (MQ=255)
                                           gcTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCag                             <  1:4727638/76‑1 (MQ=255)
                                            cTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCaga                            >  2:6419540/1‑76 (MQ=255)
                                            cTAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCAACGCGCCGCCGACATACCACCGCCagc                            >  2:541777/1‑75 (MQ=255)
                                             tAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCagag                           >  1:839333/1‑76 (MQ=255)
                                             tAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCagag                           >  1:1665458/1‑76 (MQ=255)
                                             tAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCagag                           <  1:4050255/76‑1 (MQ=255)
                                             tAATGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCagag                           <  1:5499305/76‑1 (MQ=255)
                                              aaTGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGc                          <  1:6583637/76‑1 (MQ=255)
                                              aaTGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGc                          >  1:1446108/1‑76 (MQ=255)
                                              aaTGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGc                          <  2:1441138/76‑1 (MQ=255)
                                              aaTGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGc                          >  1:4273625/1‑76 (MQ=255)
                                              aaTGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGc                          >  2:2732433/1‑76 (MQ=255)
                                               aTGAGATTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCa                         <  1:4133447/76‑1 (MQ=255)
                                               aTGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGCTCCCACCACCAGAGCa                         >  1:5177750/1‑76 (MQ=255)
                                                tGAGAGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAg                        >  2:4221993/1‑76 (MQ=255)
                                                 gagaGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAgg                       <  2:6558443/76‑1 (MQ=255)
                                                 gagaGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAgg                       >  1:182190/1‑76 (MQ=255)
                                                 gagaGTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAgg                       <  1:3839854/76‑1 (MQ=255)
                                                  agagTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGt                      >  2:5893977/1‑76 (MQ=255)
                                                   gagTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGtt                     <  2:3659879/76‑1 (MQ=255)
                                                    agTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCCCCACCAGAGCAGGTTa                    <  2:4768954/76‑1 (MQ=255)
                                                     gTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTaa                   <  1:4397640/76‑1 (MQ=255)
                                                     gTAATAACAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTaa                   <  2:5709567/76‑1 (MQ=255)
                                                     ataataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGTCGAGATACCACCACCAGAGCAGGTTaa                   <  1:3440737/75‑1 (MQ=255)
                                                      taataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTaaa                  >  1:1334702/1‑76 (MQ=255)
                                                      taataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTaaa                  >  1:5607064/1‑76 (MQ=255)
                                                       aataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAag                 >  2:772831/1‑76 (MQ=255)
                                                       aataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAag                 >  2:2322256/1‑76 (MQ=255)
                                                       aataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAag                 >  1:2750192/1‑76 (MQ=255)
                                                       aataaCAATTAGCTTACCACTACAGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAag                 <  1:5414810/76‑1 (MQ=255)
                                                        ataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAaga                <  1:2226359/76‑1 (MQ=255)
                                                        ataaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAaga                >  1:3108083/1‑76 (MQ=255)
                                                         taaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAagag               >  1:3329017/1‑76 (MQ=255)
                                                          aaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAgg              >  1:3574890/1‑76 (MQ=255)
                                                          aaCAATTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAgg              <  1:6100585/76‑1 (MQ=255)
                                                          aaCAATGAGCTTACCCCCACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAgg              <  1:5883238/76‑1 (MQ=255)
                                                             aaTTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGata           >  1:5361771/1‑76 (MQ=255)
                                                             aaTTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGata           <  1:4505082/76‑1 (MQ=255)
                                                             aaTTAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGata           <  1:294373/76‑1 (MQ=255)
                                                               ttAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATg         <  1:3723760/76‑1 (MQ=255)
                                                                tAGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGc        >  2:4354636/1‑76 (MQ=255)
                                                                 aGCTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCa       >  2:5033309/1‑76 (MQ=255)
                                                                   cTTACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATc     >  1:192396/1‑76 (MQ=255)
                                                                    ttACCACTAGCGAGGCGTTTTTCCACCTCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCa    <  2:1735152/76‑1 (MQ=255)
                                                                    ttACCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCa    <  1:6534467/76‑1 (MQ=255)
                                                                      aCCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGc  >  2:331482/1‑76 (MQ=255)
                                                                      aCCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGc  >  1:4426498/1‑76 (MQ=255)
                                                                      aCCACTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGc  >  1:4981149/1‑76 (MQ=255)
                                                                         |                                                                        
ACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTAATGAGAGTAATAACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCGCCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGC  >  REL606/3489468‑3489613

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