src/hg/makeDb/doc/felCatV17e.txt 1.7
1.7 2010/04/01 16:37:14 chinhli
Remove grep 829.2 command
Index: src/hg/makeDb/doc/felCatV17e.txt
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/doc/felCatV17e.txt,v
retrieving revision 1.6
retrieving revision 1.7
diff -b -B -U 4 -r1.6 -r1.7
--- src/hg/makeDb/doc/felCatV17e.txt 31 Mar 2010 21:32:36 -0000 1.6
+++ src/hg/makeDb/doc/felCatV17e.txt 1 Apr 2010 16:37:14 -0000 1.7
@@ -72,19 +72,19 @@
# For unplaced scalfolds, named them as chrUn_xxxxxxxx
# where xxxxxx is the original access id as: chrUn_ACBE01000381.1
# and put it into chrUn.* files
# copy all the comment lines (start with #)
+S=Primary_Assembly/unplaced_scaffolds
zcat ${S}/AGP/unplaced.scaf.agp.gz | grep "^#" > ucscChr/chrUn.agp
# append the gap records
zcat ${S}/AGP/unplaced.scaf.agp.gz | grep -v "^#" \
| sed -e "s/^/chrUn_/" >> ucscChr/chrUn.agp
gzip ucscChr/chrUn.agp &
-S=Primary_Assembly/unplaced_scaffolds
zcat ${S}/FASTA/unplaced.scaf.fa.gz \
| sed -e "s#^>.*|gb|#>chrUn_#; s#|.*##" | gzip > ucscChr/chrUn.fa.gz
# about 104034 sequences in the unplaced
-zcat chrUn.fa.gz | grep "^>" | wc
+zcat ucscChr/chrUn.fa.gz | grep "^>" | wc
# 104034 104034 2275648
# Check them with faSize
faSize Primary_Assembly/unplaced_scaffolds/FASTA/unplaced.scaf.fa.gz
@@ -96,8 +96,21 @@
##########################################################################
# Initial genome build (Done - 2010-03-03 -Chin)
+
+# decide the right order by
+ hgsql -e "select orderKey,name from dbDb order by orderKey;" \
+ hgcentraltest | less
+ # 168 cavPor3
+ # 169 cavPor2
+ # 189 oryCun2
+ # 190 oryCun1
+ # 191 ochPri2
+ # 200 eriEur1
+ # 217 felCat3
+
+
cd /hive/data/genomes/felCatV17e
cat << '_EOF_' > felCatV17e.config.ra
# Config parameters for makeGenomeDb.pl: