This is a test version of Biostars. For the public version, visit https://www.biostars.org.
fastp error in adapter or primer trimming

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

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

0 answers

No answers yet.

Log in to answer this question.