Trimmomatic Kseq and Bowtie --dovetail

Hello,

I am trying to analyze CutnRun data. It recommends suing Kseq tool for trimming (to trim up to 6-bp adapters from the 3′end) and Dovetail alignment in Bowtie2. Does anyone have any experience in setting this up in shell command script.

chip-seq trimmomatic

• 1,745 views • link

updated 4.8 years ago by drosophila2020 •

• written 6.3 years ago by KT •

I am now doing the same - did you get an answer to your question?

• 0 views • link

written 5.8 years ago by blur

Hello! I am on the same path here...did you get any answers? I am mostly concerned about the Kseq commands. Thank you!

• 0 views • link

written 4.8 years ago by drosophila2020 •

0 answers

No answers yet.

More posts like this

how to know what adapter sequences to trim for RNA-seq?
written 3.4 years ago by lunarskye222 •

Hello, I'm very new to RNA-seq analysis and am currently stuck on the trimming step. The libraries I am trying to analyze were built using …
Bowtie2 and subprocess issue
written 4.0 years ago by robjohn70000

Hi, I tried running bowtie2 in python using subprocess without any success. Codes: REFGENOME='/path-to-genome/genome.fasta' P1='/path-to-fastq/P1.fastq.bz2' P2='/path-to-fastq/P2.fastq.bz2' cmd=['bowtie2','--quiet','-p','10','--dovetail','--no-head','--no-sq','-x',REFGENOME,'-1',P1,'-2',P2,'-S', 'out.sam'] subprocess.Popen(cmd, shell=True) Error: No index, query, or …
Single end read adapter trimming via flexbar
written 4.7 years ago by iamsmor •

Hi guys I am beginner in sequence analysis.I am trying to trimming adapter because ı have seen that my data has adapters end of sequences …
paired-end bam clip bases outside insert range
written 5.6 years ago by pau.rodriguez.sodupe •

Hello! I have paired-end data with adapters which contain random sequences (UMIs). The problem is the 3' read trimming when the insert size is smaller …
Why homer makeTagDirectory output an empty TagDirectory
written 5.7 years ago by xiaoyonf •

I am trying to analyze Hi-C data using Homer4.11. I am stuck in the step of makeTagDirectory. After trim and bowtie2 alignment, I used the …
Bowtie not aligning shorter fragments from ATAC-seq
written 6.2 years ago by eyleen.corrales •

Hi! I'm trying to do the alignment for an ATAC-seq experiment using bowtie2. We did paired-end 100 pb sequencing. By checking with FastQC, of course …
How to process Illumina MiSeq paired-end reads for analysis
written 6.9 years ago by snow_seq •

Hello, I am trying to pre-process paired-end reads for immunoglobulin repertoire analysis. I have R1 and R2 fastq files from 2x250bp MiSeq. I am not …
Should adapters for RNASeq be removed before alignment?
written 10.3 years ago by ddzhangzz

I have paired end RNASeq Data of 8 samples (Illumina Hi - Seq). My question is whether I should trim off the adapters from these …
Exome Seq Alignment - What to do with unpaired reads post trimming
written 10.5 years ago by andrew.j.skelton73

Hi, I have some paired end exome sequencing, and after checking out the samples in fastqc, I found there was an enrichment for Illumina adapters. …
Trimming Adapters For Paired-End Sequences
written 13.5 years ago by deepthithomaskannan

<p>Hi all,</p> <p>I got illumina paired end fastq files. They told me to trim read 2 at the beginning for ~20 to 30 bp due …