Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,495,556 A→C intergenic (‑96/+32) yzgL ← / ← glgP hypothetical protein; b3427_1/glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,495,5560AC99.1% 914.7 / ‑2.0 213intergenic (‑96/+32)yzgL/glgPhypothetical protein; b3427_1/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (114/97);  total (114/99)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 219 total aligned reads displayed.

TGGCGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCGGCGT  >  REL606/3495483‑3495628
                                                                         |                                                                        
tGGCGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGccc                                                                        <  1:1781757/76‑1 (MQ=255)
  gCGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCt                                                                      >  2:3476655/1‑76 (MQ=255)
   cGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCtt                                                                     >  1:2477321/1‑76 (MQ=255)
   cGAAATGATTCTCCCCGAGTCCCCGTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCtt                                                                     <  2:7471276/76‑1 (MQ=255)
    gAAATTATCCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCttt                                                                    <  1:3022814/76‑1 (MQ=255)
    gAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCttt                                                                    >  2:394140/1‑76 (MQ=255)
    gAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCttt                                                                    <  2:2537345/76‑1 (MQ=255)
       aTTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGc                                                                 <  1:1469500/76‑1 (MQ=255)
         tataCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGccc                                                               >  2:850678/1‑76 (MQ=255)
         tataCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGccc                                                               <  1:3703072/76‑1 (MQ=255)
         tataCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGccc                                                               >  1:802942/1‑76 (MQ=255)
         tataCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGccc                                                               >  1:3018072/1‑76 (MQ=255)
         tataCTCCCCGAGTCCCCTTGCCCATTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGccc                                                               >  2:6154862/1‑76 (MQ=255)
          atccTCCCCGAGTCCCCTTGCCCCTGCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGcccc                                                              <  2:2437909/73‑1 (MQ=255)
          ataCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGcccc                                                              >  1:2360323/1‑76 (MQ=255)
          ataCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGcccc                                                              >  2:505529/1‑76 (MQ=255)
          ataCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGcccc                                                              <  1:3369364/76‑1 (MQ=255)
            ccTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGt                                                            <  1:2892064/75‑1 (MQ=255)
            aCTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGt                                                            >  1:1109762/1‑76 (MQ=255)
             ccccccGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGtt                                                           <  1:1081236/74‑1 (MQ=255)
             cTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGtt                                                           <  2:3279514/76‑1 (MQ=255)
             cTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGtt                                                           <  1:2731239/76‑1 (MQ=255)
             cTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGtt                                                           >  2:3181915/1‑76 (MQ=255)
             cTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGtt                                                           <  2:5606528/76‑1 (MQ=255)
              tCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTc                                                          >  2:5408058/1‑76 (MQ=255)
               ccccGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCt                                                         <  2:4026379/76‑1 (MQ=255)
               ccccGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCt                                                         >  1:1210629/1‑76 (MQ=255)
               ccccGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCt                                                         >  1:6354307/1‑76 (MQ=255)
               ccccGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCt                                                         >  2:1462560/1‑76 (MQ=255)
                cccGATTCCCCTTGCCCCTTCGGGACACTTTTCCAAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCta                                                        <  2:5776832/76‑1 (MQ=255)
                cccGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGTGAAAAAAACGGGCCCCTTTGGCCCCGTTCta                                                        <  2:5916319/76‑1 (MQ=255)
                  cGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtctt                                                      >  2:6521976/1‑76 (MQ=255)
                  cGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtatt                                                      >  2:2839730/1‑76 (MQ=255)
                  cGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtatt                                                      >  1:3533317/1‑76 (MQ=255)
                   gAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattt                                                     <  2:4955168/76‑1 (MQ=255)
                   gAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattt                                                     <  2:2523844/76‑1 (MQ=255)
                   gAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattt                                                     >  1:2220613/1‑76 (MQ=255)
                   gAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattt                                                     <  2:4804069/76‑1 (MQ=255)
                   gAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattt                                                     >  1:6269357/1‑76 (MQ=255)
                   gAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattt                                                     >  2:5243702/1‑76 (MQ=255)
                   gAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattt                                                     <  1:4575598/76‑1 (MQ=255)
                    aGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtattta                                                    >  1:1611965/1‑76 (MQ=255)
                     gTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCtatttat                                                   <  2:6830454/76‑1 (MQ=255)
                       ccccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATttatta                                                 >  1:7305467/1‑76 (MQ=255)
                       ccccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATttatta                                                 <  2:2090565/76‑1 (MQ=255)
                        cccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAg                                                <  1:2326737/76‑1 (MQ=255)
                        cccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAg                                                <  2:5205904/76‑1 (MQ=255)
                         ccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGt                                               <  2:4581654/76‑1 (MQ=255)
                         ccTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGt                                               >  1:7343303/1‑76 (MQ=255)
                          cTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCACCTTTGGCCCCGTTCTATTTATTAGTg                                              >  2:3804757/1‑76 (MQ=255)
                           ttGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGa                                             >  1:3393801/1‑76 (MQ=255)
                           ttGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGa                                             >  2:2612898/1‑76 (MQ=255)
                             gCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGaaa                                           >  2:1976079/1‑76 (MQ=255)
                             gCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGaaa                                           >  1:6743780/1‑76 (MQ=255)
                             gCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGaaa                                           >  2:3025577/1‑76 (MQ=255)
                               cccTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAAtt                                         <  1:5386916/76‑1 (MQ=255)
                                ccTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTa                                        <  2:6298657/76‑1 (MQ=255)
                                 cTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTAc                                       >  2:3316942/1‑76 (MQ=255)
                                  ttCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACa                                      <  1:2141186/76‑1 (MQ=255)
                                  ttCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACa                                      >  2:5934175/1‑76 (MQ=255)
                                    cTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAAt                                    <  1:3803366/76‑1 (MQ=255)
                                    cTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAAt                                    >  2:1059966/1‑76 (MQ=255)
                                     tGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAAtc                                   >  1:6012855/1‑76 (MQ=255)
                                      ggACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAAtct                                  <  2:408627/76‑1 (MQ=255)
                                      ggACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAAtct                                  >  1:5631178/1‑76 (MQ=255)
                                      ggACACTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAAtct                                  >  2:1588388/1‑76 (MQ=255)
                                        acacTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCa                                <  2:1720149/76‑1 (MQ=255)
                                        acacTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCa                                >  2:2729651/1‑76 (MQ=255)
                                        acacTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCa                                <  1:3310235/76‑1 (MQ=255)
                                          acTTTTCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCAcc                              <  2:5516091/76‑1 (MQ=255)
                                            ttttCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCgg                            <  2:5122316/76‑1 (MQ=255)
                                             tttCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGa                           >  2:4926615/1‑76 (MQ=255)
                                             tttCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGa                           <  2:1597317/76‑1 (MQ=255)
                                              ttCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGAt                          <  2:4457524/76‑1 (MQ=255)
                                              ttCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGAt                          <  2:1222883/76‑1 (MQ=255)
                                               tCCGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATc                         <  2:3373022/76‑1 (MQ=255)
                                                 cGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCaa                       >  2:1814903/1‑76 (MQ=255)
                                                 cGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCaa                       >  1:3201190/1‑76 (MQ=255)
                                                 cGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCaa                       >  1:646302/1‑76 (MQ=255)
                                                 cGAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCaa                       >  1:6943755/1‑76 (MQ=255)
                                                  gAAATGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAat                      >  1:5578098/1‑76 (MQ=255)
                                                   aaaTGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAata                     >  2:650982/1‑76 (MQ=255)
                                                    aaTGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAatat                    <  2:6640629/76‑1 (MQ=255)
                                                    aaTGATGGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAatat                    <  1:5117761/76‑1 (MQ=255)
                                                     atgatgGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAacatc                   >  2:7644932/1‑75 (MQ=255)
                                                     atgatgGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATg                   >  1:4679155/1‑76 (MQ=255)
                                                      tgatgGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGc                  >  2:754225/1‑76 (MQ=255)
                                                       gatgGCGAAAAAACCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGcc                 <  1:5815895/76‑1 (MQ=255)
                                                       gatgGCGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGcc                 >  1:7107803/1‑76 (MQ=255)
                                                         tgGCAAAAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAg               <  1:4250073/76‑1 (MQ=255)
                                                            cGAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGata            <  1:4348587/76‑1 (MQ=255)
                                                             gAAAAAAACGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGatat           <  2:294843/76‑1 (MQ=255)
                                                              aaaaaaaCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATg          >  2:1150150/1‑76 (MQ=255)
                                                                 aaaaCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATc       <  2:4697037/76‑1 (MQ=255)
                                                                  aaaCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCg      >  2:291595/1‑76 (MQ=255)
                                                                  aaaCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCg      >  2:7206849/1‑76 (MQ=255)
                                                                  aaaCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCg      <  1:6539776/76‑1 (MQ=255)
                                                                   aaCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCgg     <  1:7581760/76‑1 (MQ=255)
                                                                   aaCGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCgg     >  1:7455353/1‑76 (MQ=255)
                                                                     cGGGCCCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCGGCg   >  2:1008289/1‑76 (MQ=255)
                                                                         |                                                                        
TGGCGAAATTATACTCCCCGAGTCCCCTTGCCCCTTCTGGACACTTTTCCGAAATGATGGCGAAAAAAACGGGACCCTTTGGCCCCGTTCTATTTATTAGTGAAATTACAATCTCACCGGATCAATATGCCAGATATGATCGGCGT  >  REL606/3495483‑3495628

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