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

How does adapter.fasta work? #7

Open
akunln opened this issue Dec 8, 2024 · 0 comments
Open

How does adapter.fasta work? #7

akunln opened this issue Dec 8, 2024 · 0 comments

Comments

@akunln
Copy link

akunln commented Dec 8, 2024

Hello,

I tried using fastplong to strip barcodes from long reads. I added all barcodes to a fasta file (adapter.fasta). Fastplong seems to ignore the --adapter_fasta as far as I could tell and still detects the adapter (here barcode+flanking sequences):
Command:
fastplong --in in.fastq.gz --out out.fastq.gz --adapter_fasta adapter.fasta ...

Terminal output:
Trying to detect adapter sequence at read start
Detected: AAGGTTAACACAAAGACACCGACAACTTTCTTCAGCACCT
Trying to detect adapter sequence at read end
Detected: AGGTGCTGAAGAAAGTTGTCGGTGTCTTTGTGTTAACCTTAGCAAT
...

Adapter.fasta:

NB01_Kit_NB96_front
CACAAAGACACCGACAACTTTCTTCAGCACCT
NB01_Kit_NB96_back
GGTGCTGAAGAAAGTTGTCGGTGTCTTTGTG
NB02_Kit_NB96_front
ACAGACGACTACAAACGGAATCGACAGCACCT
NB02_Kit_NB96_back
...

Also I would prefer if it would be possible to specify a fasta file for the -s and -e parameters so I could define multiple sequences per position. Using only the common flanking sequence in my case remove way to many bases.

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

1 participant