paired-end fastq file

Hello,

I am trying to find tools to calculate the read depth of paired-end fastq files. I am trying to downsample fastq files to 5x. I know I can subset using seqtk but I want to get a specific read depth (in my case 5x) for fastq file.

sequencing

• 1,151 views • link

written 7.1 years ago by evelyn

Read-depth for a fastq file is an odd concept. I assume you mean to get 1/5th of data or 5x based on genome size?

Have you tried reformat.sh from BBMap suite. Check the sampling parameters for in-line help.

If you are trying to normalize the data (rather than just downsample) then use bbnorm.sh. A guide for that is here.

• 0 views • link

written 7.1 years ago by GenoMax

Thanks! I want to get 5x based on genome size. Like samtools depth can be used for bam/sam file, is there any tool for depth for fastq file. I have not used reformat.sh. Can you help provide the code.

• 0 views • link

written 7.1 years ago by evelyn

reformat.sh is part of BBMap suite of tools. You can download them here.

You would likely use the following parameter:

samplebasestarget=0     (sbt) Exact number of OUTPUT bases desired.

reformat.sh in1=your_R1.fq.gz in2=your_R2.fq.gz out1=sample_R1.fq.gz out2=sample_R2.fq.gz samplebasestarget=N

N = 5 x your genome size.

• 0 views • link

written 7.1 years ago by GenoMax

Thank you!!!!!!!!!

• 0 views • link

written 7.1 years ago by evelyn

0 answers

No answers yet.

More posts like this

Tools for Subsampling Paired-End Reads to 25X Coverage ?
written 1.7 years ago by Sony Nguyen •

Hi everyone, I have sequencing data for 20 rice accessions (paired-end). When I map the reads to the reference genome and calculate the coverage using …
Downsampling long-read BAM files
written 2.3 years ago by eebloom •

Im looking for the best way to downsample nanopore long-reads. Ideally I would want to downsample to a specific average read depth e.g. 5x I …
Downsampling paired-end FASTQ files with a precise number of bases (not reads) - How to do it prope…
written 6.2 years ago by mikizu •

Hi ! Firstly, I am only a Bioinformatics student, I am sorry if what I am saying isn't clear. I was asked to downsample one …
rseqc bed file for library type paired end RNA fastq files
written 6.4 years ago by evelyn

Hello, I am attempting transcript quantification using `Salmon`. I have indexed the genome. I am trying to figure out library type for the paired end …
FastQC on multiple paired end files
written 6.4 years ago by evelyn

Hello, I want to run FastQC on multiple FastQ files using an array. fastqc -o ${OUT_DIR}/${SAMPLE}.fastqc.out -f ${INPUT_DIR}/${SAMPLE}.fastq However, I could not find how to …
Base count NOT read count - Downsampling to specific number of bases not reads?
written 6.8 years ago by jnowacki

How to I down-sample to say 13 million bases total? **Quick summary:** I'm trying to get even depth of coverage and each fastq file has …
samtools depth for multiple files
written 7.1 years ago by evelyn

Can `samtools depth` be used for finding out sequence depths for individual `sorted bam` files? I am not trying to get one value for all …
sorted bam file
written 7.1 years ago by evelyn

Hi, I am trying to calculate the sequencing depth for `sorted bam` files made using `bowtie`. I am using: samtools depth bowtie_sorted.bam | awk '{sum+=$3} …
Checking bam quality
written 8.5 years ago by Gene_MMP8

I have multiple bam files from cancer datasets and I want to check their quality. What I did was, I converted the bam to fastq …
Extracting randomly subset of fastq reads from a huge file
written 9.0 years ago by KVC_bioinfo

Hello, I am trying to extract a random subset of paired end fastq reads from a huge 5Gb file. I tried using seqtk and zcat. …