Identifying adapter sequence to trim for Illumina IIx data
1
0
Entering edit mode
2.2 years ago
shpak.max ▴ 50

I have two question concerning using trimmomatic to trim adapter sequences.

The first is which adapter sequences to use. My sequence data is Illumina Genome Analyzer IIx, and I'm not sure which (if any) of the sequences is applicable:

https://github.com/timflutre/trimmomatic/tree/master/adapters

The example in the user guide is TruSeq3-PE.fa, while Illumina's website suggests AGAT...sequences (though I'm not sure if this applies to the IIx sequencing platform vs something more recent):

https://support.illumina.com/bulletins/2016/12/what-sequences-do-i-use-for-adapter-trimming.html

My second question deals with how trimmomatic handles orders of operations. If I apply the following commands, will trimmomatic first trim the adapters and then crop the first 50 remaining bases after the adapters are clipped (assuming I'm using the correct one, the one there now is just a placeholder), or do I need to first trim the adapters, save the new fastq file, and then apply trimmomatic with CROP?

for infile in *_1.fastq
    do
    base=$(basename ${infile} _1.fastq)
    ./trimmomatic PE ${infile} ${base}_2.fastq \
        ${base}_1.trim.fastq ${base}_1un.trim.fastq \
        ${base}_2.trim.fastq ${base}_2un.trim.fastq \
    ILLUMINACLIP:TruSeq3-PE.fa:2:30:10:2:True LEADING:3 TRAILING:3 MINLEN:36\
   CROP:50
  done
trimmomatic • 1.0k views
ADD COMMENT
0
Entering edit mode
ADD COMMENT
0
Entering edit mode

I saw the discussion that you linked to before posting and didn't find that that the suggestions there applied in my case. I ran FastQC and the Adapter Content section shows trace frequencies (<< 1%) for SOLID Small RNA Adapter. I'm not sure if this makes sense for the type of data set that I have, nor do I know which (if any) of the adapter sequences in this list it corresponds to:

https://github.com/timflutre/trimmomatic/tree/master/adapters

ADD REPLY
0
Entering edit mode

If it's just trace, why worry about it? It it makes no sense, and it's just trace, it's probably just some random sequence that just happens to kind of match an adapter in the database.

ADD REPLY

Login before adding your answer.

Traffic: 1869 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6