why I download .sra file throught the SraRunInfo.csv but got files with extensions .2
1
0
Entering edit mode
3.3 years ago
mujiangxielu ▴ 10

I want to reproduce the experimental results of this article:https://www.nature.com/articles/nprot.2013.099#Abs1

but after download sra files with this code↓:

sri = read.csv("SraRunInfo.csv", stringsAsFactors=FALSE)
keep = grep("CG8144|Untreated-", sri$ LibraryName)
sri = sri[keep,]

fs = basename(sri$download_path)
for(i in 1:nrow(sri))
download.file(sri$download_path[i], fs[i])

I got some files with extensions .2 such as SRR031708.2

By the way, my SraRunInfo.csv was download by send to→file→runinfo in the NCBI page. Please tell me how to fix this problem. Thank you very much.

enter image description here

first picture is the SraRunInfo.csv in the article

second one is what in the SraRunInfo.csv I download enter image description here

NCBI Sra • 1.8k views
ADD COMMENT
0
Entering edit mode

Which SRA accession numbers are you working with?

ADD REPLY
0
Entering edit mode

SRA accession is SRP001537

ADD REPLY
0
Entering edit mode
3.3 years ago

That is how the R basename function works:

basename("a/b.2")
[1] "b.2"

from the help

basename removes all of the path up to and including the last path separator (if any).

ADD COMMENT
0
Entering edit mode

Thank you very much, but my download path in SraRunInfo.csv which download from NCBI is different from the article, I think it is the problem

ADD REPLY
0
Entering edit mode

SRA data has since moved to cloud so the links above are likely referring to the cloud data location. Those appear to be .sra format files. You can then dump reads out using fastq-dump.

ADD REPLY
0
Entering edit mode

Thank you very much!

ADD REPLY

Login before adding your answer.

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