Question: Problem in processing SRA file from NCBI
gravatar for majeedaasim
21 months ago by
United States
majeedaasim40 wrote:

I downloaded the SRA file from NCBI for my organism of interest. It is Illumina sequenced paired end RNA data. Normally to create an assembly forward and reverse reads are required by Trinity. However the downloaded file has no separate fprward and reverse reads. It appears to be a merged file. The _1 and _2 suffixes suggest that.

I wonder how can I split the forward and reverse reads if it is merged. Is there any other way to get such data.


sra ncbi • 825 views
ADD COMMENTlink modified 11 months ago by gtrwst90 • written 21 months ago by majeedaasim40

How did you process the SRA file? Did you use the --split-files and -F options with fastq-dump to split the two read files and recover original Illumina fastq headers? Post the SRA # if you want someone to check on it.

ADD REPLYlink modified 21 months ago • written 21 months ago by genomax70k

I have not processed it yet. I just downloaded the SRA file through galaxy. On viewing the file it looks like this

ADD REPLYlink modified 21 months ago by genomax70k • written 21 months ago by majeedaasim40

If you are limited to working in Galaxy then I don't know the option you should use off the top of my head but make sure to choose split-files if that is available. Otherwise this is simple to take care of using from BBMap suite but that will have to be done on the command line. in=SRA.fq out1=R1.fq out2=R2.fq

ADD REPLYlink modified 21 months ago • written 21 months ago by genomax70k
gravatar for Sean Davis
21 months ago by
Sean Davis25k
National Institutes of Health, Bethesda, MD
Sean Davis25k wrote:

The --split-files argument to fastq-dump is needed. That should produce two separate files if, indeed, the original sequencing was paired-end.

ADD COMMENTlink written 21 months ago by Sean Davis25k
gravatar for gtrwst9
11 months ago by
gtrwst90 wrote:

fastq-dump --split-3 SRR...

This does NOT redownload everything.

ADD COMMENTlink written 11 months ago by gtrwst90
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1333 users visited in the last hour