Question: How to know which illumina adapters trimm using BCL files
0
gravatar for patrick-douglas
2.2 years ago by
Brazil/Bragança/IFPA
patrick-douglas0 wrote:

I have BCL2 files of transcriptome sequenced in the Hi Seq illumina platform, but I don't know the Illumina library prep Kit. The is any way to know which adapters to trimm directly from BCL or FASTQ files?

ADD COMMENTlink modified 2.0 years ago by chen1.9k • written 2.2 years ago by patrick-douglas0
1

Unless you have the full data folder available converting then BCL files to fastq will not be trivial. Can you ask your sequence provider to do the conversion for you? It would be significantly simpler to manage the fastq files from that point onwards.

You could use a file with all commercial adapters (e.g. adapters.fa file as supplied with BBMap suite). It may cause a bit of overtrimming but does not require you to know the exact adapter kit in advance. If you have paired end reads that overlap then it is possible to identify adapters using bbmerge.sh program from BBMap suite. FastQC contains the relevant adapter references and by running that your fastq files you may be able to identify specific illumina kit used (as long as they are present in enough numbers).

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by genomax75k

If you have the samplesheet, and the adapters are standard Illumina adapters, you can trim them on the fly while converting bcl to fastq with bcl2fastq2.

ADD REPLYlink written 2.2 years ago by h.mon28k
0
gravatar for chen
2.0 years ago by
chen1.9k
OpenGene
chen1.9k wrote:

You can use fastp to trim adapters for Illumina sequencing data, without the need of knowing the adapter sequences.

Just download fastp and run:

fastp -i in.fq -o out.fq

And then everything is done, the adapters are trimmed in out.fq

For paired end data, the command is like:

fastp -i in1.fq -o out1.fq -I in2.fq -O out2.fq

Gzip is supported for both input and output.

ADD COMMENTlink written 2.0 years ago by chen1.9k
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: 1349 users visited in the last hour