Question: fast-dump of 10x genomics SRA data
gravatar for linjc.xmu
16 months ago by
linjc.xmu10 wrote:

Dear all, I used fastq-dump to extract single cell SRA data of 10X genomics. But the sequences ID are seemed not paired between read1 and read2. Did anyone meet with this problem?

> fastq-dump --split-files --origfmt SRR8315908.sra

Rejected 105369008 READS because READLEN < 1
Read 105369008 spots for SRR8315908.sra
Written 105369008 spots for SRR8315908.sra
> head SRR8315908_1.fastq


> head SRR8315908_2.fastq

sequence • 1.3k views
ADD COMMENTlink modified 16 months ago by genomax87k • written 16 months ago by linjc.xmu10

The shorter reads are cell barcodes/UMI and longer reads as actual sequence data.

That said ENA seems to contain only one fastq file for this accession number. That must be incorrect.

ADD REPLYlink modified 16 months ago • written 16 months ago by genomax87k

I think seq ID should be equal for one cell barcode/UMI and one sequence in read 2. Is there any way to repair the pair?

ADD REPLYlink written 16 months ago by linjc.xmu10

That is correct. You could try to rewrite the headers in file 2 programmatically to match file 1 but a manipulation like that is risky for original data.

Since you are unable to generate the original Illumina headers using the right command option you may want to contact the authors to let them know. They may not have uploaded the data properly. They may be willing to share a copy with you directly.

ADD REPLYlink modified 16 months ago • written 16 months ago by genomax87k
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: 1658 users visited in the last hour