Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC REL606 3,894,997 Δ6,934 bp 100% IS150‑mediated rbsD[yieO] rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3894997 3901930 6934 52 [2] [1] 42 rbsD–[yieO] rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 3894996NA (NA)36 (0.680) 26/64 NT 97.3% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 1 (0.020)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

CTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAG..................................  >  REL606/3894960‑3894996
...................................CAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGG  >  REL606/3901931‑3901965
                                                                        
   GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGG                                   <  1:4511524/36‑1 (MQ=1)
    TGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGT                                  >  1:5941901/1‑36 (MQ=1)
    TGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGT                                  <  1:2502327/36‑1 (MQ=1)
     GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                 <  1:415321/36‑1 (MQ=1)
      TTTAACTGTCCAACTTTTTGGGGTCAGTACATGTTG                                >  1:4306690/1‑36 (MQ=38)
      TTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTG                                <  1:3296812/36‑1 (MQ=1)
      TTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTG                                >  1:7273294/1‑36 (MQ=1)
       TTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGC                               >  1:6106442/1‑36 (MQ=1)
        TAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCC                              >  1:55872/1‑36 (MQ=1)
        TAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCC                              <  1:636296/36‑1 (MQ=1)
        TAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCC                              <  1:3356077/36‑1 (MQ=1)
        TAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCC                              >  1:1668192/1‑36 (MQ=1)
           CTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAG                           <  1:1606942/36‑1 (MQ=1)
            TGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGG                          >  1:2628626/1‑36 (MQ=1)
            TGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGG                          >  1:1605999/1‑36 (MQ=1)
            TGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGG                          <  1:1390497/36‑1 (MQ=1)
            TGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGG                          <  1:4993825/36‑1 (MQ=1)
             GTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGT                         >  1:4396828/1‑36 (MQ=32)
             GTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGT                         >  1:5756695/1‑36 (MQ=32)
              TCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTG                        <  1:6739700/36‑1 (MQ=35)
                CAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGAC                      >  1:2479516/1‑36 (MQ=38)
                 AACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACC                     <  1:1571106/36‑1 (MQ=38)
                   CTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAG                   <  1:1894995/36‑1 (MQ=38)
                        TGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGC              >  1:1674917/1‑36 (MQ=38)
                         GGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCC             <  1:1074735/36‑1 (MQ=38)
                          GGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCG            <  1:5467170/36‑1 (MQ=38)
                          GGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCG            >  1:4436011/1‑36 (MQ=38)
                          GGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCG            >  1:1219124/1‑36 (MQ=38)
                            GTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCC          <  1:6077468/36‑1 (MQ=38)
                            GTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCC          <  1:7169086/36‑1 (MQ=38)
                             TCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCA         >  1:3528745/1‑36 (MQ=38)
                              CAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAA        >  1:4592622/1‑36 (MQ=38)
                              CAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAA        >  1:4504166/1‑36 (MQ=38)
                               AGTACATGTTGCCCAGGTGACCAGCACGCCGCCAAG       <  1:6438483/36‑1 (MQ=255)
                               AGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAG       >  1:5726451/1‑36 (MQ=38)
                               AGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAG       <  1:5785792/36‑1 (MQ=38)
                                                                        
CTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAG..................................  >  REL606/3894960‑3894996
...................................CAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGG  >  REL606/3901931‑3901965

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