This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SRA to fastq

Hi everyone,

I tried to convert SRA to fastq using SRATOOLKIT 2.6.3, fastq-dump. The archive is from a ChIP-seq (single reads). This is my command line:

./fastq-dump --split-3 $DIR/SRR3722567.sra > $DIR/SRR3722567.fastq

It produces a fastq file which is empty (Written 34892751 spots for $DIR/SRR3722567.sra). I also tried the option --split-files but it produced error. How can I do ?

Thanks a lot for your help !

Lenha

sequencing
• 3,643 views • link
written by lmobuchon • lmobuchon

Did you try without redirecting it into a new file?

• 0 views • link
written by venu venu

Thanks you very much ! I used ./fastq-dump SRR3722567.sra and it works ! :)

• 0 views • link
written by lmobuchon • lmobuchon

Simply use the identifier SRR3722567. fastq-dump is quite lenient with cutting off parts of the identifier it doesn't need, though

• 0 views • link
written by Michael Dondrup Michael Dondrup

2 answers

fastq-dump duesn't write the fastq files to standard out, instead it creates files in the working dir by default (and creates cache files in you home under ncbi and doesn't remove them either). Look for a file named SRR3722567.fastq in your working dir. See the help for more information.

• 0 views • link
written by Michael Dondrup Michael Dondrup

Why are you using "--split-3"? It's a single end sequencing https://trace.ncbi.nlm.nih.gov/Traces/sra/?run=SRR3722567

Simply run 

$ fastq-dump  SRR3722567

and it should create SRR3722567.fastq in the current dir.

• 0 views • link
written by Santosh Anand Santosh Anand

Log in to answer this question.

More posts like this