download SRR*.lite.sra files using fastq-dump command
1
0
Entering edit mode
xhe764 • 0
@xhe764-11560
Last seen 8.2 years ago

I'm trying to download some chip-seq datasets from NCBI Gene Expression Omnibus using fastq-dump command like

fastq-dump SRR499718.lite.sra 

and I got error message like below:

2016-09-27T16:33:36 fastq-dump.2.7 err: item not found while constructing within virtual database module - the path 'SRR499718.lite.sra' cannot be opened as database or table.

However, running the command like following

fastq-dump SRR499718

worked fine.

When checking the GEO or SRA websites, I couldn't find any files with extension .lite.sra. My questions are 1) what are the differences between SRR*.lite.sra file and SRR*.sra file? 2) where can I download files with the extension .lite.sra? 

Any suggestions will be highly appreciated.

By the way, the problem rose when I was practising the following chip-seq workflow: 

source("http://bioconductor.org/workflows.R")
workflowInstall("chipseqDB")
sra.numbers <- c("SRR499718", "SRR499719", "SRR499720", "SRR499721", 
    "SRR499734", "SRR499735", "SRR499736", "SRR499737", "SRR499738")
grouping <- c("proB-8113", "proB-8113", "proB-8108", "proB-8108", 
    "matureB-8059", "matureB-8059", "matureB-8059", "matureB-8059", "matureB-8086")
all.sra <- paste0(sra.numbers, ".lite.sra")
data.frame(SRA=all.sra, Condition=grouping)
for (sra in all.sra) {
    code <- system(paste("fastq-dump", sra))
    stopifnot(code==0L)
}

After the for loop, I got the same error message as shown above.

 

sra chipseqDB fastq-dump .lite.sra • 3.8k views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 16 hours ago
United States

I think the lite.sra files have been phased out, so you need to do

for (sra in all.sra) {
    code <- system(paste("fastq-dump", sra.numbers))
    stopifnot(code==0L)
}

instead.

ADD COMMENT
0
Entering edit mode

Thank you so much. This explains everything.

ADD REPLY
0
Entering edit mode

I would also point out that using Aspera to download is really the way to go, and using fastq-dump directly is probably much slower than using prefetch first (which will use Aspera if installed correctly) and then using fastq-dump to convert.

ADD REPLY

Login before adding your answer.

Traffic: 916 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