src/hg/makeDb/doc/h1n1.txt 1.6
1.6 2009/05/18 20:43:42 fanhsu
Added dnaSeq table.
Index: src/hg/makeDb/doc/h1n1.txt
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/doc/h1n1.txt,v
retrieving revision 1.5
retrieving revision 1.6
diff -b -B -U 4 -r1.5 -r1.6
--- src/hg/makeDb/doc/h1n1.txt 15 May 2009 22:46:08 -0000 1.5
+++ src/hg/makeDb/doc/h1n1.txt 18 May 2009 20:43:42 -0000 1.6
@@ -224,8 +224,25 @@
hgsql h1n1 -e 'delete from h1n1SeqXref'
hgsql h1n1 -e 'load data local infile "h1n1SeqXref.tab" into table h1n1SeqXref'
+# create dna sequences table
+
+ cp -p ~/h1n1/download/0514/epiflu_0514_dna_sequence.fa j1
+
+ ~/src/utils/faToTab/faToTab.pl /dev/null j1 >j2.tab
+ cut -f 1 j2.tab >jj1
+ cut -f 2 j2.tab >jseq
+
+ cat jj1 |sed -e 's/ | /\t/g'|sed -e 's/>//' >jj
+ cut -f 1 jj > jDnaId
+
+ paste jDnaId jseq >epiflu_0514_dna_sequence.tab
+
+ hgLoadSqlTab h1n1 dnaSeq ~/src/hg/lib/dnaSeq.sql epiflu_0514_dna_sequence.tab
+
+ rm j1 jj jj1 jseq jDnaId j2.tab
+
# create protein sequences file with appropriate IDs
cp -p ~/h1n1/download/0514/epiflu_0514_protein_sequence.fa j1