Question: Cutadapt list of adapter
0
gravatar for pingEde
3.7 years ago by
pingEde30
European Union
pingEde30 wrote:

Hello,

I am using Cutadapt, but I have a problem: It doesn't find adapter, but with a banal research with function "find" in the file I find adapters, so I am sure that they are in my file.

Probabily the problem is in the parameters that I wrote. In particular I have a list of adapterand I have to trim at 5'. I wrote this:

cutadapt -g adapter.fasta -o output.fastq myfile.fastq

What is the problem? Please could you help me?

Thank you in advance.

 

cutadapt trim list of adapter • 1.4k views
ADD COMMENTlink modified 3.7 years ago by Pierre Lindenbaum120k • written 3.7 years ago by pingEde30

No, I have not paired reads and I am using 454 Roche Junior

ADD REPLYlink written 3.7 years ago by pingEde30

Try with -b, to see if it finds the sequences 

-b ADAPTER, --anywhere=ADAPTER

ADD REPLYlink written 3.7 years ago by Sukhdeep Singh9.7k
0
gravatar for Brian Bushnell
3.7 years ago by
Walnut Creek, USA
Brian Bushnell16k wrote:

If you have paired reads, you can find the adapter sequences using the BBMap package like this:

bbmerge.sh in1=read1.fq in2=read2.fq outa=adapters.fa

But as long as you have standard Illumina adapters, you don't really need to do that...  they will already be in bbmap/resources/adapters.fa

ADD COMMENTlink written 3.7 years ago by Brian Bushnell16k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1034 users visited in the last hour