Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,504,889 (A)8→7 intergenic (+633/‑92) insB‑25 → / → ECB_04162 b4563(b4576); IS1 protein insB/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,504,8820A.96.6% 1281.0 / 34.0 268intergenic (+626/‑99)insB‑25/ECB_04162b4563(b4576); IS1 protein insB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base . (109/150);  total (113/155)
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.69e-01

Only 100 of 301 total aligned reads displayed.

AGATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCAAAATA  >  REL606/4504807‑4504953
                                                                           |                                                                       
aGATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTcaca                                                                         >  1:2188520/1‑76 (MQ=255)
  aTAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaa                                                                       >  1:3181787/1‑76 (MQ=255)
   tAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaaa                                                                      >  1:10210062/1‑76 (MQ=255)
   tAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaaa                                                                      >  1:4400823/1‑76 (MQ=255)
   tAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaaa                                                                      >  2:8061263/1‑76 (MQ=255)
     aTGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaaaaa                                                                    <  1:8502505/76‑1 (MQ=255)
      tGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACaaaaaaa                                                                   <  1:3655824/76‑1 (MQ=255)
           ttAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCAtt                                                             <  1:7650953/76‑1 (MQ=255)
           ttAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCAtt                                                             <  1:2590136/76‑1 (MQ=255)
            tAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTa                                                            <  1:106591/76‑1 (MQ=255)
            tAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTa                                                            <  2:5221508/76‑1 (MQ=255)
             aGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAg                                                           >  2:1538/1‑76 (MQ=255)
             aGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAg                                                           >  2:3062388/1‑76 (MQ=255)
             aGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAg                                                           >  2:8773589/1‑76 (MQ=255)
              gTAAAAAAGAGTTGATTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGt                                                          >  1:7452201/1‑76 (MQ=255)
                aaaaaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCt                                                        >  1:5727507/1‑76 (MQ=255)
                 aaaaaGAGTTGACTGAATATAATTCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTg                                                       <  2:7725504/76‑1 (MQ=255)
                 aaaaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTg                                                       <  2:9097371/76‑1 (MQ=255)
                 aaaaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTg                                                       <  2:4932219/76‑1 (MQ=255)
                   aaaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGc                                                      >  2:5056190/1‑76 (MQ=255)
                    aaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCtt                                                    <  1:8529377/76‑1 (MQ=255)
                    aaGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCtt                                                    >  2:5277348/1‑76 (MQ=255)
                       agTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCt                                                 <  2:2671170/76‑1 (MQ=255)
                         ttGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGt                                               >  2:8473012/1‑76 (MQ=255)
                           gACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTc                                             >  1:9471348/1‑76 (MQ=255)
                           gACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTc                                             >  1:5938197/1‑76 (MQ=255)
                           gACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTc                                             >  2:5259702/1‑76 (MQ=255)
                             cTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCt                                            <  2:6122936/76‑1 (MQ=255)
                             cTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCtt                                           <  1:7199989/76‑1 (MQ=255)
                              tGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTa                                          >  1:7088092/1‑76 (MQ=255)
                              tGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTa                                          <  2:6884183/76‑1 (MQ=255)
                              tGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTa                                          <  1:1981426/76‑1 (MQ=255)
                               gAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTAc                                         <  2:9307745/76‑1 (MQ=255)
                               gAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTAc                                         <  1:13558/76‑1 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACa                                        >  2:5975451/1‑76 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACa                                        >  1:7342190/1‑76 (MQ=255)
                                aaTATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACa                                        <  1:10008348/76‑1 (MQ=255)
                                 atatAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaa                                       <  1:9909062/76‑1 (MQ=255)
                                 atatAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaa                                       <  2:3559112/76‑1 (MQ=255)
                                   ataataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaata                                     <  1:1606994/76‑1 (MQ=255)
                                   ataataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaata                                     >  1:2950550/1‑76 (MQ=255)
                                    taataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataa                                    <  1:795611/76‑1 (MQ=255)
                                    taataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataa                                    <  2:455176/76‑1 (MQ=255)
                                    taataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataa                                    >  1:7759771/1‑76 (MQ=255)
                                    taataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataa                                    >  1:3420593/1‑76 (MQ=255)
                                     aataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACaataat                                   >  1:5450222/1‑76 (MQ=255)
                                      ataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATc                                  <  2:4622523/76‑1 (MQ=255)
                                      ataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATc                                  <  2:9539667/76‑1 (MQ=255)
                                      ataCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATc                                  >  2:3206576/1‑76 (MQ=255)
                                       taCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCa                                 <  2:10091009/76‑1 (MQ=255)
                                       taCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCa                                 <  1:8908268/76‑1 (MQ=255)
                                        aCAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAg                                <  2:7158616/76‑1 (MQ=255)
                                         cAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAACCCATTAGTCTGCTTCCTGTTCTTACAATAATCAGt                               >  2:10215757/1‑76 (MQ=255)
                                         cAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGt                               <  1:6015779/76‑1 (MQ=255)
                                            ccTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatt                            <  2:3899253/76‑1 (MQ=255)
                                            ccTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatt                            <  1:2363544/76‑1 (MQ=255)
                                            ccTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatt                            <  1:9994929/76‑1 (MQ=255)
                                             ctctTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtatta                           <  2:4458334/76‑1 (MQ=255)
                                              tctTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtattat                          >  2:1456650/1‑76 (MQ=255)
                                              tctTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGtattat                          >  2:5733557/1‑76 (MQ=255)
                                               ctTCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATg                         <  2:7604936/76‑1 (MQ=255)
                                                ttCCTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGt                        >  1:6627021/1‑76 (MQ=255)
                                                  ccTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtat                      >  2:127744/1‑76 (MQ=255)
                                                  ccTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtat                      >  1:6195442/1‑76 (MQ=255)
                                                   cTACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatt                     >  2:6721500/1‑76 (MQ=255)
                                                    tACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattt                    <  2:2472103/76‑1 (MQ=255)
                                                    tACATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattt                    >  2:171682/1‑76 (MQ=255)
                                                     aCATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttt                   <  2:4300142/76‑1 (MQ=255)
                                                     aCATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttt                   <  2:1940562/76‑1 (MQ=255)
                                                     aCATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttt                   <  1:8618828/76‑1 (MQ=255)
                                                     aCATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttt                   >  1:4131327/1‑76 (MQ=255)
                                                     aCATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttt                   >  1:9661282/1‑76 (MQ=255)
                                                      cATAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtatttta                  <  2:6989199/76‑1 (MQ=255)
                                                       aTAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttat                 <  2:9463441/76‑1 (MQ=255)
                                                       aTAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttat                 >  1:3832401/1‑76 (MQ=255)
                                                       aTAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttat                 <  2:5503743/76‑1 (MQ=255)
                                                        tAAGAAGGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttatt                <  1:4636259/76‑1 (MQ=255)
                                                         aagaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttattt               >  2:1575424/1‑76 (MQ=255)
                                                           gaagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTttatttat             <  1:1137363/76‑1 (MQ=255)
                                                            aagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGtattttatt                <  1:5289245/72‑1 (MQ=255)
                                                            aagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtata            >  2:6455164/1‑76 (MQ=255)
                                                            aagGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtata            >  2:7997721/1‑76 (MQ=255)
                                                             agGAATTTGCTCACATAAAATCGATTAGTTTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTGTTtatat           >  2:3732496/1‑76 (MQ=255)
                                                             agGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtatat           <  2:4125943/76‑1 (MQ=255)
                                                             agGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtatat           <  1:3442376/76‑1 (MQ=255)
                                                             agGAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtatat           >  2:4764451/1‑76 (MQ=255)
                                                              ggAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtatata          <  1:4813591/76‑1 (MQ=255)
                                                              ggAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTtatata          <  2:9758285/76‑1 (MQ=255)
                                                               gAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAg         <  2:7248705/76‑1 (MQ=255)
                                                               gAATGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAg         <  1:5036202/76‑1 (MQ=255)
                                                                aaTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGc        >  1:6156193/1‑76 (MQ=255)
                                                                aaTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGc        <  2:6168648/76‑1 (MQ=255)
                                                                aaTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGc        <  1:9839721/76‑1 (MQ=255)
                                                                 aTGTGCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCa       <  1:3027173/76‑1 (MQ=255)
                                                                  tgtgCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCaa      <  1:9156571/76‑1 (MQ=255)
                                                                  tgtgCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCaa      <  2:2188520/76‑1 (MQ=255)
                                                                   gtgCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCaaa     <  1:499214/76‑1 (MQ=255)
                                                                   gtgCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCaaa     <  1:6951249/76‑1 (MQ=255)
                                                                   gtgCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCaaa     >  2:3462519/1‑76 (MQ=255)
                                                                    tgCTCACAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCaaaa    >  2:8533762/1‑76 (MQ=255)
                                                                           |                                                                       
AGATAATGTTTTTAGTAAAAAAGAGTTGACTGAATATAATACAGCCTCTTCCTACATAAGAAGGAATGTGCTCACAAAAAAAACCATTAGTCTGCTTCCTGTTCTTACAATAATCAGTATTATGTATTTTATTTATATAGCAAAATA  >  REL606/4504807‑4504953

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