def generate_argparser():
"""Generate argparse parser
"""
parser = MvfArgumentParser()
parser.addarg_mvf()
parser.addarg_out()
parser.addarg_sample_indices(nmin=2)
parser.addarg_sample_labels(nmin=2)
parser.addarg_windowsize()
parser.addarg_mincoverage()
parser.add_argument("--data-type",
"--datatype",
choices=("dna", "prot"),
help=("Data type to compare."
"(This option is only needed for codon "
" MVF files, others will default.)"))
parser.add_argument("--ambig",
choices=("random2", "random3"),
help=("By default, ambiguous nucleotides are "
"excluded. This option will include "
"sets of ambiguous characters by "
"randomly choosing one of the options "
"for: RYMKWS ('random2') or "
"RYMKWS+BDHV ('random3')"))
parser.add_argument(
"--emit-counts",
action="store_true",
help=("output additional file that presents "
"the raw counts of pairwise patterns for "
"each sample pair tested for each window"))
return parser
def generate_argparser():
"""Generate argparse parser
"""
parser = MvfArgumentParser()
parser.addarg_mvf()
parser.addarg_out()
parser.addarg_sample_indices()
parser.addarg_sample_labels()
parser.addarg_windowsize()
parser.addarg_mincoverage()
return parser
def generate_argparser():
"""Generate argparse parser
"""
parser = MvfArgumentParser()
parser.addarg_mvf()
parser.addarg_out()
parser.addarg_sample_indices()
parser.addarg_sample_labels()
parser.addarg_windowsize()
parser.addarg_mincoverage()
parser.add_argument("--output-lists", action="store_true")
return parser
def generate_argparser():
"""Generate argparse parser
"""
parser = MvfArgumentParser()
parser.addarg_mvf()
parser.addarg_out()
parser.addarg_contig_ids()
parser.addarg_contig_labels()
parser.addarg_sample_indices()
parser.addarg_sample_labels()
parser.addarg_mincoverage()
parser.addarg_windowsize()
parser.add_argument(
"--base-match",
"--basematch",
help="String of bases to match (i.e. numerator).")
parser.add_argument(
"--base-total",
"--basetotal",
help="String of bases for total (i.e. denominator).")
return parser
def generate_argparser():
"""Generate argparse parser
"""
pallette = Pallette()
parser = MvfArgumentParser()
parser.addarg_mvf()
parser.add_argument("--out-prefix",
"--outprefix",
help="Output prefix (not required).")
parser.addarg_sample_indices(nmin=3)
parser.addarg_sample_labels(nmin=3)
parser.addarg_outgroup_indices(nmin=1)
parser.addarg_outgroup_labels(nmin=1)
parser.addarg_windowsize()
parser.add_argument(
"--contig-labels",
"--contiglabels",
nargs=1,
help=("Enter the ids of one or more contigs in the "
"order they will appear in the chromoplot (as "
"comma-separated list)"
"(defaults to all ids in order present in MVF)"))
parser.add_argument(
"--contig-ids",
"--contigids",
"--contigs",
nargs=1,
help=("Enter the labels of one or more contigs in the "
"order they will appear in the chromoplot (as "
"comma-separated list)"
"(defaults to all ids in order present in MVF)"))
parser.add_argument(
"--majority",
action="store_true",
help=("Plot only 100% shading in the majority track "
" rather than shaded proportions in all tracks."))
parser.add_argument(
"--info-track",
"--infotrack",
action="store_true",
help=("Include an additional coverage information "
"track that will show empty, uninformative, "
"and informative loci. (Useful for "
"ranscriptomes/RAD or other reduced sampling."))
parser.add_argument("--empty-mask",
"--emptymask",
choices=pallette.colornames,
default="none",
help="Mask empty regions with this color.")
parser.add_argument(
"--yscale",
default=20,
type=int,
help=("Height (in number of pixels) for each track"))
parser.add_argument(
"--xscale",
default=1,
type=int,
help="Width (in number of pixels) for each window")
parser.add_argument("--colors",
nargs=3,
choices=pallette.colornames,
help="three colors to use for chromoplot")
parser.add_argument(
"--plot-type",
"--plottype",
choices=["graph", "image"],
default="image",
help=("PNG image (default) or graph via matplotlib "
"(experimental)"))
return parser
def generate_argparser():
"""Generate argparse parser
"""
parser = MvfArgumentParser()
parser.addarg_mvf()
parser.addarg_out()
parser.addarg_sample_indices()
parser.addarg_sample_labels()
parser.addarg_contig_ids()
parser.addarg_contig_labels()
parser.addarg_windowsize()
parser.add_argument("--raxml-outgroups",
"--raxmloutgroups",
help=("Comma-separated list of outgroup "
"taxon labels to use in RAxML."))
parser.add_argument("--root-with",
"--rootwith",
help=("Comma-separated list of taxon labels "
"to root trees with after RAxML"))
parser.add_argument("--output-contig-labels",
"--outputcontiglabels",
action="store_true",
help=("Output will use contig labels "
"instead of id numbers."))
parser.add_argument("--output-empty",
"--outputempty",
action="store_true",
help=("Include entries of windows "
"with no data in output."))
parser.add_argument(
"--choose-allele",
"--chooseallele",
"--hapmode",
default="none",
dest="choose_allele",
choices=[
"none", "randomone", "randomboth", "major", "minor",
"majorminor"
],
help=("Chooses how heterozygous alleles are "
"handled. (none=no splitting (default); "
"randomone=pick one allele randomly "
"(recommended); randomboth=pick two alleles "
"randomly, but keep both; major=pick the "
"more common allele; minor=pick the less "
"common allele; majorminor= pick the major in "
"'a' and minor in 'b'"))
parser.add_argument("--min-sites",
"--minsites",
type=int,
default=100,
help="minimum number of sites ")
parser.add_argument(
"--min-seq-coverage",
"--minseqcoverage",
type=float,
default=0.1,
help=("proportion of total alignment a sequence"
"must cover to be retianed [0.1]"))
parser.add_argument("--min-depth",
"--mindepth",
type=int,
default=4,
help=("minimum number of alleles per site"))
parser.add_argument(
"--bootstrap",
type=int,
help=("turn on rapid bootstrapping for RAxML and "
"perform specified number of replicates"))
parser.add_argument("--raxml-model",
"--raxmlmodel",
default="GTRGAMMA",
help=("choose RAxML model"))
parser.add_argument("--raxml-path",
"--raxmlpath",
default="raxml",
help="RAxML path for manual specification.")
parser.add_argument(
"--raxml-opts",
"--raxmlopts",
default="",
help=("specify additional RAxML arguments as a "
"double-quotes encased string"))
parser.add_argument(
"--duplicate-seq",
"--duplicateseq",
default="dontuse",
choices=["dontuse", "keep", "remove"],
help=("dontuse=remove duplicate sequences prior to "
"RAxML tree inference, then add them to the "
"tree manually as zero-branch-length sister "
"taxa; keep=keep in for RAxML tree inference "
"(may cause errors for RAxML); "
"remove=remove entirely from alignment"))
parser.add_argument("--temp-dir",
"--tempdir",
default='./raxmltemp',
type=os.path.abspath,
help=("Temporary directory path"))
parser.add_argument("--temp-prefix",
"--tempprefix",
default="mvftree",
help=("Temporary file prefix"))
return parser
def generate_argparser():
"""Generate argparse parser
"""
parser = MvfArgumentParser()
parser.addarg_mvf()
parser.addarg_gff()
parser.addarg_out()
parser.addarg_samples()
parser.addarg_windowsize()
parser.addarg_mincoverage()
parser.add_argument(
"--allele-groups",
"--allelegroups",
nargs='*',
help=("GROUP1:LABEL,LABEL GROUP2:LABEL,LABEL "))
parser.add_argument("--species-groups",
"--speciesgroups",
nargs='*')
parser.add_argument(
"--chi-test",
"--chitest",
help=("Input two number values for expected "
"Nonsynonymous and Synonymous expected values."))
parser.add_argument(
"--target",
nargs="*",
help=("Specify the taxa labels that define the "
"target lineage-specific branch to be tested."))
parser.add_argument("--num-target-species",
"--targetspec",
type=int,
default=1,
help=("Specify the minimum number of "
"taxa in the target set "
"that are required to conduct analysis"))
parser.add_argument("--output-align",
"--outputalign",
help=("Output alignment to this file path in "
"phylip format."))
parser.add_argument(
"--outgroup",
help=("Specify sample name with which to root trees."))
parser.add_argument(
"--use-labels",
"--uselabels",
action="store_true",
help="Use contig labels instead of IDs in output.")
parser.add_argument("--codeml-path",
"--codemlpath",
default="codeml",
type=os.path.abspath,
help="Full path for PAML codeml executable.")
parser.add_argument("--raxml-path",
"--raxmlpath",
type=os.path.abspath,
default="raxml",
help="Full path to RAxML program executable.")
parser.add_argument("--start-contig",
"--startcontig",
type=int,
default=0,
help="Numerical ID for the starting contig.")
parser.add_argument("--end-contig",
"--endcontig",
type=int,
default=100000000,
help="Numerical id for the ending contig.")
parser.add_argument(
"--paml-tmp",
"--pamltmp",
default="pamltmp",
type=os.path.abspath,
help="path for temporary folder for PAML output files")
parser.add_argument(
"--all-sample-trees",
"--allsampletrees",
action="store_true",
help=("Makes trees from all samples instead of "
"only the most complete sequence from "
"each species"))
parser.add_argument(
"--branch-lrt",
"--branchlrt",
type=os.path.abspath,
help=("Specify the output file for and turn on the "
"RAxML-PAML format LRT test scan for "
"selection on the target branch in addition "
"to the basic patterns scan"))
parser.add_argument("--verbose",
action="store_true",
help="additional screen output")
return parser
