Question: counting the number of reads in FASTQ files
0
gravatar for Ric
10 months ago by
Ric300
Australia
Ric300 wrote:

Hi, I would like to get the read coverage and I found here C = LN / G whereas:

  • C stands for coverage
  • G is the haploid genome length
  • L is the read length
  • N is the number of reads

How is it possible to count the number of reads in FASTQ files?

Thank you in advance

sequencing coverage next-gen • 587 views
ADD COMMENTlink modified 10 months ago by swbarnes28.6k • written 10 months ago by Ric300

How to count fastq reads

on a different note, I know you didn't ask but there's software to count read depth and coverage.

ADD REPLYlink modified 10 months ago • written 10 months ago by Mark800

The same question has been asked a number of times, e.g.

How to count fastq reads

Sequence Number Count In Fastq.Gz File

ADD REPLYlink written 10 months ago by h.mon31k
1
gravatar for swbarnes2
10 months ago by
swbarnes28.6k
United States
swbarnes28.6k wrote:
wc -l mysample.fastq

or

zcat mysample.fastq.gz | wc -l

Then divide by 4.

But I strongly doubt every read in the fastq is going to align to your genome. So I don't see how this number helps.

ADD COMMENTlink written 10 months ago by swbarnes28.6k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 712 users visited in the last hour