Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Adapter sequences not removed from short SE data #577

Open
MafGal opened this issue Sep 26, 2024 · 1 comment
Open

Adapter sequences not removed from short SE data #577

MafGal opened this issue Sep 26, 2024 · 1 comment

Comments

@MafGal
Copy link

MafGal commented Sep 26, 2024

Dear fastp team,
First of all thank you for the tool, very nice!

I am running fastp on SE data from a smallRNAseq dataset of 50nt reads, provided the default adapter sequence "TGGAATTCTCGGGTGCCAAGGC" and minimum length 15.

Although fastp finds many of the adapter sequences, it is not removing them.
Screen Shot 2024-09-26 at 11 40 31

Why is this happening / how to get files with all adapter sequences/subsequences trimmed?
Screen Shot 2024-09-26 at 00 22 22

This is my code:
`Adapter2remov="TGGAATTCTCGGGTGCCAAGGC"

threads=4

fastp -w $threads -i "$file" -o "${file/.fastq.gz/_fastpadapt30_2.fastq.gz}" -h "${file/.fastq.gz/_fastpadapt30_2.html}" -j "${file/.fastq.gz/_fastpadapt30_2.json}" --low_complexity_filter --complexity_threshold 30 --average_qual 30 --qualified_quality_phred 30 -5 -3 -r --trim_poly_x --poly_x_min_len 6 --trim_poly_g --poly_g_min_len 6 -z 4 --adapter_sequence "$Adapter2remov"`

Thanks in advance for the help,
Mafalda.

@Roleren
Copy link

Roleren commented Dec 3, 2024

Just a check, can you rerun with direct input of adapter string and "equal sign" included:

fastp -w $threads -i "$file" -o "${file/.fastq.gz/_fastpadapt30_2.fastq.gz}" -h "${file/.fastq.gz/_fastpadapt30_2.html}" -j "${file/.fastq.gz/_fastpadapt30_2.json}" --low_complexity_filter --complexity_threshold 30 --average_qual 30 --qualified_quality_phred 30 -5 -3 -r --trim_poly_x --poly_x_min_len 6 --trim_poly_g --poly_g_min_len 6 -z 4 --adapter_sequence="TGGAATTCTCGGGTGCC"`

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants