| 277 | } |
| 278 | |
| 279 | void wrong_input() |
| 280 | { |
| 281 | printf("Usage: dialign-t [OPTIONS] <conf-directory> <fasta-file> [<fasta-out-file>]\n"); |
| 282 | printf("\n -d\tDebug-Mode [DEFAULT 0]\n"); |
| 283 | printf(" \t\t 0 no debug statements\n"); |
| 284 | printf(" \t\t 1 debugs the current phase of the processing\n"); |
| 285 | printf(" \t\t 2 very loquacious debugging\n"); |
| 286 | printf(" \t\t 5 hardcore debugging\n"); |
| 287 | printf(" -s\tmaximum amount of input sequences [DEFAULT 5000]\n"); |
| 288 | printf(" -a\tmaximum number of characters per line in a FASTA file [DEFAULT 100]\n"); |
| 289 | printf(" -c\tmaximum amount of characters per line when printing a sequence\n \t[DEFAULT 80]\n"); |
| 290 | printf(" -l\tsensitivity mode, the higher the level the less likely\n"); |
| 291 | printf(" \tspurious random fragments are aligned in local alignments \n \t[DEFAULT 0]\n"); |
| 292 | printf(" \t\t 0 switched off \n"); |
| 293 | printf(" \t\t 1 level-1, reduced sensitivity\n"); |
| 294 | printf(" \t\t 2 level-2, strongly reduced sensitivity\n"); |
| 295 | printf(" -m\tscore matrix file name (in the configuration directory)\n \t\t[DEFAULT PROTEIN: BLOSUM.scr]\n \t\t[DEFAULT DNA: dna_matrix.scr]\n"); |
| 296 | printf(" -w\tdefines the minimum weight when the weight formula is changed\n \tto 1-pow(1-prob, factor) [DEFAULT 0.000000065]\n"); |
| 297 | printf(" -p\tprobability distribution file name (in the configuration\n \tdirectory) \n \t\t[DEFAULT PROTEIN: BLOSUM.diag_prob_t10]\n\t\t[DEFAULT DNA: dna_diag_prob_100_exp_550000]\n"); |
| 298 | printf(" -v\tadd to each score (to prevent negative values) [DEFAULT 0]\n"); |
| 299 | printf(" -t\t\"even\" threshold for low score for sequences alignment \n \t\t[DEFAULT PROTEIN: 4]\n\t\t[DEFAULT DNA: 0]\n"); |
| 300 | printf(" -n\tmaximum number of consecutive positions for window containing\n \tlow scoring positions \n \t\t[DEFAULT PROTEIN: 4]\n\t\t[DEFAULT DNA: 4]\n"); |
| 301 | printf(" -g\tglobal minimum fragment length for stop criterion \n \t\t[DEFAULT PROTEIN: 40] \n\t\t[DEFAULT DNA: 40]\n"); |
| 302 | printf(" -m\tminimal allowed average score in frag window containing low \n \tscoring positions \n \t\t[DEFAULT PROTEIN: 4.0]\n\t\t[DEFAULT DNA: 0.25]\n"); |
| 303 | printf(" -o\twhether overlap weights are calculated or not [DEFAULT 0]\n"); |
| 304 | printf(" -f\tminimum fragment length [DEFAULT 1]\n"); |
| 305 | printf(" -r\tthreshold weight to consider the fragment at all [DEFAULT 0.0]\n"); |
| 306 | printf(" -u\t[DEFAULT 0]\n"); |
| 307 | printf(" \t\t1: only use a sqrt(amount_of_seqs) stripe of neighbour\n \t\t sequences to calculate pairwise alignments (increase performance)\n"); |
| 308 | printf(" \t\t0: all pairwise alignments will be calculated\n"); |
| 309 | printf(" -A\toptional anchor file [DEFAULT none]\n"); |
| 310 | |
| 311 | printf(" -D\tinput is DNA-sequence\n"); |
| 312 | printf(" -T\ttranslate DNA into aminoacids from begin to end (length will be cut to mod 3 = 0)\n\tWARNING: Do not use -D with this option \n\t(Default values for PROTEIN input will be loaded)\n"); |
| 313 | printf(" -L\tcompare only longest Open Reading Frame\n\tWARNING: Do not use -D with this option \n\t(Default values for PROTEIN input will be loaded)\n"); |
| 314 | printf(" -O\ttranslate DNA to aminoacids, reading frame for each sequence calculated due to its longest ORF\n\tWARNING: Do not use -D with this option \n\t(Default values for PROTEIN input will be loaded)\n"); |
| 315 | printf(" -P\toutput in aminoacids, no retranslation of DNA sequences\n\t[DEFAULT: input = output]\n"); |
| 316 | printf(" -F\tfast mode (implies -l0, since it already significantly reduces sensitivity)\n"); |
| 317 | printf(" -C\tgenerate probability table saved in <config_dir>/prob_table and exit\n"); |
| 318 | printf(" -H -h\tprint this message\n\n"); |
| 319 | exit(1); |
| 320 | } |
| 321 | |
| 322 | void parameters(int argc, char** argv) |
| 323 | { |