fastq-dump unable to separate fastq files from scRNA-seq SRA files
1
0
Entering edit mode
13 days ago
ravihansa82 ▴ 110

Hi all, I recently downloaded a scRNA-seq file (SRR12661516) using sratoolkit. Then I needed to get the fastq files as I plan to upload the resulting fastq files to 10X genomics cloud for a new analysis. I used fastq-dump --split-3 SRR12661516 command to get the fastq files. I expected to get 3 fastq file like SRR12661516_S1_L001_I1_001.fastq, SRR12661516_S1_L001_R1_001.fastq SRR12661516_S1_L001_R2_001.fastq. Unfortunately, I got a fastq file named, SRR12661516_pass.fastq. Can someone explain why I get this output? Also, I would appreciate it if you could help me to get the required e fastq files. Thank you

SRA sratoolkit fastq-dump scRNA-seq • 270 views
ADD COMMENT
0
Entering edit mode

Your best bet is to download the original fastq files that were submitted. You will find those under the Data Access tab of each sample record. These files are in amazon cloud. Example below.

https://trace.ncbi.nlm.nih.gov/Traces/?view=run_browser&acc=SRR12661516&display=data-access

ADD REPLY
0
Entering edit mode
13 days ago
GenoMax 125k

I was able to get three expected files by doing (sra-toolkit v.3.0.1)

$ fastq-dump -F --split-files SRR12661516

$ head -4 SRR12661516*
==> SRR12661516_1.fastq <==
@NS500422:619:HM2FNBGX5:1:11101:9269:1047
GCACAATG
+NS500422:619:HM2FNBGX5:1:11101:9269:1047
AA/AAEAE

==> SRR12661516_2.fastq <==
@NS500422:619:HM2FNBGX5:1:11101:9269:1047
AAGTCNGGTAGTAGTAAAAACAGGTA
+NS500422:619:HM2FNBGX5:1:11101:9269:1047
AAAAA#EEEEEEEEEEEEEEEEEEEE

==> SRR12661516_3.fastq <==
@NS500422:619:HM2FNBGX5:1:11101:9269:1047
NACAAGTCAANTANTACAATNANATATNATACTNNTNANATGANTATTNNNAATANNN
+NS500422:619:HM2FNBGX5:1:11101:9269:1047
#AAA//////#/6#/6////#/#////#///6/##/#/#//6/#<A6/###/<</###

You will need to rename the files in the format cellranger expects (I1,R1, R2) manually.

ADD COMMENT
0
Entering edit mode

Hi GenoMax, Thanks for your help. I also found suggestions from this link and was able to get three files as you have also shown here. Yeah, I feel the same that I have to rename the sample according to cellranger format.

ADD REPLY

Login before adding your answer.

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