Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,941,586 G→A V396M (GTG→ATG)  wecF → putative enterobacterial common antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,941,5860GA99.4% 690.2 / ‑6.2 164V396M (GTG→ATG) wecFputative enterobacterial common antigen polymerase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (90/73);  total (91/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

CAACCGCGAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTC  >  REL606/3941513‑3941658
                                                                         |                                                                        
cAACCGCGAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATg                                                                        <  2:244998/76‑1 (MQ=255)
 aaCCGCGAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGc                                                                       <  2:3702383/76‑1 (MQ=255)
 aaCCGCGAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGc                                                                       <  2:4970837/76‑1 (MQ=255)
   ccGCGAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTg                                                                     <  2:2684212/76‑1 (MQ=255)
   ccGCGAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTg                                                                     <  1:978327/76‑1 (MQ=255)
    cgcgAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTgg                                                                    <  1:5015861/76‑1 (MQ=255)
      cgAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTggag                                                                  >  2:4371194/1‑74 (MQ=255)
         gaCTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGt                                                               >  2:3909897/1‑76 (MQ=255)
         gaCTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGt                                                               <  2:325668/76‑1 (MQ=255)
         gaCTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGt                                                               >  1:1364067/1‑76 (MQ=255)
          aCTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTg                                                              <  1:1700972/76‑1 (MQ=255)
           cTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGa                                                             >  2:1247612/1‑76 (MQ=255)
            tAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGaa                                                            <  1:4936749/76‑1 (MQ=255)
            tAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGaa                                                            <  1:4356055/76‑1 (MQ=255)
            tAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGaa                                                            >  2:3814648/1‑76 (MQ=255)
            tAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGaa                                                            <  2:1529706/76‑1 (MQ=255)
             aaTCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAg                                                           <  1:3090544/76‑1 (MQ=255)
             aaTCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAg                                                           <  2:695900/76‑1 (MQ=255)
              aTCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAgg                                                          <  1:1876887/76‑1 (MQ=255)
              aTCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAgg                                                          >  2:3432852/1‑76 (MQ=255)
                cGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGc                                                        >  1:2709680/1‑76 (MQ=255)
                cGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGc                                                        >  2:4031086/1‑76 (MQ=255)
                 gCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCt                                                       <  2:1486057/76‑1 (MQ=255)
                 gCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCt                                                       <  1:3728674/76‑1 (MQ=255)
                 gCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCt                                                       >  2:2838348/1‑76 (MQ=255)
                  cTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTg                                                      >  1:399729/1‑76 (MQ=255)
                  cTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTg                                                      <  2:3025130/76‑1 (MQ=255)
                  cTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTg                                                      <  2:2669298/76‑1 (MQ=255)
                  cTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTg                                                      <  1:3481290/76‑1 (MQ=255)
                  cTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTg                                                      >  1:3162498/1‑76 (MQ=255)
                  cTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTg                                                      >  1:287484/1‑76 (MQ=255)
                   tataAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTgg                                                     <  2:239962/76‑1 (MQ=255)
                   tataAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTgg                                                     >  1:616032/1‑76 (MQ=255)
                   tataAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTgg                                                     >  1:4894828/1‑76 (MQ=255)
                   tataAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTgg                                                     <  1:4320095/76‑1 (MQ=255)
                    ataAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGa                                                    <  2:3817392/76‑1 (MQ=255)
                    ataAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGa                                                    >  1:3880635/1‑76 (MQ=255)
                    ataAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGa                                                    >  2:4401655/1‑76 (MQ=255)
                     taAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGAt                                                   >  1:4158648/1‑76 (MQ=255)
                     taAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGAt                                                   >  2:4224225/1‑76 (MQ=255)
                      aaGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGGtt                                                  >  1:431314/1‑76 (MQ=255)
                       aGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTc                                                 >  1:4680721/1‑76 (MQ=255)
                       aGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTc                                                 <  1:1052940/76‑1 (MQ=255)
                       aGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTc                                                 <  1:3621516/76‑1 (MQ=255)
                       aGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTc                                                 <  1:2745273/76‑1 (MQ=255)
                          cTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGtt                                              >  1:2816884/1‑76 (MQ=255)
                          cTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGtt                                              <  2:1185525/76‑1 (MQ=255)
                          cTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGtt                                              <  2:2397785/76‑1 (MQ=255)
                           tGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGttt                                             >  1:82195/1‑76 (MQ=255)
                            gCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTg                                            <  2:2240819/76‑1 (MQ=255)
                             cGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGt                                           >  2:3482370/1‑76 (MQ=255)
                             cGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGt                                           >  2:1460366/1‑76 (MQ=255)
                               atatTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGtct                                         >  1:3455067/1‑76 (MQ=255)
                               atatTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGtct                                         <  2:5016568/76‑1 (MQ=255)
                               atatTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGtct                                         >  1:1849446/1‑76 (MQ=255)
                                tatTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGtctc                                        <  1:2272595/76‑1 (MQ=255)
                                 atTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCa                                       <  2:775152/76‑1 (MQ=255)
                                  ttACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAc                                      <  2:4508399/76‑1 (MQ=255)
                                  ttACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAc                                      <  2:4237571/76‑1 (MQ=255)
                                   tACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAcg                                     <  1:74601/76‑1 (MQ=255)
                                    acacAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAcgc                                    >  2:3081705/1‑76 (MQ=255)
                                    acacAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAcgc                                    >  2:4728655/1‑76 (MQ=255)
                                    acacAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAcgc                                    >  1:2148954/1‑76 (MQ=255)
                                    acacAGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAcgc                                    >  2:4657447/1‑76 (MQ=255)
                                     cacaGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCAcgcg                                   <  2:129401/76‑1 (MQ=255)
                                        aGTTTCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTgg                                >  2:4498483/1‑76 (MQ=255)
                                          tttCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTc                              >  1:3392659/1‑76 (MQ=255)
                                           ttCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCt                             >  1:2635064/1‑76 (MQ=255)
                                           ttCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCt                             >  1:1805413/1‑76 (MQ=255)
                                           ttCTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCt                             <  1:1736620/76‑1 (MQ=255)
                                             cTGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCttt                           >  1:958988/1‑76 (MQ=255)
                                              tGCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGGCtttt                          >  1:1326276/1‑76 (MQ=255)
                                               gCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCttttt                         >  1:2117193/1‑76 (MQ=255)
                                               gCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCttttt                         >  2:4972825/1‑76 (MQ=255)
                                               gCTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGCCttttt                         >  1:4540972/1‑76 (MQ=255)
                                                cTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCtttttt                        <  2:3647609/76‑1 (MQ=255)
                                                cTTCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCtttttt                        <  1:1104105/76‑1 (MQ=255)
                                                  tCGGGGCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTAt                      >  2:3854787/1‑76 (MQ=255)
                                                       gCGATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTgg                 >  2:2209438/1‑76 (MQ=255)
                                                         gATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTc               >  2:943843/1‑76 (MQ=255)
                                                         gATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTc               <  1:1578041/76‑1 (MQ=255)
                                                         gATCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTc               >  2:1783773/1‑76 (MQ=255)
                                                          aTCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCt              >  2:4399082/1‑76 (MQ=255)
                                                          aTCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCt              >  1:2940499/1‑76 (MQ=255)
                                                          aTCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCt              >  1:4661707/1‑76 (MQ=255)
                                                           tCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCtt             <  2:1723739/76‑1 (MQ=255)
                                                           tCTTCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCtt             <  1:1154716/76‑1 (MQ=255)
                                                             ttCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCg           <  2:150844/76‑1 (MQ=255)
                                                             ttCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCg           >  1:2564779/1‑76 (MQ=255)
                                                             ttCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCg           >  1:1243065/1‑76 (MQ=255)
                                                              tCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCgg          <  1:4783501/76‑1 (MQ=255)
                                                              tCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCgg          >  1:2893292/1‑76 (MQ=255)
                                                              tCAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCgg          <  1:1306213/76‑1 (MQ=255)
                                                               cAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGgc         >  2:360904/1‑76 (MQ=255)
                                                               cAATATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCATGGTCTTCGgc         >  2:2954857/1‑76 (MQ=255)
                                                                aaTATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGgcg        >  1:4550759/1‑76 (MQ=255)
                                                                aaTATGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGgcg        >  1:4209075/1‑76 (MQ=255)
                                                                  tatGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCTCa      <  1:2045147/76‑1 (MQ=255)
                                                                  tatGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCa      >  1:4236813/1‑76 (MQ=255)
                                                                  tatGATCATGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCa      <  2:2209089/76‑1 (MQ=255)
                                                                         |                                                                        
CAACCGCGAGACTAATCGCTATAAGGCTGCGATATTACACAGTTTCTGCTTCGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTC  >  REL606/3941513‑3941658

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