How to see what percentage of reads remain after quality filter to a certain phred score

Hi all,

I ran FASTQC on my raw reads to check out some quality stats. However, I would like to see how many reads will remain if i set the phred quality score to 30. Is there a way / a tool to do this?

Thanks for the help! Nikelle

phred fastqc rna-seq quality control

• 2,462 views • link

written 9.1 years ago by nikelle.petrillo

Count before/after: How to count fastq reads

• 0 views • link

written 9.1 years ago by Pierre Lindenbaum

1 answer

I'd suggest using BBDuk with various trimq thresholds; it tells you exactly how many bases are trimmed and how many reads are removed as a result:

bbduk.sh in=foo.fq out=trimmed.fq qtrim=r trimq=30

Input:                          20000 reads             3011442 bases.
QTrimmed:                       16939 reads (84.70%)    1007270 bases (33.45%)
Total Removed:                  34 reads (0.17%)        1007270 bases (33.45%)
Result:                         19966 reads (99.83%)    2004172 bases (66.55%)

• 0 views • link

written 9.1 years ago by Brian Bushnell

Omit out=filename if you just want to see what the result/stats look like.

• 0 views • link

written 9.1 years ago by GenoMax

More posts like this

How to extract reads passing a threshold for Mean Sequence Quality
written 7.9 years ago by newbinf •

I have some RNA-seq data that I trimmed and then ran through FastQC to check the quality. For the R2 file for my read, there …
EdgeR Error, 0 gene counts for 1 sample
written 10.2 years ago by nikelle.petrillo

Hi everyone, I have performed de novo transcriptome assembly through Trinity. I am now trying to perform DE in EdgeR. However, I have one sample …
Understanding Trimmomatic Sliding Window Approach
written 10.2 years ago by nikelle.petrillo

Hello all, I am performing de novo transcriptome assembly. I have used Trimmomatic to quality filter my reads. I used the argument: SLIDINGWINDOW:4:30 Can someone …
Determining read count
written 10.3 years ago by nikelle.petrillo

Hello everyone, I have run my reads through trimmomatic. I would like to know the number of reads in my raw reads and then the …
Using featurecounts for EdgeR input?
written 10.4 years ago by nikelle.petrillo

Hello, I have paired end reads and I am performing de novo RNA seq. I used tophat to align my reads to a reference transcript …
How to decipher FPKM Values?
written 10.4 years ago by nikelle.petrillo

Hello everyone, I ran count_features_given_MIN_FPKM_threshold.pl to obtain some FPKM values on my transcript. http://postimg.org/image/4tkh3pj13/ I have read through some articles in order to figure out …
How to find/decipher Trimmomatic Output
written 10.4 years ago by nikelle.petrillo

Hello all, I have ran this script. ``` java -jar trimmomatic-0.35.jar \ PE \ -phred33 \ -trimlog /home/richardsonlab/AMMA_transcripts/trimmo30trimlog.txt \ /home/richardsonlab/AMMA_transcripts/R1V22_B9F8St_ACTTGA_L002_R1_001.fastq \ /home/richardsonlab/AMMA_transcripts/R1V22_B9F8St_ACTTGA_L002_R2_001.fastq \ pairedoutput1 \ …
de novo transcriptome assembly work flow, paired end reads
written 10.5 years ago by nikelle.petrillo

Hello! I am new to RNA seq. I am trying to assemble a transcriptome, de novo. I have 100 bp paired end reads. I am …
de novo transcriptome assembly, read length?
written 10.5 years ago by nikelle.petrillo

Hello I have quality filtered my reads and I have been advised to throw out any reads with a length of 45 nucleotides of less. …
Generating Read Length?
written 10.5 years ago by nikelle.petrillo

Hello! I am new to RNA seq. I quality trimmed my fastq sequences via fastX toolkit using a phred score of 30. I would like …