This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SRA: fastq-dump gives different number of sequences

I have downloaded read sequences using fastq-dump with split file option and SRR id for paired sequences. But splitted files have different number of sequence reads. As per my understanding, since these are paired-end reads these should have equal number of sequences.

$fastq-dump -I --split-files SRR390728

$grep -c '>' SRR7716545_1.fastq

694067

$grep -c '>' SRR7716545_2.fastq

1026976

Please correct me if I am wrong.

sra sequence
• 2,223 views • link
updated by Biostar Biostar • written by jeetsahu • jeetsahu

1 answer

Both files have the same number of reads. You have to grep for '^@', because @ is the fastq header prefix. > is fasta.

ls *.fastq | parallel "echo {} && grep -c '^@' {}"
SRR7716545_1.fastq
5644111
SRR7716545_2.fastq
5644111
• 1 views • link
written by ATpoint ATpoint

Thanks, I grepped different symbol. One quick question - Does fastq-dump gives latest dataset used for assembly? if yes how can I get old datasets?

• 1 views • link
written by jeetsahu • jeetsahu

fastq-dump gives the fastq based on the input SRR you give it. I have no detail knowledge about your SRR.

• 1 views • link
written by ATpoint ATpoint

Hello jeetsahu ,

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one if they work.

Upvote|Bookmark|Accept

• 1 views • link
written by finswimmer finswimmer

Log in to answer this question.

More posts like this