Perfect! Thank you!
• 0 views
•
link
written
by
Darrenjdd •
How to remove duplicate reads after ONT basecalling from fastq.gz files
Hi
I ran into an error when my PC was halfway through basecalling, it errored on a fast5 file - number 577. I restarted basecalling from 578, and then when finished I basecalled 577 on its own. It looks like originally it failed after processing about 1000 reads. I can create a text file with a list of read IDs that were completed the first time, but I need something to search through the fastq.gz files for read IDs in that list and remove them. Does anyone know how I can do that?
Thanks.
• 1,467 views
•
link
written
by
Darrenjdd •
1 answer
You should be able to use filterbyname.sh from BBTools. Provide it with the names of reads you want to filter.
Log in to answer this question.
More posts like this
-
Does Guppy have no --fast_out parameters?
written by Domi •
-
Why does Guppy output different sequence data (same model)?
written by Ludwig Kian Soon •
-
Processing Fast5 files from ONT for direct RNA
written by Shagnik •
-
Bowtie2 - Snakemake
written by Mohammed •
-
BWA-MEM2 and samtools using for loop
written by rodolfo.peacewalker •
-
Tutorial: Dealing with nanopore fast5 files compressed with vbz
written by dariober
-
Read at offset 59230953264 does not begin with "@".
written by Justine •
-
CyVerse HISAT2 output
written by phamh •
-
Dealing with multiple read groups, from BWA to sorted .bam
written by jamiedm •
-
Gene_short_name missing in CuffDiff output (or how to retrieve it?)
written by nevev
