Used https://sra-explorer.info/ (SRA Explorer) to get tsv file (attached) listing SRR numbers for RNA-Seq data files associated with the experiment.

Next, we can use fasterq-dump (installed on the Charlotte HPC system) to download the files from the SRA and simultaneously convert them to fastq format.

See: this comment from IGBF-2984 "Investigate tools for detecting strandedness in RNA-Seq" for the correct options to use.

Used cut, uniq, and sort to retrieve the 29 SRR identifiers, in order, for the RNA-Seq datasets listed in the attached tsv file.

Using faster-dump, a utility part of SRA tools (released from NCBI), we can retrieve data using the following bash commands.

Dividing into 6 sections of consecutive numbers, with bash commands to run:

fasterq-dump -S -P -p SRR759123{2..4}
SRR7591232
SRR7591233
SRR7591234

fasterq-dump -S -P -p SRR759124{3..5}
SRR7591243
SRR7591244
SRR7591245

fasterq-dump -S -P -p SRR7591249
SRR7591249

fasterq-dump -S -P -p SRR75912{58..66}
SRR7591258
SRR7591259
SRR7591260
SRR7591261
SRR7591262
SRR7591263
SRR7591264
SRR7591265
SRR7591266

fasterq-dump -S -P -p SRR75912{69..73}
SRR7591269
SRR7591270
SRR7591271
SRR7591272
SRR7591273

fasterq-dump -S -P -p SRR759128{2..9}
SRR7591282
SRR7591283
SRR7591284
SRR7591285
SRR7591286
SRR7591287
SRR7591288
SRR7591289

Nowlan Freese - is the above information enough to enable you to create a first draft of the sample sheet file?

I have attached the RNA-Seq sample sheet (samplesheet_RNA-Seq.csv).
Note that the file has a header line.

A few additional notes:

• Not all samples were sequenced with paired end (we did not pay for paired end, but some samples were placed on flow cells that were running paired end samples for other clients).
• The Agami salt treated Az treated shoot samples only had two biological replicates (2 and 3). My recollection is that something had happened to the biological replicate 1 sample during library preparation and we were unable to sequence it.
• Some samples (the Az samples) were sequenced multiple times. These samples have paired end reads and an additional single end sequencing run that was done on a separate date. All of the fastq files were submitted to the SRA (the initial paired fastq files as well as the additional single end fastq file). However, it is unclear to me how those files are stored by the SRA. It could be that the SRA has combined the two R1 fastq submissions, and this is why there is a discrepancy in the spot lengths (see comment below). We may need to look at this more closely when aligning the Az samples.
• The following samples had abnormal average spot lengths according to the SRA Run Selector: SRR7591233, SRR7591234, SRR7591249, SRR7591266, SRR7591286. These samples are the Az treated samples and are the same samples that were resequenced. This may be due to an issue with the sequencing itself or it may be due to the SRA combining R1 fastq files from two files (see above comment).
• I double-checked the forward/reverse logic and it does appear correct that these samples would be labeled as reverse for strandedness.
• The read lengths all appear to be 125bp.

Adding sample sheet to the repository. Moving to Done.

Note that the sample sheet attached to this ticket is outdated. The newest version of the sample sheet can be found in the bseq_rice repository.