Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,698,011 G→A S3L (TCG→TTG)  slyA ← transcriptional regulator SlyA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,698,0110GA98.8% 1045.1 / 1.1 247S3L (TCG→TTG) slyAtranscriptional regulator SlyA
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (89/155);  total (91/156)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 248 total aligned reads displayed.

CAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTG  >  REL606/1697938‑1698083
                                                                         |                                                                        
cAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAAt                                                                        <  2:1993347/76‑1 (MQ=255)
 aGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAAtt                                                                       >  2:4425680/1‑76 (MQ=255)
 aGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAAtt                                                                       <  1:4248067/76‑1 (MQ=255)
  gCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTc                                                                      <  1:6052980/76‑1 (MQ=255)
  gCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTc                                                                      <  2:431459/76‑1 (MQ=255)
    ggTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCa                                                                    <  1:2227496/76‑1 (MQ=255)
    ggTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCa                                                                    <  1:2880179/76‑1 (MQ=255)
       ttCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAAtt                                                                 >  1:1640909/1‑76 (MQ=255)
       ttCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAAtt                                                                 <  1:2163321/76‑1 (MQ=255)
       ttCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAAtt                                                                 >  1:4699724/1‑76 (MQ=255)
        tCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAAttt                                                                <  2:333923/76‑1 (MQ=255)
        tCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAAttt                                                                >  1:2074078/1‑76 (MQ=255)
           ggCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCAt                                                             <  1:5262896/76‑1 (MQ=255)
           ggCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCAt                                                             <  2:3651974/76‑1 (MQ=255)
            gCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCAtc                                                            <  1:6092357/76‑1 (MQ=255)
            gCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCAtc                                                            <  2:4243731/76‑1 (MQ=255)
             cGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCAtct                                                           >  2:3705573/1‑76 (MQ=255)
             cGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCAtct                                                           >  2:68675/1‑76 (MQ=255)
               aTGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTcc                                                         >  2:3838898/1‑76 (MQ=255)
                 ggTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCtt                                                       >  1:5665091/1‑76 (MQ=255)
                  gTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTta                                                      <  2:7598444/76‑1 (MQ=255)
                   tCTATCAGAGCACGCCATATGCGCACCAACCGTTCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTtat                                                     <  1:3175333/76‑1 (MQ=255)
                   tCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTtat                                                     <  2:1033044/76‑1 (MQ=255)
                   tCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTtat                                                     >  2:2373228/1‑76 (MQ=255)
                   tCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTtat                                                     <  1:5115841/76‑1 (MQ=255)
                    cTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTtata                                                    <  1:2953261/76‑1 (MQ=255)
                    cTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTtata                                                    >  1:1405440/1‑76 (MQ=255)
                     tATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATaa                                                   <  2:1483133/76‑1 (MQ=255)
                     tATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATaa                                                   <  2:7519301/76‑1 (MQ=255)
                     tATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATaa                                                   >  2:1942328/1‑76 (MQ=255)
                     tATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATaa                                                   <  1:4580744/76‑1 (MQ=255)
                       tCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAAtt                                                 <  2:5210412/76‑1 (MQ=255)
                       tCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAAtt                                                 >  2:6331340/1‑76 (MQ=255)
                       tCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAAtt                                                 <  1:120226/76‑1 (MQ=255)
                        cAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTa                                                >  1:7520480/1‑76 (MQ=255)
                        cAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTa                                                >  2:2585132/1‑76 (MQ=255)
                         agagCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAg                                               <  1:3940932/76‑1 (MQ=255)
                           agCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCt                                             <  1:645446/76‑1 (MQ=255)
                                gCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTa                                        <  2:1372756/76‑1 (MQ=255)
                                gCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTa                                        <  1:3465691/76‑1 (MQ=255)
                                gCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTa                                        <  1:7362585/76‑1 (MQ=255)
                                 ccATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTaa                                       <  1:3092919/76‑1 (MQ=255)
                                  cATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAg                                      >  2:2826064/1‑76 (MQ=255)
                                    tatGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGat                                    <  1:7092727/76‑1 (MQ=255)
                                    tatGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGat                                    <  2:2245299/76‑1 (MQ=255)
                                    tatGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGat                                    <  1:3788368/76‑1 (MQ=255)
                                     atGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGata                                   >  1:1261293/1‑76 (MQ=255)
                                     atGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGata                                   <  1:4647202/76‑1 (MQ=255)
                                      tGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGatat                                  <  2:487312/76‑1 (MQ=255)
                                       gcgcACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGAtatt                                 <  2:821217/76‑1 (MQ=255)
                                       gcgcACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGAtatt                                 <  2:5132201/76‑1 (MQ=255)
                                       gcgcACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGAtatt                                 <  2:5396148/76‑1 (MQ=255)
                                        cgcACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGAtatta                                >  2:6757233/1‑76 (MQ=255)
                                        cgcACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGAtatta                                <  1:80889/76‑1 (MQ=255)
                                         gcACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGAtattat                               <  1:1687283/76‑1 (MQ=255)
                                         gcACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGAtattat                               <  2:6274985/76‑1 (MQ=255)
                                          cACCAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATg                              >  2:3191647/1‑76 (MQ=255)
                                            ccAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCg                            >  1:6983056/1‑76 (MQ=255)
                                             cAACCGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCgg                           >  2:4151580/1‑76 (MQ=255)
                                              aaCCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGc                          >  2:6482737/1‑76 (MQ=255)
                                                ccGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCtt                        <  1:4877554/76‑1 (MQ=255)
                                                 cGTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCttt                       <  1:6583924/76‑1 (MQ=255)
                                                  gTGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCtttt                      <  2:911017/76‑1 (MQ=255)
                                                   tGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTa                     <  2:5529711/76‑1 (MQ=255)
                                                   tGCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTa                     <  2:2810074/76‑1 (MQ=255)
                                                    gCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAg                    <  1:2479727/76‑1 (MQ=255)
                                                    gCCAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAg                    <  1:7382353/76‑1 (MQ=255)
                                                      cAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGaa                  <  2:1846694/76‑1 (MQ=255)
                                                      cAGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGaa                  <  2:608152/76‑1 (MQ=255)
                                                       aGATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAAt                 <  1:2101164/76‑1 (MQ=255)
                                                        gATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATa                <  1:5018581/76‑1 (MQ=255)
                                                        gATCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATa                <  2:4697207/76‑1 (MQ=255)
                                                         aTCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAg               <  1:432378/76‑1 (MQ=255)
                                                         aTCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAg               <  2:3979230/76‑1 (MQ=255)
                                                          tCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAgt              <  1:4108860/76‑1 (MQ=255)
                                                          tCAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAgt              >  2:2188508/1‑76 (MQ=255)
                                                           cAGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAgtg             >  2:2148755/1‑76 (MQ=255)
                                                            aGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAgtgt            <  1:350007/76‑1 (MQ=255)
                                                            aGAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAgtgt            >  2:7085306/1‑76 (MQ=255)
                                                             gAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGtgtg           <  2:4567321/76‑1 (MQ=255)
                                                             gAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGtgtg           <  2:5759020/76‑1 (MQ=255)
                                                             gAACCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGtgtg           >  2:2736889/1‑76 (MQ=255)
                                                              aaCCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGtgtgg          >  2:6031117/1‑75 (MQ=255)
                                                              aaCCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgc          >  1:4162649/1‑76 (MQ=255)
                                                               aCCTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgca         <  1:2608621/76‑1 (MQ=255)
                                                                ccTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag        <  2:4472914/76‑1 (MQ=255)
                                                                ccTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag        >  1:327278/1‑76 (MQ=255)
                                                                ccTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag        >  2:1738487/1‑76 (MQ=255)
                                                                 cTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagc       >  1:2807106/1‑76 (MQ=255)
                                                                 cTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagc       >  2:5696557/1‑76 (MQ=255)
                                                                 cTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagc       <  1:549368/76‑1 (MQ=255)
                                                                 cTAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagc       >  2:1153555/1‑76 (MQ=255)
                                                                  tAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagca      <  1:1174256/76‑1 (MQ=255)
                                                                  tAGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagca      <  1:419300/76‑1 (MQ=255)
                                                                   tgTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCaa     >  2:6599894/2‑76 (MQ=255)
                                                                   aGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCaa     >  2:4858620/1‑76 (MQ=255)
                                                                   aGTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCaa     <  2:7170261/76‑1 (MQ=255)
                                                                    gTGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAAt    <  2:125031/76‑1 (MQ=255)
                                                                     tGGCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAAtt   <  1:1101083/76‑1 (MQ=255)
                                                                      ggCAATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTg  >  1:4900638/1‑76 (MQ=255)
                                                                         |                                                                        
CAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTG  >  REL606/1697938‑1698083

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