Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC REL606 3,894,997 Δ6,138 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 3893584–3894993 3901134 6142–7551 56 [51] [0] 227 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 = 3894996NA (NA)133 (1.440) 64/78 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901135 = 0 (0.000)coding (1285/1428 nt) yieO predicted multidrug or homocysteine efflux system

GGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATA  >  REL606/3893533‑3893601
                                                  |                  
ggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGaa                    <  1:5870255/51‑1 (MQ=255)
ggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGaa                    <  1:10662479/51‑1 (MQ=255)
ggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGaa                    <  1:73816/51‑1 (MQ=255)
ggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGaa                    <  1:1196949/51‑1 (MQ=255)
ggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGaa                    >  1:8067677/1‑51 (MQ=255)
 gTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGaaa                   <  1:6232133/51‑1 (MQ=255)
  tCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAAt                  <  1:3402796/51‑1 (MQ=255)
  tCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAAt                  >  1:1611691/1‑51 (MQ=255)
  tCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAAt                  >  1:4764188/1‑51 (MQ=255)
  tCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAAt                  >  1:7241292/1‑51 (MQ=255)
  tCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAAt                  >  1:3884355/1‑51 (MQ=255)
   tgAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATg                 <  1:4644043/50‑1 (MQ=255)
   cGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATg                 <  1:1657444/51‑1 (MQ=255)
   cGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATg                 <  1:5284079/51‑1 (MQ=255)
    gAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATgg                >  1:8266121/1‑51 (MQ=255)
    gAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATgg                <  1:352128/51‑1 (MQ=255)
    gAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATgg                <  1:10766570/51‑1 (MQ=255)
     aaaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGa               >  1:307292/1‑51 (MQ=255)
     aaaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGa               <  1:10257247/51‑1 (MQ=255)
      aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa              <  1:2194596/51‑1 (MQ=255)
      aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa              >  1:7376165/1‑51 (MQ=255)
       aTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt             >  1:4350688/1‑51 (MQ=255)
       aTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt             >  1:5283302/1‑51 (MQ=255)
       aTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt             <  1:7582427/51‑1 (MQ=255)
       aTCTAACGCCAGACTGTACTGCACCCATTTGTTGGACGATGAAATGGAATa            <  1:9248119/51‑1 (MQ=255)
        tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2063118/51‑1 (MQ=255)
        tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            >  1:7973901/1‑51 (MQ=255)
        tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:5499956/51‑1 (MQ=255)
        tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2203008/51‑1 (MQ=255)
        tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2270214/51‑1 (MQ=255)
         cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg           >  1:4635512/1‑51 (MQ=255)
         cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg           >  1:8019158/1‑51 (MQ=255)
          tAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGc          >  1:7513431/1‑51 (MQ=255)
           aaCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2925144/48‑1 (MQ=255)
            aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        >  1:4461796/1‑51 (MQ=255)
            aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        >  1:4292986/1‑51 (MQ=255)
            aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        <  1:7191698/51‑1 (MQ=255)
            aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        <  1:7007078/51‑1 (MQ=255)
             cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc       >  1:6720716/1‑51 (MQ=255)
             cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc       >  1:6529167/1‑51 (MQ=255)
             cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc       >  1:6687293/1‑51 (MQ=255)
              gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt      <  1:7438542/51‑1 (MQ=255)
              gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt      <  1:5765465/51‑1 (MQ=255)
              gCCAGACTGTACTGCACCCATTTTGTTAGACGATGAAATGGAATAGCCCCt      <  1:7731718/51‑1 (MQ=255)
               ccAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa     >  1:4956105/1‑51 (MQ=40)
               ccAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa     >  1:4020073/1‑51 (MQ=40)
                cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa    <  1:894067/51‑1 (MQ=38)
                cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa    <  1:5529550/51‑1 (MQ=38)
                cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa    >  1:5837419/1‑51 (MQ=38)
                 aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   <  1:776138/51‑1 (MQ=35)
                 aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   <  1:2958245/51‑1 (MQ=35)
                 aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   >  1:6644406/1‑51 (MQ=35)
                 aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   <  1:10913548/51‑1 (MQ=35)
                 aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   <  1:2664828/51‑1 (MQ=35)
                  gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata  <  1:5119266/51‑1 (MQ=33)
                  gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata  <  1:2816290/51‑1 (MQ=33)
                                                  |                  
GGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATA  >  REL606/3893533‑3893601

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: