Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,121,304 Δ1 bp 10.0% intergenic (+86/‑30) ECB_02011 → / → ECB_02012 protein similar to bacteriophage P2 C protein/integrase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,121,3040A→.10.0% 5.1 ‑1.3 ‑0.030intergenic (+86/‑30)ECB_02011/ECB_02012protein similar to bacteriophage P2 C protein/integrase
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (16/11):  total (19/11)
Fisher's exact test strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.93e-01

TTTCAGCCACCACGAAAAAAGCGATTAGTAGTAGCAAAAAAAAGTACCACTCGGAGGGTTTTCTGATGGCA  >  REL606/2121269‑2121339
                                   |                                   
tttCAGCCACCACGAAAAAAGCGATTAGTAGTAGCa                                     <  1:2744602/36‑1 (MQ=255)
 ttCAGCCACCACGAAAAAAGCGATTAGTAGTAGCaa                                    <  1:5502886/36‑1 (MQ=255)
 ttCAGCCACCACGAAAAAAGCGATTAGTAGTAGCaa                                    <  1:4934106/36‑1 (MQ=255)
  tCATCCACCACGAAAAAAGCGATTAGTAGTAGCaaa                                   <  1:6599405/36‑1 (MQ=255)
  tCAGCCACCACGAAAAAAGCGATTAGTAGTAGCaaa                                   <  1:4515800/36‑1 (MQ=255)
   cAGCCACCACGAAAAAAGCGATTAGTAGTAGCaaaa                                  >  1:78624/1‑36 (MQ=255)
    aGCCAACACGAAAACATCGCTTAGTAGTAGCaaaaa                                 <  1:6620687/36‑1 (MQ=19)
     gCCACCACGAAAAAAGCGATTAGTAGTAGCaaaaaa                                >  1:8263665/1‑36 (MQ=255)
     gCCACCACGAAAAAAGCGATTAGTAGTAGCaaaaaa                                >  1:4888747/1‑36 (MQ=255)
      ccaccaCGAAAAAAGCGATTAGTAGTAGCaaaaaaa                               <  1:133444/36‑1 (MQ=255)
      ccacaaCGAAAAAAGCGATTAGTAGTAGCaaaaaaa                               <  1:2192238/36‑1 (MQ=255)
       caccacGAAAAAAGCGATTAGTAGTAGCaaaaaaaa                              >  1:7214972/1‑36 (MQ=255)
       caccacGAAAAAAGCGATTAGTAGTAGCaaaaaaaa                              >  1:6009264/1‑36 (MQ=255)
       caccacGAAAAAAGCGATTAGTAGTAGCaaaaaaaa                              >  1:1710750/1‑36 (MQ=255)
       caccacGAAAAAAGCGATTAGTAGTAGCaaaaaaaa                              <  1:2142965/36‑1 (MQ=255)
           acGAAAAAAGCGATTAGTAGTAGCAAAAAAAAGTAc                          >  1:1923155/1‑36 (MQ=255)
             gAAAAAAGCGATTAGTAGTAGCAAAAAAAAGTACCa                        >  1:3064698/1‑36 (MQ=255)
             gAAAAAAGCGATTAGTAGTAGCAAAAAAAAGTACCa                        >  1:7478234/1‑36 (MQ=255)
             gAAAAAAGCGATTAGTAGTAGCAAAAAAAAGTACCa                        >  1:6147477/1‑36 (MQ=255)
             gAAAAAAGCGATTAGTAGTAGCAAAAAAAAGTACCa                        <  1:1984208/36‑1 (MQ=255)
               aaaaaGCGATTAGTAGTAGCAAAAAAAAGTACCACt                      <  1:7409724/36‑1 (MQ=255)
               aaaaaGCGATTAGTAGTAGCAAAAAAAAGTACCACt                      >  1:5945594/1‑36 (MQ=255)
                aaaaGCGATTAGTAGTAGCAAAAAAAAGTACCACTc                     <  1:6976891/36‑1 (MQ=255)
                aaaaGCGATTAGCAGAAGCAAAAAAAAGTACCACTc                     >  1:4846399/1‑36 (MQ=37)
                  anGCGATTAGTAGTANCAAANAAAAGTACCACTCgg                   <  1:2461462/36‑1 (MQ=255)
                    gCGATTAGTAGTAGCAAAAAAACGTACCACTCGGAg                 <  1:7475376/36‑1 (MQ=255)
                    gCGATTAGGAGGAGCAAAAAAAAGAACCACTCggtg                 >  1:7848593/1‑34 (MQ=21)
                      gATTAGTAGTAGCAAAAAAAAGTACCACTCGGAggg               >  1:4590037/1‑36 (MQ=255)
                      gATTAGTAGTAGCAAAAAAAAGTACCACTCGGAggg               >  1:6746602/1‑36 (MQ=255)
                        ttAGTAGTAGCAAAAAAAAGTACCACTCGGAGGGtt             <  1:8245748/36‑1 (MQ=255)
                         tagtagtagCAAAAAAAAGTACCACTCGGAGGGttt            <  1:6289995/36‑1 (MQ=255)
                         tagtagtagCAAAAAAAAGTACCACTCGGAGGGttt            >  1:6610414/1‑36 (MQ=255)
                         tagtagtagC.AAAAAAAGTACCACTCGGAGGGtttt           >  1:8214720/1‑36 (MQ=38)
                          agtagtagCAAAAAAAAGTACCACTCGGAGGGtttt           <  1:3633843/36‑1 (MQ=255)
                          agtagtagC.AAAAAAAGTACCACTCGGAGGGTTTTc          >  1:2175140/1‑36 (MQ=38)
                           gtagtagCAAAAAAAAGTACCACTCGGAGGGTTTTc          >  1:6926124/1‑36 (MQ=255)
                            tagtagCAAAAAAAAGTACCACTCGGAGGGTTTTCt         >  1:1810523/1‑36 (MQ=255)
                            tagtagC.AAAAAAAGTACCACTCGGAGGGTTTTCTg        >  1:1972455/1‑36 (MQ=38)
                             agtagCAAAAAAAAGTACCACTCGGAGGGTTTTCTg        >  1:7069410/1‑36 (MQ=255)
                             agtagCAAAAAAAAGTACCACTCGGAGGGTTTTCTg        <  1:6625940/36‑1 (MQ=255)
                              gtagCAAAAAAAAGTACCACTCGGAGGGTTTTCTGa       <  1:6518267/36‑1 (MQ=255)
                              gtagCAAAAAAAAGTACCACTCGGAGGGTTTTCTGa       >  1:8315675/1‑36 (MQ=255)
                               tagCAAAAAAAAGTACCACTCGGAGGGTTTTCTGAt      >  1:4285790/1‑36 (MQ=255)
                               tagCAAAAAAAAGTACCACTCGGAGGGTTTTCTGAt      <  1:1413441/36‑1 (MQ=255)
                                 gCAAAAAAAAGTACCACTCGGAGGGTTTTCTGATgg    >  1:8431587/1‑36 (MQ=255)
                                   aaaaaaaaGTACCACTCGGAGGGTTTTCTGATGGCa  >  1:413967/1‑36 (MQ=255)
                                   |                                   
TTTCAGCCACCACGAAAAAAGCGATTAGTAGTAGCAAAAAAAAGTACCACTCGGAGGGTTTTCTGATGGCA  >  REL606/2121269‑2121339

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG