This is a test version of Biostars. For the public version, visit https://www.biostars.org.
scRNA-seq: Kallisto processing of bioproject data (fastq-dump)

Hi,

I was trying to get started with scRNA seq analysis. I downloaded a test dataset from bioproject (NCBI). However, each sample is in three fastq files (e.g. SRR123_1.fastq.gz, SRR123_2.fastq.gz, SRR123_3.fastq.gz).

How do I process these in Kallisto? Do I need to combine all of these files (how?), or split each file into forward and reverse reads before combining?

Else, what fastq-dump command do I need to issue to download the forward and reverse reads as separate files? My current command is:

fastq-dump -I --split-files SRR123

thanks!

scrna-seq rna-seq kallisto bioproject fastq-dump

could you post the command of :

gzip -dc SRR123_1.fastq.gz | head -n 4
gzip -dc SRR123_2.fastq.gz | head -n 4
gzip -dc SRR123_3.fastq.gz | head -n 4
gzip -dc SRR123_1.fastq.gz | head -n 4
@NB501328:163:HK2GVBGX5:2:11101:14937:1054
ACGAGCCANTGTACCTGTGATGGAAC
+NB501328:163:HK2GVBGX5:2:11101:14937:1054
AAAAAEEE#EEEEEEEEEEEEEEEEE

gzip -dc SRR123_2.fastq.gz | head -n 4
@NB501328:163:HK2GVBGX5:2:11101:14937:1054
TATCTAAAATNAANGTNGTNAAAAGTTATNTNNCTGTGTTNTTACNNTNNTTAANANTGTNNNATTNNNNTCCNNCANTNNTNANNNNTNNNNNNNAT
+NB501328:163:HK2GVBGX5:2:11101:14937:1054
AAAAAAAEEE#6E#EE#EE#EEEEEEEEE#E##EEEEEEE#EEEA##E##EEAE#6#EEE###/EE####EE<##/A#/##/#<####/#######/<

gzip -dc SRR123_3.fastq.gz | head -n 4
@NB501328:163:HK2GVBGX5:2:11101:14937:1054
GAAACCCT
+NB501328:163:HK2GVBGX5:2:11101:14937:1054
AAAAAEEE

ok so I guess that the SRR123_3.fastq.gz is the sequencing barcode ( to multiplex multiple samples). Could you maybe post the link to the bioproject please ?

I got the fastq files by executing: fastq-dump --split-files --gzip SRR8611970

0 answers

No answers yet.

Log in to answer this question.