1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
#!/bin/bash -ve
if [ -e top100k.Left.fq.gz ] && ! [ -e top100k.Left.fq ]
then
gunzip -c top100k.Left.fq.gz > top100k.Left.fq
fi
if [ -e top100k.Right.fq.gz ] && ! [ -e top100k.Right.fq ]
then
gunzip -c top100k.Right.fq.gz > top100k.Right.fq
fi
if [ -e top100k.genome.gz ] && ! [ -e top100k.genome ]
then
gunzip -c top100k.genome.gz > top100k.genome
fi
../../Trinity --left top100k.Left.fq --right top100k.Right.fq --jaccard_clip --genome top100k.genome --genome_guided_max_intron 1000 --JM 1G --seqType fq --SS_lib_type RF --output test_GG_Trinity_outdir
|
