Min length of gene families about pirate HOT 8 CLOSED

Hi Tim,

There isn't a parameter that can be accessed from the main command currently. If you want a quick hack the you can modify the line below in run_PIRATE.pl to include "-t 75"

	# extract sequences.
	`cat $pirate_dir/genome_list.txt | parallel -k -j $threads \"perl $script_path/extract_feature_sequences.pl -s {} -d $pirate_dir -o $pirate_dir/genome_multifastas/{}.fasta $e_args >> $pirate_dir/gff_parser_log.txt 2>> $pirate_dir/gff_parser_log.txt\"`;

e.g.

	# extract sequences.
	`cat $pirate_dir/genome_list.txt | parallel -k -j $threads \"perl $script_path/extract_feature_sequences.pl -t 75 -s {} -d $pirate_dir -o $pirate_dir/genome_multifastas/{}.fasta $e_args >> $pirate_dir/gff_parser_log.txt 2>> $pirate_dir/gff_parser_log.txt\"`;

I am happy to add a branch with this option to the main script and include it in the next release.

All the best,
Sion

from pirate.

Hi Sion - thanks for the super quick reply. Will use this hack and also would be interested in having a branch available in the next release. Tim

from pirate.

When I get a spare minute I will set one up for you :)

from pirate.

I don't see the nohup warnings you were getting when I created a branch with the same option. Are you still seeing them and are the outputs okay?

S

from pirate.

I found an error message with link_clysters.pl but it looks like it maybe an issue of rerunning PIRATE in the same directory with different settings. When I sent outputs to a fresh directory it didn't occur. I will push to master shortly. Thank you for your help and testing.

S

from pirate.

Hi @SionBayliss!

I noticed the PR for this feature was merged then reverted (#55). Are there still plans to implement this in the future?

Thanks!

from pirate.

I have merged this with another commit that I wanted to push together. The revert was due to reading about the nohup errors after merging and then slightly panicking.

I have pushed them both to master now. The master branch has all of the changes that I previously reverted.

from pirate.

I should have waited! Thanks Sion!

from pirate.