@@ -388,45 +388,10 @@ \subsubsection{Expression Counts}
388388\ begin{lstlisting}
389389mkdir -p ../../samples
390390cd ../../samples
391- cat <<EOF >samples.tsv
392- # Samples of University Toronto, HiSeq2000, paper, 10/2011
393- SRR352202 Root.full_L1 Cannabis sativa Purple Kush Root polyA RNA-Seq library PK-RT
394- SRR352200 Shoot.full_L1 Cannabis sativa Purple Kush Shoot polyA RNA-Seq library PK-SHT
395- SRR352198 Flower.pre_L1 Cannabis sativa Purple Kush Pre-Flower polyA RNA-Seq library PK-PFLW
396- SRR352196 Flower.early_L1 Cannabis sativa Purple Kush Early-stage Flower polyA RNA-Seq library PK-EFLW
397- SRR352195 Flower.mid_L1 Cannabis sativa Purple Kush Mid-stage Flower polyA RNA-Seq library PK-MFLW
398- # Samples of Michigan State, Illumina Iix, 09/2011
399- SRR306886 Leaf.mature_L1 Cannabis sativa mature leaf (fully expanded) RNA-Seq (CSA_AR)
400- SRR306885 Leaf.mature_L2 Cannabis sativa mature leaf (fully expanded) RNA-Seq (CSA_AQ)
401- SRR306884 Flower.buds_L1 Cannabis sativa flower buds RNA-Seq (CSA_AP)
402- SRR306877 Stem.primary_L1 Cannabis sativa primary stem (entire) RNA-Seq (CSA_AO)
403- SRR306875 Leaf.mature_L3 Cannabis sativa mature leaf (fully expanded) RNA-Seq (CSA_AN)
404- SRR306874 Leaf.young_L1 Cannabis sativa young leaf (<25% expanded) RNA-Seq (CSA_AM)
405- SRR306872 Leaf.young_L2 Cannabis sativa young leaf (<25% expanded) RNA-Seq (CSA_AL)
406- SRR306871 Leaf.young_L3 Cannabis sativa young leaf (<25% expanded) RNA-Seq (CSA_AK)
407- SRR306870 Flower.mature_L1 Cannabis sativa mature flower (fully expanded) RNA-Seq (CSA_AJ)
408- SRR306869 Flower.mature_L2 Cannabis sativa mature flower (fully expanded) RNA-Seq (CSA_AI)
409- SRR306868 Flower.mature_L3 Cannabis sativa mature flower (fully expanded) RNA-Seq (CSA_AH)
410- SRR306867 Flower.buds_L2 Cannabis sativa flower-buds RNA-Seq (CSA_AG)
411- SRR306866 Flower.buds_L3 Cannabis sativa flower- buds RNA-Seq (CSA_AF)
412- SRR306865 Stem.petiole_L1 Cannabis sativa stem- petioles (entire) RNA-Seq (CSA_AE)
413- SRR306864 Stem.petiole_L2 Cannabis sativa stem- petioles (entire) RNA-Seq (CSA_AD)
414- SRR306863 Root.entire_L1 Cannabis sativa entire root RNA-Seq (CSA_AC)
415- SRR306862 Root.entire_L2 Cannabis sativa entire root RNA-Seq (CSA_AB)
416- SRR306861 Root.entire_L3 Cannabis sativa entire root RNA-Seq (CSA_AA)
417- SRR192373 Stem.primary_L2 Cannabis sativa Primary Stem PE RNA-Seq (CSA_RI)
418- EOF
419-
420- # tested with sratoolkit version 2.4.4-ubuntu64
421- for i in $(grep -v "^#" samples.tsv | cut -f1)
422- do
423- /path/to/sratoolkit/bin/fastq-dump -v --split-files $1 >$1.sra.log 2>$1.sra.err
424- done
425391
392+ /path/to/sratoolkit/bin/fastq-dump *
426393\end {lstlisting }
427- This dumps the fastq files for the SRAs in the samples directory.
428- Attention: The sra files are cached in the home directory if not configured otherwise.
429- In the samples directory we now have a \texttt {\_ 1.fastq } file for each downloaded
394+ In the samples directory we now have a \texttt {.fq } file for each downloaded
430395\texttt {.sra } file. The SRA files are no longer required so you can delete them
431396to save some space. The next step is the quantification by mapping the reads
432397onto the transcripts. This quantification is done separately for each sample in
0 commit comments