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 47 [4] [2] 31 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)29 (0.530) 21/64 NT 93.5% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 2 (0.030)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

TCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCCTTTCTTCATAAGGGGGCGTTTTTG  >  REL606/3894961‑3895029
                                   |                                 
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTTCA                                   >  1:7544149‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:4083742‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   >  1:3358992‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:4575894‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:472971‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:2648630‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:2508416‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   >  1:4757615‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:5828185‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:6407107‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:6650538‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   >  1:194678‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:132634‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                   <  1:7938968‑M1/36‑1 (MQ=255)
  gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAgg                                 <  1:3124728‑M1/36‑3 (MQ=255)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAggt                                >  1:6676501‑M1/1‑33 (MQ=255)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAggt                                >  1:4333004‑M1/1‑33 (MQ=255)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc                                <  1:6376532/36‑1 (MQ=38)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc                                <  1:6723155/36‑1 (MQ=38)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc                                <  1:7355919/36‑1 (MQ=38)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc                                <  1:108894/36‑1 (MQ=38)
       tAACTGTCCGACTTTTTGGGGTCAGTACAggttgcc                            <  1:2257480‑M1/36‑8 (MQ=255)
        aaCTGTCCAACTTTTTGGGGTCAGTACAggttgccc                           <  1:4300724‑M1/36‑9 (MQ=255)
        aaCTGTCCAACTTTTTGGGGTCAGTACAggttgccc                           <  1:6885048‑M1/36‑9 (MQ=255)
         aCTGTCCAACTTTTTGGGGTCAGTACAggttgccca                          >  1:3241439‑M1/1‑27 (MQ=255)
          cTGTCCAACTTTTTGGGGTCAGTACAggttgcccag                         <  1:1501402‑M1/36‑11 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAggttgcccagg                        >  1:5260205‑M1/1‑25 (MQ=255)
            gTCCAACTTTTTGGGGTCAGTACAggttgcccaggt                       <  1:7182572‑M1/36‑13 (MQ=255)
                aaCTTTTTGGGGTCAGTACAggttgcccaggtgacc                   >  1:2295320‑M1/1‑20 (MQ=255)
                aaCTTTTTGGGGTCAGTACAggttgcccaggtgacc                   <  1:3944526‑M1/36‑17 (MQ=255)
                aaCTTTTTGGGGTCAGTACAggttgcccaggcgacc                   >  1:5208128‑M1/1‑20 (MQ=255)
                aaCTTTTTCGGGTCAGTACAggttgcccaggtgacc                   >  1:2802335‑M1/1‑20 (MQ=255)
                  cTTTTTGGGGTCAGTACAggttgcccaggtgaccag                 >  1:1798994‑M1/1‑18 (MQ=255)
                   tttttGGGGTCAGTACAggttgcccaggtgaccagc                <  1:2402067‑M1/36‑20 (MQ=255)
                     tttGGGGTCAGTACAggttgcccaggtgaccagcac              <  1:5389028‑M1/36‑22 (MQ=255)
                        ggggTCAGTACAggttgcccaggtgaccagcacgcc           <  1:5473568‑M1/36‑25 (MQ=255)
                        ggggTCAGTACAggttgcccaggtgaccagcacgcc           <  1:5727704‑M1/36‑25 (MQ=255)
                          ggTCAGTACAggttgcccaggtgaccagcacgccgc         <  1:7057901‑M1/36‑27 (MQ=255)
                           gTCAGTACAggttgcccaggtgaccagcacgccgcc        <  1:1546030‑M1/36‑28 (MQ=255)
                           gTCAGTACAggttgcccaggtgaccagcacgccgcc        >  1:7214063‑M1/1‑9 (MQ=255)
                           gTCAGTACAggttgcccaggtgaccagcacgccgcc        <  1:2322056‑M1/36‑28 (MQ=255)
                            tCAGTACAggttgcccaggtgaccagcacgccgcca       <  1:3049105‑M1/36‑29 (MQ=255)
                             cAGTACAggttgcccaggtgaccagcacgccgccaa      >  1:633399‑M1/1‑7 (MQ=255)
                             cAGTACAggttgcccaggtgaccagcacgccgccaa      <  1:5746792‑M1/36‑30 (MQ=255)
                                 aCAggttgcccaggtgaccagcacgccgccaagaac  <  1:2402568‑M1/36‑34 (MQ=255)
                                 aCAggttgcccaggtgaccagcacgccgccaagaac  <  1:2415721‑M1/36‑34 (MQ=255)
                                 aCAggttgcccaggtgaccagcacgccgccaagaac  <  1:94320‑M1/36‑34 (MQ=255)
                                   |                                 
TCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCCTTTCTTCATAAGGGGGCGTTTTTG  >  REL606/3894961‑3895029

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