This is a test version of Biostars. For the public version, visit https://www.biostars.org.
single cell - generating cellranger counts from FASTQ files

Hi Guys,

I am trying to run cellranger count go generate count matrix for the sequenced fastq on NextSeq 500. My fastqs are named like:

SRR12345671.fastq.gz 
SRR12345672.fastq.gz 
SRR12345673.fastq.gz 
SRR12345674.fastq.gz 
SRR12345675.fastq.gz

I am using following command:

cellranger count --id=SRR1234567 --fastqs=/home/user/cellranger_count/fastqs
--sample=SRR1234567  --transcriptome=/mnt/home/user/run-cellranger/refdata-cellranger-GRCh38-3.0.0

but getting the following error as:

ERROR: No input FASTQs were found for the requested parameters.

If your files came from bcl2fastq or mkfastq:

  • Make sure you are specifying the correct --sample(s), i.e. matching the sample sheet
  • Make sure your files follow the correct naming convention, e.g. SampleName_S1_L001_R1_001.fastq.gz (and the R2 version)
  • Make sure your --fastqs points to the correct location.
  • Make sure your --lanes, if any, are correctly specified.

Refer to the "Specifying Input FASTQs" page at https://support.10xgenomics.com/ for more details.

Please help me to fix this.

10x cellranger
• 3,761 views • link
updated by seidel seidel • written by shivangi.agarwal800 shivangi.agarwal800

1 answer

See the last section here: https://support.10xgenomics.com/de-novo-assembly/software/pipelines/latest/using/fastq-input

You need to rename your files so they follow a convention cellranger expects.

[Sample Name]_S1_L00[Lane Number]_[Read Type]_001.fastq.gz

Where Read Type is one of:

    I1: Sample index read (optional)
    I2: Sample index read (optional)
    R1: Read 1
    R2: Read 2
• 0 views • link
written by GenoMax GenoMax

Hi,

Thanks. But, I see did not get the fastqs as R1, R2 and I1. There is only single fastq for each sample (5 samples = 5 FASTQ)

I expect that these are already merged. But, not sure. In this case, how to do the naming convention and how to check whether these are merged (R1, R2 and I1). Please suggest.

Thanks

• 0 views • link
written by shivangi.agarwal800 shivangi.agarwal800

Hi GenoMax,

I tried directly downloading from https://sra-pub-run-odp.s3.amazonaws.com/sra/SRR12345675/SRR12345675 given on "Data access" tab and it downloads single file. I could not get option to download R1, R2 and I1. Also, I tried to grep "R2" from the downloaded file, but could not find any. Then, I tried following commands:

fastq-dump SRR12345675

fastq-dump SRR12345675 --split-3

The second command split SRR12345675 into SRR12345675_1.fastq and SRR12345675_2.fastq. Are these original read 1 and read 2 files ? And what about I1 index file ?

Please suggest.

Thanks

• 0 views • link
updated by GenoMax GenoMax • written by shivangi.agarwal800 shivangi.agarwal800

Hopefully those are the right files. You don't need I1 file since these samples are already demultiplexed. You will need to rename the files as SRR12345675_S1_L001_R1_001.fastq and SRR12345675_S1_L001_R2_001.fastq.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this
  • cellranger count: an extremely low rate of correct barcodes was observed
    written by Alex • Alex
    I am new to cellranger. And I tried to run cellranger count for a fastq.gz file. My code is something like this: (**** is just …
  • Cellranger count with error information "... which wasn't expected, or isn't valid in this context"
    written by diqixiaoyaoer • diqixiaoyaoer
    Hello everyone, I just practiced the scRNAseq upstream analysis from cellranger official tutorial. But I met a small problem but it stopped me from going …
  • CellRanger - No input FASTQs were found for the requested parameters
    written by shuaizh117 • shuaizh117
    Hello! Can anyone help with the problem of "No input FASTQs were found for the requested parameters" using CellRanger 6.1.2? My FASTQ files are saved …
  • Cellranger count.srun issues
    written by gonzalezb549 • gonzalezb549
    Hi very new to bioinformatics and am doing an analysis of single-cell sequencing fastq data from a previously published paper. I downloaded the FASTQ files …
  • cellranger count error: No input FASTQs were found for the requested parameters.
    written by tiancaigg • tiancaigg
    `/mnt/Donghui_shared/Ann_SC_2020/X401SC21012937-Z01-F001/raw_data/JAL_1h » ls -a` JAL_1h-AK3563_S2_L001_HG2VYCCX2_I1_001.fastq.gz JAL_1h-AK3564_S3_L001_HG2VYCCX2_R2_001.fastq.gz .. JAL_1h-AK3563_S2_L001_HG2VYCCX2_R1_001.fastq.gz JAL_1h-AK3565_S4_L001_HG2VYCCX2_I1_001.fastq.gz JAL_1h-AK3562_S1_L001_HG2VYCCX2_I1_001.fastq.gz JAL_1h-AK3563_S2_L001_HG2VYCCX2_R2_001.fastq.gz JAL_1h-AK3565_S4_L001_HG2VYCCX2_R1_001.fastq.gz JAL_1h-AK3562_S1_L001_HG2VYCCX2_R1_001.fastq.gz JAL_1h-AK3564_S3_L001_HG2VYCCX2_I1_001.fastq.gz JAL_1h-AK3565_S4_L001_HG2VYCCX2_R2_001.fastq.gz JAL_1h-AK3562_S1_L001_HG2VYCCX2_R2_001.fastq.gz JAL_1h-AK3564_S3_L001_HG2VYCCX2_R1_001.fastq.gz But running cellranger count always has error: /mnt/Donghui_shared/Ann_SC_2020/X401SC21012937-Z01-F001/raw_data/JAL_1h …
  • Not able to read 10x data
    written by sidrah.maryam • sidrah.maryam
    Hello everyone, I have a few samples of 10x data which I am running using cellranger pipeline. All the samples got converted easily, however, one …
  • 10x Cell Ranger 'count' function error:
    written by miyagi • miyagi
    Dear all, I am trying to use CellRanger 'count' function on the 10x single-cell data deposited here (https://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-8077/). I tried it two ways: **1) While …
  • cellranger count help
    written by jsl • jsl
    Hello, I am trying to analyze the public dataset https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE126030 I've downloaded the fastq files onto my cluster, and would like to proceed with cellranger …
  • no matrix.mtx.gz in cellranger output
    written by yueli7 yueli7
    Hello, I followed link to run cellranger count https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/using/tutorial_ct It said in the folder `filtered_feature_bc_matrix `should have three files, I have `barcodes.tsv.gz`, `features.tsv.gz`, but I …
  • Cellranger count argument for 10xgenomic RNA+surface protein
    written by Stemcellguy • Stemcellguy
    I can't get cellranger count work. Dataset from here: https://support.10xgenomics.com/single-cell-gene-expression/datasets/3.0.0/pbmc_1k_protein_v3. The datasets has two FASTQ files, one for RNA and another for surface protein. Here …