fastq-dump

[osris]

Hello:
I got paired-end reads from http://www.ncbi.nlm.nih.gov/sra/?term=ERR032965
, and decompress it with fastq dump --split-files, --split-spot .
But only Single-end appear. Why ??
Please Help me !
thanks!

sra-stat --xml --statistics ERR032965 shows this run is not paired.
<Statistics nreads="2" nspots="34127049">
<Read index="0" count="34127049" average="101" stdev="0" />
<Read index="1" count="0" average="0" stdev="0" />
</Statistics>
second read (index=1) is empty. the metadata is incorrect.

[dpryan]

They only uploaded single-end data (101 base reads). My guess is that they just clicked the wrong box when they submitted the data. The files that GenoMax linked to are the originals, so what's there is what they actually submitted (rather than what they say they submitted).

BTW, you'll need to trim that dataset as a lot of the reads are really crappy.

[transforu]

Hi dpryan and GenoMax

Thanks for your help ! and i am return result

I have tried various argument for fastq-dump, it does not work too.

GenoMax teach me download it from ebi ftp,
But it is single-end reads too. and spots less then that files from NCBI !?

Thank you very much!

[GenoMax]

ftp://ftp.sra.ebi.ac.uk/vol1/fastq/ERR032/ERR032965/

[dpryan]

Try --split-3. Also, you're guaranteed that samples with ERR... accession numbers are available in fastq format from ENA. That's usually easier than dealing with SRA.

[transforu]

Hello:
I got paired-end reads from http://www.ncbi.nlm.nih.gov/sra/?term=ERR032965
, and decompress it with fastq dump --split-files, --split-spot .
But only Single-end appear. Why ??
Please Help me !
thanks!