Commit 0f1e495e authored by Mathieu Giraud's avatar Mathieu Giraud

vidjil.cpp: the experimental '-g' options takes the path of germlines as argument

This option is now undocumented, its use will change in a future release.
parent 479944e6
......@@ -107,10 +107,10 @@ void MultiGermline::insert(Germline *germline)
germlines.push_back(germline);
}
void MultiGermline::load_default_set()
void MultiGermline::load_default_set(string path)
{
germlines.push_back(new Germline("TRG", 'G', "germline/TRGV.fa", "", "germline/TRGJ.fa", "#####-#####", -10, 20));
germlines.push_back(new Germline("IGH", 'H', "germline/IGHV.fa", "germline/IGHD.fa", "germline/IGHJ.fa", "######-######", 0, 80));
germlines.push_back(new Germline("TRG", 'G', path + "/TRGV.fa", "", path + "/TRGJ.fa", "#####-#####", -10, 20));
germlines.push_back(new Germline("IGH", 'H', path + "/IGHV.fa", path + "/IGHD.fa", path + "/IGHJ.fa", "######-######", 0, 80));
}
......
......@@ -67,7 +67,7 @@ class MultiGermline {
~MultiGermline();
void insert(Germline *germline);
void load_default_set();
void load_default_set(string path);
void out_stats(ostream &out);
};
......
......@@ -147,7 +147,7 @@ void usage(char *progname)
<< " -D <file> D germline multi-fasta file (automatically implies -d)" << endl
<< " -J <file> J germline multi-fasta file" << endl
<< " -G <prefix> prefix for V (D) and J repertoires (shortcut for -V <prefix>V.fa -D <prefix>D.fa -J <prefix>J.fa)" << endl
<< " -g multiple germlines (experimental)" << endl
<< " -g <path> multiple germlines (experimental)" << endl
<< endl
<< "Window prediction" << endl
......@@ -279,6 +279,7 @@ int main (int argc, char **argv)
bool output_segmented = false;
bool output_unsegmented = false;
bool multi_germline = false;
string multi_germline_file = "";
string forced_edges = "" ;
......@@ -294,7 +295,7 @@ int main (int argc, char **argv)
//$$ options: getopt
while ((c = getopt(argc, argv, "AhagG:V:D:J:k:r:vw:e:C:t:l:dc:m:M:N:s:p:Sn:o:L%:Z:y:z:uU")) != EOF)
while ((c = getopt(argc, argv, "Ahag:G:V:D:J:k:r:vw:e:C:t:l:dc:m:M:N:s:p:Sn:o:L%:Z:y:z:uU")) != EOF)
switch (c)
{
......@@ -355,7 +356,8 @@ int main (int argc, char **argv)
break;
case 'g':
multi_germline = true ;
multi_germline = true;
multi_germline_file = string(optarg);
break;
case 'G':
......@@ -785,7 +787,7 @@ int main (int argc, char **argv)
if (multi_germline)
{
multigermline->load_default_set();
multigermline->load_default_set(multi_germline_file);
}
else
{
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment