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

TCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCCTTTCTTCATAAGGGGGCGTTTTTGTT  >  REL606/3894961‑3895031
                                   |                                   
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     <  1:2718408‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:4095071‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:2337659‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     <  1:2005026‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:189383‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     <  1:4495357‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:5541653‑M1/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     <  1:1426309‑M1/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:7187602‑M1/1‑36 (MQ=255)
 cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAg                                    <  1:5154731‑M1/36‑2 (MQ=255)
 cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAg                                    <  1:2033233‑M1/36‑2 (MQ=255)
  gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAgg                                   <  1:4511524‑M1/36‑3 (MQ=255)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAggt                                  >  1:5941901‑M1/1‑33 (MQ=255)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAggt                                  <  1:2502327‑M1/36‑4 (MQ=255)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc                                  <  1:2577437/36‑1 (MQ=38)
    gtTTAACTGTCCAACTTTTTGGGGTCAGTACAggtt                                 <  1:415321‑M1/36‑5 (MQ=255)
    gtTTAACTGTCCAACTTTTTGGGGTCAGTACAGAct                                 <  1:4922983/36‑2 (MQ=25)
     tttAACTGTCCAACTTTTTGGGGTCAGTACAtgttg                                >  1:4306690‑M1/1‑31 (MQ=255)
     tttAACTGTCCAACTTTTTGGGGTCAGTACAggttg                                >  1:7273294‑M1/1‑31 (MQ=255)
     tttAACTGTCCAACTTTTTGGGGTCAGTACAggttg                                <  1:3296812‑M1/36‑6 (MQ=255)
      ttAACTGTCCAACTTTTTGGGGTCAGTACAggttgc                               >  1:6106442‑M1/1‑30 (MQ=255)
       tAACTGTCCAACTTTTTGGGGTCAGTACAggttgcc                              <  1:3356077‑M1/36‑8 (MQ=255)
       tAACTGTCCAACTTTTTGGGGTCAGTACAggttgcc                              <  1:636296‑M1/36‑8 (MQ=255)
       tAACTGTCCAACTTTTTGGGGTCAGTACAggttgcc                              >  1:55872‑M1/1‑29 (MQ=255)
       tAACTGTCCAACTTTTTGGGGTCAGTACAggttgcc                              >  1:1668192‑M1/1‑29 (MQ=255)
          cTGTCCAACTTTTTGGGGTCAGTACAggttgcccag                           <  1:1606942‑M1/36‑11 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAggttgcccagg                          <  1:1390497‑M1/36‑12 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAggttgcccagg                          >  1:1605999‑M1/1‑25 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAggttgcccagg                          <  1:4993825‑M1/36‑12 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAggttgcccagg                          >  1:2628626‑M1/1‑25 (MQ=255)
            gTCCAACTTTTTGGGGTCAGTACAggttgcccaggt                         >  1:5756695‑M1/1‑24 (MQ=255)
            gTCCAACTTTTTGGGGTCAGTACAggttgcccaggt                         >  1:4396828‑M1/1‑24 (MQ=255)
             tCCAACTTTTTGGGGTCAGTACAggttgcccaggtg                        <  1:6739700‑M1/36‑14 (MQ=255)
               cAACTTTTTGGGGTCAGTACAggttgcccaggtgac                      >  1:2479516‑M1/1‑21 (MQ=255)
                aaCTTTTTGGGGTCAGTACAggttgcccaggtgacc                     <  1:1571106‑M1/36‑17 (MQ=255)
                  cTTTTTGGGGTCAGTACAggttgcccaggtgaccag                   <  1:1894995‑M1/36‑19 (MQ=255)
                       tGGGGTCAGTACAggttgcccaggtgaccagcacgc              >  1:1674917‑M1/1‑13 (MQ=255)
                        ggggTCAGTACAggttgcccaggtgaccagcacgcc             <  1:1074735‑M1/36‑25 (MQ=255)
                         gggTCAGTACAggttgcccaggtgaccagcacgccg            <  1:5467170‑M1/36‑26 (MQ=255)
                         gggTCAGTACAggttgcccaggtgaccagcacgccg            >  1:4436011‑M1/1‑11 (MQ=255)
                         gggTCAGTACAggttgcccaggtgaccagcacgccg            >  1:1219124‑M1/1‑11 (MQ=255)
                           gTCAGTACAggttgcccaggtgaccagcacgccgcc          <  1:6077468‑M1/36‑28 (MQ=255)
                           gTCAGTACAggttgcccaggtgaccagcacgccgcc          <  1:7169086‑M1/36‑28 (MQ=255)
                            tCAGTACAggttgcccaggtgaccagcacgccgcca         >  1:3528745‑M1/1‑8 (MQ=255)
                             cAGTACAggttgcccaggtgaccagcacgccgccaa        >  1:4592622‑M1/1‑7 (MQ=255)
                             cAGTACAggttgcccaggtgaccagcacgccgccaa        >  1:4504166‑M1/1‑7 (MQ=255)
                              aGTACAtgttgcccaggtgaccagcacgccgccaag       <  1:6438483‑M1/36‑31 (MQ=255)
                              aGTACAggttgcccaggtgaccagcacgccgccaag       >  1:5726451‑M1/1‑6 (MQ=255)
                              aGTACAggttgcccaggtgaccagcacgccgccaag       <  1:5785792‑M1/36‑31 (MQ=255)
                                  cAggttgcccaggtgaccagcacgccgccaagaacg   >  1:606162‑M1/1‑2 (MQ=255)
                                   aggttgcccaggtgaccagcacgccgccaagaacgg  >  1:4132375‑M1/1‑1 (MQ=255)
                                   aggttgcccaggtgaccagcacgccgccaagaacgg  >  1:3616706‑M1/1‑1 (MQ=255)
                                   |                                   
TCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCCTTTCTTCATAAGGGGGCGTTTTTGTT  >  REL606/3894961‑3895031

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