Skip to content

Commit c227f19

Browse files
committed
add script to run EID on HIV data
This commit includes a PBS script that runs EpitopeID on the HIV data FASTQ files with the appropriate EpiDB that includes the HIV genome as a FASTA_Tag and hg19 as the main genomic sequence. Also included is an appropriate update to the `.gitignore` file.
1 parent 77b1741 commit c227f19

2 files changed

Lines changed: 27 additions & 0 deletions

File tree

paper/.gitignore

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -16,6 +16,7 @@ ENCODEdata-eGFP/results/ID
1616
HIV_samples/logs/*.err-*
1717
HIV_samples/logs/*.out-*
1818
HIV_samples/results/FASTQ
19+
HIV_samples/results/ID
1920
SyntheticDeletion/synthetic_genome/
2021
SyntheticDeletion/logs/*.err-*
2122
SyntheticDeletion/logs/*.out-*
Lines changed: 26 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,26 @@
1+
#!/bin/bash
2+
#PBS -l nodes=1:ppn=4
3+
#PBS -l pmem=16gb
4+
#PBS -l walltime=03:00:00
5+
#PBS -A open
6+
#PBS -o logs/eid.hiv.log.out
7+
#PBS -e logs/eid.hiv.log.err
8+
9+
module load gcc
10+
module load samtools
11+
module load bwa
12+
module load bedtools
13+
module load anaconda3
14+
source activate genopipe
15+
16+
WRK=/path/to/GenoPipe/paper/HIV_samples
17+
cd $WRK
18+
19+
[ -d logs ] || mkdir logs
20+
[ -d results/ID ] || mkdir -p results/ID
21+
22+
DB=$WRK/../db/hiv_EpiDB
23+
24+
EPITOPEID=$WRK/../../EpitopeID
25+
cd $EPITOPEID
26+
bash identify-Epitope.sh -i $WRK/results/FASTQ/ -o $WRK/results/ID/ -d $DB -t 4

0 commit comments

Comments
 (0)