How to Split FASTQ File by UMI Indices into Multiple Files
Hello Biostars Community,
I have a FASTQ file that contains sequences with UMI indices in the headers. I have already counted the UMI indices and saved a list of the top 1000 UMIs in a text file named top_UMIs.txt.
I would like to split my initial FASTQ file into 1000 separate FASTQ files, where each file contains the sequences corresponding to each UMI.
Could anyone provide guidance or code to accomplish this task?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Large FASTQ file processing failure: MetaPhlAn 3.0.7 + HUMAnN on 16-core 64GB system
written by chatterjeejayeeta100 •
-
How to get host and virus UMI count per cell with Bacterial genome data?
written by San •
-
Pre-processing duplicates in SMARTer-seq data
written by kiki •
-
Bash solution to replace part of FASTA headers witch contain matching values in another file
written by clmattson •
-
Dropseq, GSE124872 and structure of the reads in the FastQ files
written by tmms •
-
BCL2Fastq when I5 index is UMI
written by ccagg
-
Split multifasta-file into multiple files according to a pattern in header
written by 6schulte •
-
fastq-dump returns I1 and R1 files instead of R1 and R2
written by C_sinensis •
-
Processing corresponding fasta records from separate multi-fasta files
written by Digsby •
-
bedtool getfasta only reads the first line of my bedfile
written by ann-katrin.llarena •
You left out a critical bit of info. Where are these UMI currently? Inside sequences, in fastq headers? Show examples if possible.
thank you GenoMax! they are in the header
Can you show an example?
It may be best to start with original data (if you moved the UMI to header) and use
reaper(LINK), which was recommended by UMI-tools author in a prior thread.Reaper works on barcodes in the sequence, not in the header.
I don't know if any tool that can do this out of the box. Might need to be a case of writing some custom code.