Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,489,365 Δ1 bp coding (689/831 nt) glpG ← predicted intramembrane serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,489,3640A.99.0% 982.5 / 2.5 194coding (690/831 nt)glpGpredicted intramembrane serine protease
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (85/107);  total (86/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 207 total aligned reads displayed.

GATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAG  >  REL606/3489289‑3489435
                                                                           |                                                                       
taTGGGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATa                                                                         <  1:2047823/75‑1 (MQ=255)
gATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATa                                                                         <  2:3611801/76‑1 (MQ=255)
 aTGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATat                                                                        <  1:5429067/76‑2 (MQ=255)
 aTGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATat                                                                        >  1:3097833/1‑75 (MQ=255)
  tgtgtgCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATatt                                                                       >  1:1631605/1‑74 (MQ=255)
   gtgtcCCCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATatta                                                                      <  2:5706688/76‑4 (MQ=255)
     gtgCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAAc                                                                   >  1:4233825/1‑76 (MQ=255)
     gtgCTCCGTTCCCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAAc                                                                   <  1:6241301/76‑1 (MQ=255)
     gtccTCCTTCCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAAc                                                                   <  2:4042568/73‑1 (MQ=255)
      tgCTCCGTCCCCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAAcc                                                                  <  1:6325585/76‑1 (MQ=255)
      tgCCCCGTTCCCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAAcc                                                                  <  1:4880240/76‑1 (MQ=255)
       gCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAAccc                                                                 >  1:1437882/1‑76 (MQ=255)
        cTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCa                                                                <  2:3027453/76‑1 (MQ=255)
        cTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCa                                                                <  2:1504103/76‑1 (MQ=255)
        cTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCa                                                                <  2:7425723/76‑1 (MQ=255)
        cTCCGTCCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCa                                                                <  2:5327495/76‑1 (MQ=255)
         tCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCAc                                                               >  1:2470653/1‑76 (MQ=255)
          ccGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACg                                                              >  2:6385311/1‑76 (MQ=255)
          ccGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACg                                                              >  1:5606285/1‑76 (MQ=255)
             ttCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTg                                                           >  1:3119386/1‑76 (MQ=255)
             ttCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTg                                                           <  2:236679/76‑1 (MQ=255)
             ttCCCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTg                                                           <  2:4926203/76‑1 (MQ=255)
               cGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCa                                                         <  2:651247/76‑1 (MQ=255)
                  cATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGt                                                      <  1:4344135/76‑1 (MQ=255)
                   aTCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTa                                                     <  2:902778/76‑1 (MQ=255)
                     cGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTaaa                                                   >  1:6037323/1‑76 (MQ=255)
                      gACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAAt                                                  <  2:5415512/76‑1 (MQ=255)
                       aCATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATg                                                 >  2:6720027/1‑76 (MQ=255)
                       aCATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATg                                                 >  1:1996516/1‑76 (MQ=255)
                         aTCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGcc                                               >  2:5457772/1‑76 (MQ=255)
                          tCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCa                                              >  2:2700505/1‑76 (MQ=255)
                          tCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCa                                              >  2:7125318/1‑76 (MQ=255)
                            ccAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACt                                            <  2:5245831/76‑1 (MQ=255)
                              aaacaaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACttt                                          >  2:4349058/1‑76 (MQ=255)
                              aaacaaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACttt                                          <  1:6996883/76‑1 (MQ=255)
                                acaaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGc                                        <  1:2260696/76‑1 (MQ=255)
                                acaaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGc                                        <  1:7313488/76‑1 (MQ=255)
                                 caaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCg                                       >  2:7352227/1‑76 (MQ=255)
                                 caaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCg                                       >  2:5924854/1‑76 (MQ=255)
                                  aaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTACAGGTAAATGCCACTTTGCgg                                      >  2:209912/1‑76 (MQ=255)
                                   aaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGa                                     <  2:331001/76‑1 (MQ=255)
                                    aTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGAt                                    <  2:7211038/76‑1 (MQ=255)
                                     tCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATc                                   <  2:7254395/76‑1 (MQ=255)
                                      cAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcg                                  >  1:6857587/1‑76 (MQ=255)
                                       aaaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgc                                 >  2:2514680/1‑76 (MQ=255)
                                       aaaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgc                                 <  2:5081854/76‑1 (MQ=255)
                                       aaaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgc                                 >  2:3907835/1‑76 (MQ=255)
                                       aaaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgc                                 >  1:3218089/1‑76 (MQ=255)
                                       aaaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgc                                 >  2:4173650/1‑76 (MQ=255)
                                       aaaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgc                                 >  2:5940698/1‑76 (MQ=255)
                                       aaaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgc                                 >  1:4450798/1‑76 (MQ=255)
                                        aaCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATcgcg                                >  1:1167375/1‑76 (MQ=255)
                                         aCCATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGt                               <  2:3387404/76‑1 (MQ=255)
                                           cATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTc                             <  2:1832765/76‑1 (MQ=255)
                                           cATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTc                             <  2:1436655/76‑1 (MQ=255)
                                           cATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTc                             >  2:4556248/1‑76 (MQ=255)
                                           cATCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTc                             <  2:3515267/76‑1 (MQ=255)
                                           cATCAGGCGACAATCCAGATCAGCGCAAGGATATTAACCCACGTTGCAGGTAAATTCCACTATGCGTATCAGGTTc                             >  2:31900/1‑76 (MQ=255)
                                            aTCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCg                            >  2:6954562/1‑76 (MQ=255)
                                            aTCCGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCg                            >  1:5811513/1‑76 (MQ=255)
                                               cGGCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCa                         >  2:262260/1‑76 (MQ=255)
                                                ggCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCAc                        >  1:2736985/1‑76 (MQ=255)
                                                 gCGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACg                       <  1:2032225/76‑1 (MQ=255)
                                                  cGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGt                      >  2:3704179/1‑76 (MQ=255)
                                                  cGACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGAATCGCGTTCGCCACgc                      >  2:7034767/1‑75 (MQ=255)
                                                   gACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTa                     <  1:2109813/76‑1 (MQ=255)
                                                   gACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTa                     <  1:4130926/76‑1 (MQ=255)
                                                   gACAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTa                     <  1:4040707/76‑1 (MQ=255)
                                                    aCAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAg                    >  1:6659731/1‑76 (MQ=255)
                                                    aCAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAg                    >  1:3817389/1‑76 (MQ=255)
                                                     cAATCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGc                   >  2:415601/1‑76 (MQ=255)
                                                      aaTCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGcc                  <  2:429744/76‑1 (MQ=255)
                                                       aTCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCa                 <  1:314610/76‑1 (MQ=255)
                                                       aTCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCa                 >  1:5374510/1‑76 (MQ=255)
                                                       aTCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCa                 <  2:2568151/76‑1 (MQ=255)
                                                        tCCAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAg                >  1:5913080/1‑76 (MQ=255)
                                                          cAGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGAc              >  1:430698/1‑76 (MQ=255)
                                                           aGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACg             <  2:7398867/76‑1 (MQ=255)
                                                           aGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACg             >  2:5521130/1‑76 (MQ=255)
                                                           aGATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACg             >  2:5782557/1‑76 (MQ=255)
                                                            gATCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGt            >  2:6060393/1‑76 (MQ=255)
                                                             aTCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTa           >  1:7096846/1‑76 (MQ=255)
                                                             aTCAGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGCAGCCACACgca           >  1:50884/1‑74 (MQ=255)
                                                              acAGCGCAAAGATATTAACCCCGTTTGCAGGTAAATGCCACTTTGCGGATCGCGTGCGCCCCGTAGCCAGACGTAg          <  2:3122416/75‑1 (MQ=255)
                                                               cAGCGCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGCAg          >  2:3937391/1‑76 (MQ=255)
                                                                aGCGCAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGcc        >  2:6814396/1‑76 (MQ=255)
                                                                 gcgcAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGccc       >  2:6144396/1‑76 (MQ=255)
                                                                 gcgcAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGccc       >  1:7395636/1‑76 (MQ=255)
                                                                 gcgcAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTGGCGGATCGCGTTCGCCACGTAGCCAGACGTAGccc       <  1:6456201/76‑1 (MQ=255)
                                                                  cgcAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCa      <  2:6638729/76‑1 (MQ=255)
                                                                  cgcAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCa      >  2:1118365/1‑76 (MQ=255)
                                                                   gcAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCAt     <  1:4813933/76‑1 (MQ=255)
                                                                   gcAAAGATATTAACCCACGTTGCAGGTAAAGGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCAt     <  1:1923683/76‑1 (MQ=255)
                                                                    cAAAGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATc    >  2:158774/1‑76 (MQ=255)
                                                                     aaaGATATTACCCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCa   <  2:4730544/76‑1 (MQ=255)
                                                                     aaaGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCa   <  2:3622077/76‑1 (MQ=255)
                                                                     aaaGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCa   <  1:4994925/76‑1 (MQ=255)
                                                                     aaaGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCa   <  2:5233436/76‑1 (MQ=255)
                                                                      aaGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAg  >  1:674372/1‑76 (MQ=255)
                                                                      aaGATATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAg  <  2:2209409/76‑1 (MQ=255)
                                                                           |                                                                       
GATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAG  >  REL606/3489289‑3489435

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