SRA to fastq
2
0
Entering edit mode
4.4 years ago
lmobuchon ▴ 40

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 • 2.1k views
ADD COMMENT
1
Entering edit mode

Did you try without redirecting it into a new file?

ADD REPLY
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode

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

ADD REPLY
0
Entering edit mode
4.4 years ago

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.

ADD COMMENT
0
Entering edit mode
4.4 years ago

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.

ADD COMMENT

Login before adding your answer.

Traffic: 2294 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6