fastp error in adapter or primer trimming
0
0
Entering edit mode
4 months ago

For processing our 16s data two primer sequences were being used for amplification. I want to remove those two sequences with fastp. Could you tell anyone how it is possible?

forward primer -  CCTACGGGNGGCWGCAG
                          ^   ^
reverse primer -  GACTACHVGGGTATCTAATCC
                        ^^

My command is the following I was writing a script with BASH. Could anyone help me to rectify it? Thanks in advance.

fastp \
  -i ${first} \
  -I ${second} \
  -F 5 \
  -T 30 \
  -n 1 \
  -x \
  --adapter_sequence=CCTACGGG[ACGT]GGC[AT]GCAG \
  --adapter_sequence_r2=GACTAC[ACT][ACG]GGGTATCTAATCC \
  -o ${first}_trimmed.fastq \
  -O ${second}_trimmed.fastq
NGS fastp trimmomatic • 265 views
ADD COMMENT
0
Entering edit mode

Apparently cutadapt will support IUPAC characters. May want to look at that as an option.

ADD REPLY

Login before adding your answer.

Traffic: 2815 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