# for emacs: -*- mode: sh; -*-
#########################################################################
# hivmn2 DATABASE BUILD (STARTED 04/15/08, DONE 4/21/08, Fan)

    ssh hiv1
    mkdir -p /cluster/store12/medical/hiv/hivmn2
    cd /cluster/store12/medical/hiv/hivmn2

    cd /gbdb
    mkdir hivmn2
    cd hivmn2
    cp -Rp ../hivmn1/* .

# CREATING DATABASE 

# Create the hivmn2 database.
    
    echo 'create database hivmn2' | hgsql hiv1

# CREATING GRP TABLE FOR TRACK GROUPING (DONE 4/15/08)
    
    echo "create table grp (PRIMARY KEY(NAME)) select * from hiv1.grp" \
      | hgsql hivmn2

# MAKE HGCENTRALHIV1 ENTRY AND TRACKDB TABLE FOR HIVMN2 

    echo 'insert into defaultDb values("HIV MN (GP120) V2.0", "hivmn2");' \
      | hgsql -h localhost hgcentralhiv1

    echo 'insert into dbDb values("hivmn2", "Oct. 2007", \
          "/gbdb/hivmn2/nib", "HIV MN (GP120) V2.0", "chr1", 1, 2030, \
    "HIV MN (GP120) V2.0","Human immunodeficiency     virus 1", \
    "/gbdb/hivmn2/html/description.html", 0, 0, " sequence as of Oct., 2007");' \
      | hgsql hgcentralhiv1 -h localhost

    echo 'insert into genomeClade values("HIV MN (GP120) V2.0", "other", 110);'\
      | hgsql hgcentralhiv1 -h localhost

# COPY OVER MYSQL TABLES FROM hivmn1

    Cd medical/hiv/hivmn2

# create do1 with the following lines:

echo processing table $1 ...

#hgsql hivmn2 -e "drop table ${1}"

getDbTableDef hivmn1 $1 >$1.sql
hgsql hivmn1 -N -e "select * from ${1}" >$1.tab
hgsql hivmn2 <$1.sql

hgsql hivmn2 -e "load data local infile '${1}.tab' into table ${1}"

# create doall with the following lines:

#hgsql hivmn1 -e 'drop database hivmn2'
hgsql hivmn1 -e 'create database hivmn2'
do1 aaSeq 
do1 chromInfo 
do1 dnaSeq 
do1 extFile 
do1 grp 
do1 gsIdXref 
do1 gsidClinicRec 
do1 gsidClinicRecWithSeq 
do1 gsidSubjInfo 
do1 hSeq 
do1 hgFindSpec 
do1 hgFindSpec_fanhsu 
do1 history 
do1 hivGene 
do1 interPro 
do1 seq 
do1 tableDescriptions 
do1 trackDb 
do1 trackDb_fanhsu 
do1 vax004 
do1 vax004AaCons 
do1 vax004Cons 
do1 vax004Msa

    chmod +x do*

    doall


# create trackDb

    cd kent/src/hg/makeDb/trackDb

# edit makefile to add hivmn2
    vi trackDb.ra

    cd hiv
    mkdir hivmn2
    cd hivmn2
    cp -p ../hivmn1/* .

    cd ../..
    make alpha DBS=hivmn2

# Ask admin to start BLAT server process for hivmn2 and then
# MAKE HGCENTRALHIV1 BLATSERVERS ENTRY FOR HIVMN2
    ssh hiv1
    echo 'insert into blatServers values("hivmn2", "hiv1", "17792", "1", "0"); \
          insert into blatServers values("hivmn2", "hiv1", "17793", "0", "0");' \
      | hgsql hgcentralhiv1 -h localhost

# COPY OVER MSA TABLES

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msa/AE
    cd /cluster/store12/medical/hiv/hivmn2/msa/AE

# get table definition
    mysqldump -d hivVax003Vax004 vax003AEMsa -u medcat -p$HGPSWD|hgsql hivmn2

# load the table   
    hgsql hivmn2 -e "insert into vax003AEMsa select * from hivVax003Vax004.vax003AEMsa"


# CREATE MAF TRACKS FOR VAX004

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msa
    cd /cluster/store12/medical/hiv/hivmn2/msa

# create a script file, doall

    hgsql hivmn2 -N -e \
    'select id from dnaSeq where id like "%U%"'\
    |sed -e 's/ss/do1 ss/g' >doall

# create one line script file, do1, with the following line in it:

    hgsql hivmn2 -N -e  "select id, seq from vax004Msa where id='${1}'"

    chmod +x do*

# run the script to get the .tab file with all MSA sequences of VAX004
    doall >mn2.tab
# convert .tab into .fa file
    tabToFa mn2

# grab the base alignment sequence
    echo ">hivmn2" >mn2.aln
    hgsql hivmn2 -N -e 'select seq from vax004Msa where id="MN"'  >> mn2.aln

# prepare an interium file, jjAll.mfa
    cat mn2.aln mn2.fa >jjAll.mfa
    echo = >>jjAll.mfa

# Run xmfaToMafMn1 to create a precursor file for the final .maf

    xmfaToMafMn1 jjAll.mfa j.out  org1=hivmn2
    cat j.out|sed -e 's/\./_/g'|sed -e 's/_chr/\.chr/g' >chr1.tmp

    rm jjAll.mfa j.out

    cat chr1.tmp |sed -e 's/ss_U/U/g' >chr1.maf

# copy .maf to /gbdb.

    mkdir -p  /gbdb/hivmn2/vax004Maf 
    cp chr1.maf /gbdb/hivmn2/vax004Maf -p
    echo before load
    hgLoadMaf hivmn2 vax004Maf

# create another copy for protein MAF.

    mkdir -p  /gbdb/hivmn2/vax004AaMaf 
    cp -p chr1.maf /gbdb/hivmn2/vax004AaMaf
    hgLoadMaf hivmn2 vax004AaMaf

# CREATE CONSERVATION TRACKS FOR AE STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/conservation/AE
    cd /cluster/store12/medical/hiv/hivmn2/conservation/AE

# create the .wig file and .fa file of the consensus sequence.
    gsidMsa hivmn2 vax003AEMsa MN_D533-gp120 166 vax003AECons.wig vax003AEConsensus.fa

# encode and load the wig file
    wigEncode vax003AECons.wig stdout vax003AECons.wib \
    | hgLoadWiggle hivmn2 vax003AECons stdin

# copy .wib file to /gbdb
    mkdir -p /gbdb/hivmn2/wib
    cp vax003AECons.wib /gbdb/hivmn2/wib

# do the same for protein conservation track

    mkdir aa
    cd aa

# create .wig file
    gsidAaMsa2 hivmn2 vax003AEMsa MN_D533-gp120 166 vax003AEAaCons.wig vax003AEAaConsensus.fa

# encode and load the .wib file   
    wigEncode vax003AEAaCons.wig stdout vax003AEAaCons.wib \
    | hgLoadWiggle hivmn2 vax003AEAaCons stdin

    cp vax003AEAaCons.wib /gbdb/hivmn2/wib

# CREATE MAF TRACKS FOR AE STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msa/AE
    cd /cluster/store12/medical/hiv/hivmn2/msa/AE

# create a script file, doall

    hgsql hivmn2 -N -e \
    'select id from dnaSeq where id like "%T%"'\
    |sed -e 's/ss/do1 ss/g' >doall

# create one line script file, do1, with the following line in it:

    hgsql hivmn2 -N -e  "select id, seq from vax003AEMsa where id='${1}'"

    chmod +x do*

# run the script to get the .tab file with all MSA sequences of VAX004
    doall >mn2.tab
# convert .tab into .fa file
    tabToFa mn2

# grab the base alignment sequence
    echo ">hivmn2" >mn2.aln
    hgsql hivmn2 -N -e 'select seq from vax003AEMsa where id="MN_D533-gp120"'  >> mn2.aln

# prepare an interium file, jjAll.mfa
    cat mn2.aln mn2.fa >jjAll.mfa
    echo = >>jjAll.mfa

# Run xmfaToMafMn2AE to create a precursor file for the final .maf

    xmfaToMafMn2AE jjAll.mfa j.out  org1=hivmn2
    cat j.out|sed -e 's/\./_/g'|sed -e 's/_chr/\.chr/g' >chr1.tmp

    rm jjAll.mfa j.out

    cat chr1.tmp |sed -e 's/ss_T/T/g' >chr1.maf

# copy .maf to /gbdb.

    mkdir -p  /gbdb/hivmn2/vax003AEMaf
    cp chr1.maf /gbdb/hivmn2/vax003AEMaf -p
    echo before load
    hgLoadMaf hivmn2 vax003AEMaf

# create another copy for protein MAF.

    mkdir -p  /gbdb/hivmn2/vax003AEAaMaf 
    cp -p chr1.maf /gbdb/hivmn2/vax003AEAaMaf
    hgLoadMaf hivmn2 vax003AEAaMaf



# COPY OVER MSA TABLES FOR B STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msa/B
    cd /cluster/store12/medical/hiv/hivmn2/msa/B

# get table definition
    mysqldump -d hivVax003Vax004 vax003BMsa -u medcat -p$HGPSWD|hgsql hivmn2

# load the table   
    hgsql hivmn2 -e "insert into vax003BMsa select * from hivVax003Vax004.vax003BMsa"

# CREATE CONSERVATION TRACKS FOR B STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/conservation/B
    cd /cluster/store12/medical/hiv/hivmn2/conservation/B

# create the .wig file and .fa file of the consensus sequence.
    gsidMsa hivmn2 vax003BMsa MN_D533-gp120 166 vax003BCons.wig vax003BConsensus.fa

# encode and load the wig file
    wigEncode vax003BCons.wig stdout vax003BCons.wib \
    | hgLoadWiggle hivmn2 vax003BCons stdin

# copy .wib file to /gbdb
    mkdir -p /gbdb/hivmn2/wib
    cp vax003BCons.wib /gbdb/hivmn2/wib

# do the same for protein conservation track

    mkdir aa
    cd aa

# create .wig file
    gsidAaMsa2 hivmn2 vax003BMsa MN_D533-gp120 166 vax003BAaCons.wig vax003BAaConsensus.fa

# encode and load the .wib file   
    wigEncode vax003BAaCons.wig stdout vax003BAaCons.wib \
    | hgLoadWiggle hivmn2 vax003BAaCons stdin

    cp vax003BAaCons.wib /gbdb/hivmn2/wib

# CREATE MAF TRACKS FOR B STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msa/B
    cd /cluster/store12/medical/hiv/hivmn2/msa/B

# create a script file, doall

    hgsql hivmn2 -N -e \
    'select id from dnaSeq where id like "%T%"'\
    |sed -e 's/ss/do1 ss/g' >doall

# create one line script file, do1, with the following line in it:

    hgsql hivmn2 -N -e  "select id, seq from vax003BMsa where id='${1}'"

    chmod +x do*

# run the script to get the .tab file with all MSA sequences of VAX004
    doall >mn2.tab
# convert .tab into .fa file
    tabToFa mn2

# grab the base alignment sequence
    echo ">hivmn2" >mn2.aln
    hgsql hivmn2 -N -e 'select seq from vax003BMsa where id="MN_D533-gp120"'  >> mn2.aln

# prepare an interium file, jjAll.mfa
    cat mn2.aln mn2.fa >jjAll.mfa
    echo = >>jjAll.mfa

# Run xmfaToMafMn2B to create a precursor file for the final .maf

    xmfaToMafMn1 jjAll.mfa j.out  org1=hivmn2
    cat j.out|sed -e 's/\./_/g'|sed -e 's/_chr/\.chr/g' >chr1.tmp

    rm jjAll.mfa j.out

    cat chr1.tmp |sed -e 's/ss_T/T/g' >chr1.maf

# copy .maf to /gbdb.

    mkdir -p  /gbdb/hivmn2/vax003BMaf
    cp chr1.maf /gbdb/hivmn2/vax003BMaf -p
    
    hgLoadMaf hivmn2 vax003BMaf

# create another copy for protein MAF.

    mkdir -p  /gbdb/hivmn2/vax003BAaMaf 
    cp -p chr1.maf /gbdb/hivmn2/vax003BAaMaf
    hgLoadMaf hivmn2 vax003BAaMaf

#########################################################################
# CREATE VAX003 TRACK	(Done 5/13/08, Fan) 

    cd ~/medical/hiv/hivmn2

# get vax003 sequences
    hgsql hivmn2 -N -e 'select * from dnaSeq where id like "%T%"' >vax003.tab

# create .fa file
    tabToFa vax003

    mkdir -p /gbdb/hivmn2/vax003
    cp -p vax003.fa /gbdb/hivmn2/vax003/vax003.fa

    hgLoadSeq -replace hivmn2 /gbdb/hivmn2/vax003/vax003.fa

# BLAT
    gfClient -minScore=200 -minIdentity=70 -nohead hiv1.soe.ucsc.edu 17793 \
    /gbdb/hivmn2/nib -out=psl -t=dna -q=dna vax003.fa vax003.psl

# count the result
    wc *.psl
    cut -f 10 vax003.psl |wc
    cut -f 10 vax003.psl |sort -u |wc

# load the psl result into vax003 table
    hgLoadPsl hivmn2 vax003.psl

# hgLoadPsl has some file permission problem.  Finish this by manually load the psl.tab file.
    hgsql hivmn2 -e 'load data local infile "psl.tab" into table vax003'

#########################################################################
# Build the gsidSubjSeq table (used by Table View).

   gsidSubjSeq hivmn2 dnaSeqId > j.dna
   gsidSubjSeq hivmn2 aaSeqId > j.aa

   cut -f 1 j.dna >j.1
   cut -f 1 j.aa  >j.2

   cut -f 2 j.dna  >j.3
   cut -f 2 j.aa   >j.4

   paste j.1 j.3 j.4> gsidSubjSeq.tab

   hgsql hivmn2 -e 'delete from gsidSubjSeq'
   hgsql hivmn2 -e \
   'load data local infile "gsidSubjSeq.tab" into table gsidSubjSeq'

   rm j.1 j.2 j.3 j.4 j.dna j.aa

#########################################################################
# RE-BUILD CONSERVATION AND MSA TRACKS FOR vax003AE

# COPY OVER shortened MSA sequences from hivVax003Vax003

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msaNew/AE
    cd /cluster/store12/medical/hiv/hivmn2/msaNew/AE

    hgsql hivmn2 -e "delete from vax003AEMsa"
    hgsql hivmn2 -e "insert into vax003AEMsa select * from hivVax003Vax003.vax003AEMsa"

# CREATE CONSERVATION TRACKS FOR AE STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/conservationNew/AE
    cd /cluster/store12/medical/hiv/hivmn2/conservationNew/AE

# create the .wig file and .fa file of the consensus sequence.
    gsidMsa hivmn2 vax003AEMsa MN_D533-gp120 166 vax003AECons.wig vax003AEConsensus.fa

# encode and load the wig file
    wigEncode vax003AECons.wig stdout vax003AECons.wib \
    | hgLoadWiggle hivmn2 vax003AECons stdin

# copy .wib file to /gbdb
    mkdir -p /gbdb/hivmn2/wib
    cp vax003AECons.wib /gbdb/hivmn2/wib

# do the same for protein conservationNew track

    mkdir aa
    cd aa

# create .wig file
    gsidAaMsa2 hivmn2 vax003AEMsa MN_D533-gp120 166 vax003AEAaCons.wig vax003AEAaConsensus.fa

# encode and load the .wib file   
    wigEncode vax003AEAaCons.wig stdout vax003AEAaCons.wib \
    | hgLoadWiggle hivmn2 vax003AEAaCons stdin

    cp vax003AEAaCons.wib /gbdb/hivmn2/wib

# CREATE MAF TRACKS FOR vax003AE STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msaNew/AE
    cd /cluster/store12/medical/hiv/hivmn2/msaNew/AE

# create a script file, doall

    hgsql hivmn2 -N -e \
    'select id from vax003AEMsa where id like "%T%"'\
    |sed -e 's/ss/do1 ss/g' >doall

# create one line script file, do1, with the following line in it:

    hgsql hivmn2 -N -e  "select id, seq from vax003AEMsa where id='${1}'"

    chmod +x do*

# run the script to get the .tab file with all MSA sequences of VAX004
    doall >mn2.tab
# convert .tab into .fa file
    tabToFa mn2

# grab the base alignment sequence
    echo ">hivmn2" >mn2.aln
    hgsql hivmn2 -N -e 'select seq from vax003AEMsa where id="MN_D533-gp120"'  >> mn2.aln

# prepare an interium file, jjAll.mfa
    cat mn2.aln mn2.fa >jjAll.mfa
    echo = >>jjAll.mfa

# Run xmfaToMafMn2AE to create a precursor file for the final .maf

    xmfaToMafMn2AE jjAll.mfa j.out  org1=hivmn2
    cat j.out|sed -e 's/\./_/g'|sed -e 's/_chr/\.chr/g' >chr1.tmp

    rm jjAll.mfa j.out

    cat chr1.tmp |sed -e 's/ss_T/T/g' >chr1.maf

# copy .maf to /gbdb.

    mkdir -p  /gbdb/hivmn2/vax003AEMaf
    cp chr1.maf /gbdb/hivmn2/vax003AEMaf -p
    echo before load
    hgLoadMaf hivmn2 vax003AEMaf

# create another copy for protein MAF.

    mkdir -p  /gbdb/hivmn2/vax003AEAaMaf 
    cp -p chr1.maf /gbdb/hivmn2/vax003AEAaMaf
    hgLoadMaf hivmn2 vax003AEAaMaf

#########################################################################
# RE-BUILD CONSERVATION AND MSA TRACKS FOR vax003B

# COPY OVER shortened MSA sequences from hivVax003Vax004

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msaNew/B
    cd /cluster/store12/medical/hiv/hivmn2/msaNew/B

    hgsql hivmn2 -e "delete from vax003BMsa"
    hgsql hivmn2 -e "insert into vax003BMsa select * from hivVax003Vax004.vax003BMsa"

# CREATE CONSERVATION TRACKS FOR B STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/conservationNew/B
    cd /cluster/store12/medical/hiv/hivmn2/conservationNew/B

# create the .wig file and .fa file of the consensus sequence.
    gsidMsa hivmn2 vax003BMsa MN_D533-gp120 166 vax003BCons.wig vax003BConsensus.fa

# encode and load the wig file
    wigEncode vax003BCons.wig stdout vax003BCons.wib \
    | hgLoadWiggle hivmn2 vax003BCons stdin

# copy .wib file to /gbdb
    mkdir -p /gbdb/hivmn2/wib
    cp vax003BCons.wib /gbdb/hivmn2/wib

# do the same for protein conservationNew track

    mkdir aa
    cd aa

# create .wig file
    gsidAaMsa2 hivmn2 vax003BMsa MN_D533-gp120 166 vax003BAaCons.wig vax003BAaConsensus.fa

# encode and load the .wib file   
    wigEncode vax003BAaCons.wig stdout vax003BAaCons.wib \
    | hgLoadWiggle hivmn2 vax003BAaCons stdin

    cp vax003BAaCons.wib /gbdb/hivmn2/wib

# CREATE MAF TRACKS FOR vax003B STRAIN

    mkdir -p /cluster/store12/medical/hiv/hivmn2/msaNew/B
    cd /cluster/store12/medical/hiv/hivmn2/msaNew/B

# create a script file, doall

    hgsql hivmn2 -N -e \
    'select id from vax003BMsa where id like "%T%"'\
    |sed -e 's/ss/do1 ss/g' >doall

# create one line script file, do1, with the following line in it:

    hgsql hivmn2 -N -e  "select id, seq from vax003BMsa where id='${1}'"

    chmod +x do*

# run the script to get the .tab file with all MSA sequences of VAX004
    doall >mn2.tab
# convert .tab into .fa file
    tabToFa mn2

# grab the base alignment sequence
    echo ">hivmn2" >mn2.aln
    hgsql hivmn2 -N -e 'select seq from vax003BMsa where id="MN_D533-gp120"'  >> mn2.aln

# prepare an interium file, jjAll.mfa
    cat mn2.aln mn2.fa >jjAll.mfa
    echo = >>jjAll.mfa

# Run xmfaToMafMn2B to create a precursor file for the final .maf

    xmfaToMafMn1 jjAll.mfa j.out  org1=hivmn2
    cat j.out|sed -e 's/\./_/g'|sed -e 's/_chr/\.chr/g' >chr1.tmp

    rm jjAll.mfa j.out

    cat chr1.tmp |sed -e 's/ss_T/T/g' >chr1.maf

# copy .maf to /gbdb.

    mkdir -p  /gbdb/hivmn2/vax003BMaf
    cp chr1.maf /gbdb/hivmn2/vax003BMaf -p
    echo before load
    hgLoadMaf hivmn2 vax003BMaf

# create another copy for protein MAF.

    mkdir -p  /gbdb/hivmn2/vax003BAaMaf 
    cp -p chr1.maf /gbdb/hivmn2/vax003BAaMaf
    hgLoadMaf hivmn2 vax003BAaMaf

##############################################################
# REBUILD THE gsidClinicRecWithSeq TABLE (DONE 11/03/08, Fan)

# See details in hivVax003Vax004.txt.
##############################################################
# Create Positive Selection tracks for VAX003 subtype B

cd /hive/groups/gsid/medical/hiv/hivmn2
mkdir posSelection
cd posSelection

# BLAT /hive/groups/gsid/medical/hiv/hiva244/posSelection/BMsaAaConsensus.fa
# against hivmn2 base genome, select psl without header option
# cut and paste the result into the file BMsa.psl

hgLoadPsl -keep -table=BMsaPsl -nobin hivmn2 BMsa.psl 

# will get the following error:

#Processing BMsa.psl
#Can't start query:
#LOAD DATA CONCURRENT  INFILE
'/cluster/hive/groups/gsid/medical/hiv/hivmn2/posSelection/BMsa.psl'  INTO
TABLE BMsaPsl

#mySQL error 13: Can't get stat of
'/cluster/hive/groups/gsid/medical/hiv/hivmn2/posSelection/BMsa.psl' (Errcode:
13)

# load manually then

hgsql hivmn2
load data local infile "BMsa.psl" into table BMsaPsl;
quit

# build the positive selection tracks for model 2 and model 8.

gsidPosSelect hivmn2  BMsaPsl posSelBuild pSelectBModel2  posSelModel2.bed
hgLoadBed hivmn2 posSelModel2 posSelModel2.bed

gsidPosSelect hivmn2 BMsaPsl posSelBuild pSelectBModel8  posSelModel8.bed
hgLoadBed hivmn2 posSelModel8 posSelModel8.bed

##########################################################################
# BUILD THE POSITIVE SELECTION TRACKS FOR VAX003 SUBTYPE AE

    ssh hiv1
    mkdir -p /hive/groups/gsid/medical/hiv/posSelection/AE/hivmn2
    cd /hive/groups/gsid/medical/hiv/posSelection/AE/hivmn2

# BLAT
# /cluster/hive/groups/gsid/medical/hiv/posSelection/AE/AEMsaAaConsensus.fa
# against hivmn2 base genome, select psl without header option
# cut and paste the result into the file AEMsa.psl

hgLoadPsl -keep -table=AEMsaPsl -nobin hivmn2 AEMsa.psl 

# will get the following error:

#Processing AEMsa.psl
#Can't start query:
#LOAD DATA CONCURRENT  INFILE
'/cluster/hive/groups/gsid/medical/hiv/hivmn2/posSelection/AEMsa.psl'  INTO
TABLE AEMsaPsl

#mySQL error 13: Can't get stat of
'/cluster/hive/groups/gsid/medical/hiv/hivmn2/posSelection/AEMsa.psl'
(Errcode: 13)

# load manually then

hgsql hivmn2
load data local infile "AEMsa.psl" into table AEMsaPsl;
quit

# build positive selection tracks for model 2 and model 8.

gsidPosSelect hivmn2  AEMsaPsl posSelBuild pSelectAEModel2  posSelAEModel2.bed
hgLoadBed hivmn2 posSelAEModel2 posSelAEModel2.bed

gsidPosSelect hivmn2 AEMsaPsl posSelBuild pSelectAEModel8  posSelAEModel8.bed
hgLoadBed hivmn2 posSelAEModel8 posSelAEModel8.bed

##########################################################################
# BUILD THE POSITIVE SELECTION TRACKS FOR VAX004 (Done Fan, 3/2/09)

# Please see the corresponding section in hivVax003Vax004.txt for details.
##########################################################################

