src/hg/makeDb/doc/mm9.txt 1.98

1.98 2009/06/11 22:28:31 hartera
Fixed an off by 1 error in the BED file loaded for miRNA track.
Index: src/hg/makeDb/doc/mm9.txt
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/doc/mm9.txt,v
retrieving revision 1.97
retrieving revision 1.98
diff -b -B -U 4 -r1.97 -r1.98
--- src/hg/makeDb/doc/mm9.txt	11 Jun 2009 00:11:49 -0000	1.97
+++ src/hg/makeDb/doc/mm9.txt	11 Jun 2009 22:28:31 -0000	1.98
@@ -9363,9 +9363,9 @@
     cat fb.mm9.chainHg19Link.txt 
     #	1013880568 bases of 2620346127 (38.693%) in intersection
 
 #############################################################################
-#  RE-BUILD miRNA TRACK (DONE, 2009-06-09-2009-06-10, hartera)
+#  RE-BUILD miRNA TRACK (DONE, 2009-06-09-2009-06-11, hartera)
     # The miRNA track from miRBase is out of date so update the track. 
     mkdir -p /hive/data/genomes/mm9/bed/miRNA-2009-06-09
     cd /hive/data/genomes/mm9/bed/miRNA-2009-06-09
     # Download GFF file of latest miRNA annotations from miRBase at the
@@ -9393,11 +9393,13 @@
     # chr1    .       miRNA   20669091        20669163        .       +
     # .       mmu-mir-206
     # use score 906 for + strand and 480 for - strand. This will show 
     # up black on the track for + strand and grey for - strand.
+    # Re-do below and re-load track as appears off by 1 compared to 
+    # Ensembl track and other miRNA resources (2009-06-11)
     awk 'BEGIN {FS="\t"} {OFS="\t"} \
-        {if ($0 !~ /#/ && $7 == "+") print $1, $4, $5, $9, 960, $7; \
-       else if ($0 !~ /#/ && $7 == "-") print $1, $4, $5, $9, 480, $7;}' \
+        {if ($0 !~ /#/ && $7 == "+") print $1, $4-1, $5, $9, 960, $7; \
+       else if ($0 !~ /#/ && $7 == "-") print $1, $4-1, $5, $9, 480, $7;}' \
         mmMirBaseFormatIdOnly.gff > mmMirBaseFormatIdOnly.bed
     # Remove previous table
     hgsql -e 'drop table miRNA' mm9
     hgLoadBed mm9 miRNA mmMirBaseFormatIdOnly.bed