Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC REL606 3,894,997 Δ6,460 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 3893577–3894996 3901456 6461–7880 51 [46] [0] 180 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)95 (0.960) 38/46 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901457 = 0 (0.000)coding (963/1428 nt) yieO predicted multidrug or homocysteine efflux system

TCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATG  >  REL606/3893541‑3893586
                                   |          
tCTAACGCCGACTGTACTGCACCCATTTTGTTGGa            >  1:7205251/1‑35 (MQ=17)
 cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa            <  1:4759200/35‑1 (MQ=255)
 cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa            <  1:9375308/35‑1 (MQ=255)
 cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa            <  1:12577486/35‑1 (MQ=255)
 cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa            <  1:7564924/35‑1 (MQ=255)
 cTAACGCCGACTGTACTGCACCCATTTTGTTGGAc           <  1:3646778/35‑1 (MQ=17)
 cTAACGCCGACTGTACTGCACCCATTTTGTTGGAc           <  1:5101731/35‑1 (MQ=17)
  tAACGCCAGACTGTACTGCACCCATTTTGTGGGAc           <  1:10046849/35‑1 (MQ=255)
   aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg          <  1:19520888/35‑1 (MQ=255)
   aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg          >  1:15502888/1‑35 (MQ=255)
   aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg          >  1:15434413/1‑35 (MQ=255)
   aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg          >  1:11718126/1‑35 (MQ=255)
   aaCGCCAGACTGTACGGCACCCATTTTGTTGGACg          >  1:11160342/1‑35 (MQ=255)
   aaCGCCGACTGTACTGCACCCATTTTGTTGGACGa         >  1:7885979/1‑35 (MQ=14)
   aaCGCCGACTGTACTGCACCCATTTTGTTGGACGa         >  1:6006317/1‑35 (MQ=14)
    aCGCCAGACTGTACTGCACCCATTTTGTTGGACGa         <  1:3365432/35‑1 (MQ=255)
    aCGCCAGACTGTACTGCACCCATTTTGTTGGACGa         >  1:678807/1‑35 (MQ=255)
    aCGCCAGACTGTACTGCACCCATTTTGTTGGACGa         >  1:7307590/1‑35 (MQ=255)
     cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt        <  1:1275610/35‑1 (MQ=255)
     cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt        >  1:19050878/1‑35 (MQ=255)
     cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt        <  1:7521486/35‑1 (MQ=255)
     cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt        <  1:15207157/35‑1 (MQ=255)
       gccgACTGTACTGCACCCATTTTGTTGGACGATGa      >  1:19070967/4‑35 (MQ=255)
       gccgACTGTACTGCACCCATTTTGTTGGACGATGa      <  1:11701875/32‑1 (MQ=255)
       gccgACTGTACTGCACCCATTTTGTTGGACGATGa      >  1:3688206/4‑35 (MQ=255)
       cccgACTGTACTGCACCCATTTTGTTGGACGATGa      <  1:11748627/32‑1 (MQ=255)
       ccAGACTGTACTGCACCCATTTTGTTGGACGATGa      >  1:11600026/1‑35 (MQ=255)
       ccAGACTGTACTGCACCCATTTTGTTGGACGATGa      >  1:11225130/1‑35 (MQ=255)
       ccAGACTGTACTGCACCCATTTTGTTGGACGATGa      <  1:13220283/35‑1 (MQ=255)
        ccgACTGTACTGCACCCATTTTGTTGGACGATGaa     >  1:18274383/3‑35 (MQ=25)
        cAGACTGTACTGCACCCATTTTGTTGGACGATGaa     >  1:3007875/1‑35 (MQ=40)
        cAGACTGTACTGCACCCATTTTGTTGGACGATGaa     >  1:18691505/1‑35 (MQ=40)
         gtaCTGTACTGCACCCATTTTGTTGGACGATGaaa    >  1:14084420/3‑35 (MQ=17)
         cgACTGTACTGCACCCATTTTGTTGGACGATGaaa    >  1:3903288/2‑35 (MQ=25)
         cgACTGTACTGCACCCATTTTGTTGGACGATGaaa    <  1:3473095/34‑1 (MQ=25)
         cgACTGTACTGCACCCATTTTGTTGGACGATGaaa    <  1:1657623/34‑1 (MQ=25)
         aaaCTGTACTGCACCCATTTTGTTGGACGATGaaa    >  1:8413227/3‑35 (MQ=17)
         aGACTGTACTGCACCCATTTTGTTGGACTATGaaa    >  1:5659756/1‑35 (MQ=255)
         aGACTGTACTGCACCCATTTTGTTGGACGATGaaa    >  1:514955/1‑35 (MQ=38)
         aGACTGTACTGCACCCATTTTGTTGGACGATGaaa    >  1:12036791/1‑35 (MQ=38)
          gACTGTACTGCACCCATTTTGTTGGACGATGAAAt   >  1:4768112/1‑35 (MQ=35)
          gACTGTACTGCACCCATTTTGTTGGACGATGAAAt   <  1:14295949/35‑1 (MQ=35)
          gACTGTACTGCACCCATTTTGTTGGACGATGAAAt   >  1:15295427/1‑35 (MQ=35)
          gACTGTACTGCACACATTTTGTTGGACGATGAAAt   <  1:5519752/35‑1 (MQ=255)
           aCTGTACTGCACCCATTTTGTTGGACGATGAAATg  <  1:993263/35‑1 (MQ=32)
           aCTGTACTGCACCCATTTTGTTGGACGATGAAATg  >  1:7352109/1‑35 (MQ=32)
           aCTGTACTGCACCCATTTTGTTGGACGATGAAATg  >  1:7309178/1‑35 (MQ=32)
           aCTGTACTGCACCCATTTTGTTGGACGATGAAATg  >  1:16830497/1‑35 (MQ=32)
           aCTGTACTGCACCCATTTTGTTGGACGATGAAATg  >  1:6324007/1‑35 (MQ=32)
           aCTGTACTGCACCCATTTTGTTGGACGATGAAATg  >  1:2436271/1‑35 (MQ=32)
           aCTGTACTGCACCCATTTTGTTGGACGATGAAATg  >  1:3816741/1‑35 (MQ=32)
                                   |          
TCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATG  >  REL606/3893541‑3893586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.