d486edad7214c55c607e4c8302ec1b0e10a1daf8 hiram Mon Mar 2 14:13:48 2020 -0800 correct one line header in usage message no redmine diff --git src/utils/cpgIslandExt/cpg_lh.c src/utils/cpgIslandExt/cpg_lh.c index 93e15fc..5baa608 100644 --- src/utils/cpgIslandExt/cpg_lh.c +++ src/utils/cpgIslandExt/cpg_lh.c @@ -48,32 +48,32 @@ had CPGSCORE 27, but allowed score to reach 0 once without reporting */ #define MALLOCBLOCK 200000 #define MAXNAMELEN 512 /*"length 0, title > MAXNAMELEN or starts with bad character." if too small */ /*------------------------------------------------------*/ int readSequence (FILE *fil, int *conv, char **seq, char **id, char **desc, int *length) ; void findspans ( int start, int end, char *seq, char *seqname ) ; void getstats ( int start, int end, char *seq, char *seqname, int *ncpg, int *ngpc, int *ng, int *nc ) ; void usage (void) { - fprintf(stderr, "cpglh - calculate CpG Island data for cpgIslandExt tracks\n"); - fprintf(stderr, "usage:\n cpglh \n") ; + fprintf(stderr, "cpg_lh - calculate CpG Island data for cpgIslandExt tracks\n"); + fprintf(stderr, "usage:\n cpg_lh \n") ; fprintf(stderr, "where is fasta sequence, must be more than\n"); fprintf(stderr, " 200 bases of legitimate sequence, not all N's\n"); fprintf(stderr, "\nTo process the output into the UCSC bed file format:\n\n" "cpglh fastaInput.fa \\\n" " | awk '{$2 = $2 - 1; width = $3 - $2;\n" " printf(\"%%s\\t%%d\\t%%s\\t%%s %%s\\t%%s\\t%%s\\t%%0.0f\\t%%0.1f\\t%%s\\t%%s\\n\",\n" " $1, $2, $3, $5, $6, width, $6, width*$7*0.01, 100.0*2*$6/width, $7, $9);}' \\\n" \ " | sort -k1,1 -k2,2n > output.bed\n"); fprintf(stderr, "\nThe original cpg.c was written by Gos Miklem from the Sanger Center.\n" "LaDeana Hillier added some modifications --> cpg_lh.c, and UCSC hass\n" "added some further modifications to cpg_lh.c, so that its expected\n" "number of CpGs in an island is calculated as described in\n" " Gardiner-Garden, M. and M. Frommer, 1987\n" " CpG islands in vertebrate genomes. J. Mol. Biol. 196:261-282\n"