...
Sample | Cell_line | GEO_Accession | Treatment | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 126 | 2.89 G | 1.55 Gb | N61311 | - | GSM1275862 | 2013-11-2616:39:00Z | GSM1275862 | Untreated | |||
2 | 126 | 2.67 G | 1.45 Gb | N61311 | - | GSM1275863 | 2013-11-2616:38:00Z | GSM1275863 | Dexamethasone | |||
3 | 126 | 2.88 G | 1.56 Gb | N61311 | - | GSM1275864 | 2013-11-2616:41:00Z | GSM1275864 | Albuterol | |||
4 | 126 | 2.76 G | 1.48 Gb | N61311 | - | GSM1275865 | 2013-11-2616:42:00Z | GSM1275865 | Albuterol_Dexamethasone | |||
5 | 126 | 3.55 G | 2.01 Gb | N052611 | - | GSM1275866 | 2013-11-2616:47:00Z | GSM1275866 | Untreated | |||
6 | 87 | 3.79 G | 2.22 Gb | N052611 | 1 | GSM1275867 | 2013-11-2617:02:00Z | GSM1275867 | Dexamethasone | |||
7 | 126 | 5.04 G | 3.03 Gb | N052611 | 1 | GSM1275868 | 2013-11-2616:49:00Z | GSM1275868 | Albuterol | |||
8 | 114 | 3.38 G | 1.86 Gb | N052611 | 1 | GSM1275869 | 2013-11-2617:03:00Z | GSM1275869 | Albuterol_Dexamethasone | |||
9 | 120 | 3.61 G | 2.10 Gb | N080611 | - | GSM1275870 | 2013-11-2617:03:00Z | GSM1275870 | Untreated | |||
10 | 126 | 4.32 G | 2.59 Gb | N080611 | - | GSM1275871 | 2013-11-2616:49:00Z | GSM1275871 | Dexamethasone | |||
11 | 126 | 3.84 G | 2.30 Gb | N080611 | - | GSM1275872 | 2013-11-2616:49:00Z | GSM1275872 | Albuterol | |||
12 | 107 | 3.38 G | 1.96 Gb | N080611 | - | GSM1275873 | 2013-11-2617:02:00Z | GSM1275873 | Albuterol_Dexamethasone | |||
13 | 101 | 3.52 G | 2.06 Gb | N061011 | 1 | GSM1275874 | 2013-11-2617:04:00Z | GSM1275874 | Untreated | |||
14 | 98 | 4.07 G | 2.39 Gb | N061011 | 1 | GSM1275875 | 2013-11-2617:08:00Z | GSM1275875 | Dexamethasone | |||
15 | 125 | 3.56 G | 1.96 Gb | N061011 | 1 | GSM1275876 | 2013-11-2617:04:00Z | GSM1275876 | Albuterol | |||
16 | 126 | 3.92 G | 2.22 Gb | N061011 | 1 | GSM1275877 | 2013-11-2616:52:00Z | GSM1275877 | Albuterol_Dexamethasone |
Fetching FASTQ files from SRA
Install the sra-tool kit if not yet installed. For example:
Code Block |
---|
#Create a conda environment
conda create -n sra -c bioconda sra-tools
conda activate sra |
Save a list of SRA accession numbers to download in a file called SraAccList.txt
Code Block |
---|
SRR1039508
SRR1039509
SRR1039510
SRR1039511
SRR1039512
SRR1039513
SRR1039514
SRR1039515
SRR1039516
SRR1039517
SRR1039518
SRR1039519
SRR1039520
SRR1039521
SRR1039522
SRR1039523 |
Use the following PBS Pro submission script to fetch the FASTQ files for all listed samples. Note, data will be downloaded to the folder where the job is submitted. Example script (fetch_SraAccList.pbs):
Code Block |
---|
#!/bin/bash -l
#PBS -N sra_fetch
#PBS -l walltime=8:00:00
#PBS -l mem=8gb
#PBS -l ncpus=4
#PBS -m bae
###PBS -M email@host
#PBS -j oe
#Usage: qsub fetch_SraAccList.pbs
cd $PBS_O_WORKDIR
for i in `cat SraAccList.txt`;
do
echo $i
prefetch $i
fastq-dump --split-files $i
done |