Hello!
I have a problem, we have done a PE run which was not demultiplexed automatically. I can see barcode sequences at the end of reverse reads, so it should be possible. Is there any tools to make this operation? I can make my own script, but it will be very ineffective and slow (my fastqs are 70+ gb). Thanks in advance for any help!
•
written
by
Валентин •
2 answers
Look at tools like fastp, cutadapt and similar ones.
These are more performant and probably offer more features some may provide exactly what you need.
Depends on the exact use case, for example, cutadapt has demultiplexing, though I haven't needed to use demultiplexing myself for quite a while now
https://cutadapt.readthedocs.io/en/stable/guide.html#demultiplexing
fastp has something called filter by index, but I only recall that approximately, there are many other tools actually out there
I can see barcode sequences at the end of reverse reads,
That means you are using barcodes that are integral part of the sequence reads. If that is the case then you need something like sabre: https://github.com/najoshi/sabre
Log in to answer this question.
-
How to transfer the cluster info(metadata) from seurat object to demultiplexed seurat object?
written by SoYeon.Kim •
-
Extract unique barcodes and find their frequency
written by kspata
-
How can I predict the time of calculations for RNAseq?
written by agata88
-
Converting PE to SE
written by Picasa
-
Inline barcodes in the reverse reads
written by Picasa
-
Adding barcodes to forward and reverse reads after sequencing
written by cmogren •
-
Demultiplexing fastq files with dual barcodes
written by Thomas Moody •
-
How to deal with demultiplexed Miseq pair-end (2*250bp) 16S data using QIIME?
written by nkuyfq
-
Download Pdf Format For Pubmed Articles
written by I am •
-
Are There Any Tools To Demultiplex An Sff (Or Fastq) File Which Has A Forward And Reverse Barcode?
written by Joseph Hughes
At best you ask the sequencing provider as they do this routinely and should take them little to no effort. If that is not possible please show an example of the data, so one can understand where the barcoes are. Are they in the header?Is it a single "unassigned" fastq file? There have been quite some posts here on that, please also use the search function. Demultiplexing fastq.gz files
Thank you for the reply, seems like I found partial solution with FastX toolkit barcode split, but only for the reverse reads which contains barcodes. Is there an easy way to extract forward reads from big fastq that are paired with properly demultiplexed reverse reads?
What application is this for? it's possible that the people who designed the protocol also made a demultiplexer too.