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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3893578–3894995 3901930 6936–8353 16 [14] [1] 40 insJ‑5–[yieO] insJ‑5, insK‑5, rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 3652533NA (NA)34 (0.650) 23/66 0.9 97.1% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 1 (0.020)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA..................................  >  REL606/3652500‑3652533
................................CAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACG  >  REL606/3901931‑3901964
                                                                    
GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGG                                  <  2:92476/36‑1 (MQ=1)
GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGG                                  <  2:3215824/36‑1 (MQ=1)
 ttTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGT                                 <  2:507901/34‑1 (MQ=1)
     TAACTGTCCAACTTTTTGGGGTCGGTACAGGTTGCC                             <  2:1390980/36‑1 (MQ=255)
      AACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCC                            >  2:187745/1‑36 (MQ=38)
      AACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCC                            >  2:3742162/1‑36 (MQ=38)
      AACTGTCCAACTTTTTGCGGTCAGTACAGGTTGCCC                            >  2:295431/1‑36 (MQ=255)
       ACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCA                           <  1:1658625/36‑1 (MQ=38)
           TCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTG                       <  1:2228018/36‑1 (MQ=38)
            CCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGA                      >  2:2707943/1‑36 (MQ=38)
              AACTTTTTGGGGTCAGTACAGGTTGCCCCGGTGACC                    >  1:3417914/1‑36 (MQ=255)
               ACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCA                   >  1:2205068/1‑36 (MQ=38)
                CTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAG                  <  1:203217/36‑1 (MQ=38)
                CTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAG                  <  1:2242740/36‑1 (MQ=38)
                  TTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCA                >  1:2638320/1‑36 (MQ=255)
                  TTTTGGGGTCAGTACAGGTTGCCCAGGAGACCAGCA                >  1:3329157/1‑36 (MQ=255)
                   TTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCAC               >  1:2738077/1‑36 (MQ=255)
                    TTGGGGTCAGTATAGGTTGCCCAGGTGACCAGCACG              >  2:3717702/1‑36 (MQ=255)
                    TTGGGGTCAGTACAGGTTGCCCAGGTGACCCGCACG              >  2:1554094/1‑36 (MQ=255)
                    TTGGGGTCAGTACAGGTTGCCCAGGTGACCCGCACG              >  2:1707951/1‑36 (MQ=255)
                     TGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGC             >  2:1938686/1‑36 (MQ=255)
                      GGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCC            <  1:768231/36‑1 (MQ=255)
                      GGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCC            <  1:2078368/36‑1 (MQ=255)
                      GGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCC            <  1:1266118/36‑1 (MQ=255)
                       GGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCG           >  1:1498292/1‑36 (MQ=255)
                         GTCAGTACAGGTTGCCCAGGTGACCAGCACGCAGCC         >  1:2400219/1‑36 (MQ=255)
                          TCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCA        >  2:1848885/1‑36 (MQ=255)
                          TCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCA        <  1:1217590/36‑1 (MQ=255)
                            AGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAG      <  2:1189242/36‑1 (MQ=255)
                            AGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAG      >  1:1819558/1‑36 (MQ=255)
                              TACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAA    <  1:2226448/36‑1 (MQ=255)
                               ACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAAC   <  1:404236/36‑1 (MQ=255)
                               ACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAAC   >  2:3720465/1‑36 (MQ=255)
                               ACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAAC   <  1:1221608/36‑1 (MQ=255)
                                                                    
GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA..................................  >  REL606/3652500‑3652533
................................CAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACG  >  REL606/3901931‑3901964

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG