a44421a79fb36cc2036fe116b97ea3bc9590cd0c braney Fri Dec 2 09:34:39 2011 -0800 removed rcsid (#295) diff --git src/utils/faFrag/faFrag.c src/utils/faFrag/faFrag.c index 8b6f9c5..ea12d4f 100644 --- src/utils/faFrag/faFrag.c +++ src/utils/faFrag/faFrag.c @@ -1,56 +1,55 @@ /* faFrag - Extract a piece of DNA from a .fa file.. */ #include "common.h" #include "dnaseq.h" #include "fa.h" #include "options.h" -static char const rcsid[] = "$Id: faFrag.c,v 1.5 2006/03/18 01:54:46 angie Exp $"; void usage() /* Explain usage and exit. */ { errAbort( "faFrag - Extract a piece of DNA from a .fa file.\n" "usage:\n" " faFrag in.fa start end out.fa\n" "options:\n" " -mixed - preserve mixed-case in FASTA file\n"); } void faFrag(char *inName, int start, int end, char *outName, boolean mixed) /* faFrag - Extract a piece of DNA from a .fa file.. */ { struct dnaSeq *seq; char name[512]; if (start >= end) usage(); if (mixed) seq = faReadAllMixed(inName); else seq = faReadAllDna(inName); if (seq->next != NULL) warn("More than one sequence in %s, just using first\n", inName); if (end > seq->size) { warn("%s only has %d bases, truncating", seq->name, seq->size); end = seq->size; if (start >= end) errAbort("Sorry, no sequence left after truncating"); } sprintf(name, "%s:%d-%d", seq->name, start, end); faWrite(outName, name, seq->dna + start, end-start); printf("Wrote %d bases to %s\n", end-start, outName); } int main(int argc, char *argv[]) /* Process command line. */ { optionHash(&argc, argv); if (argc != 5) usage(); if (!isdigit(argv[2][0]) || !isdigit(argv[3][0])) usage(); faFrag(argv[1], atoi(argv[2]), atoi(argv[3]), argv[4], optionExists("mixed")); return 0; }