Problem with creating core alignments about pirate HOT 3 CLOSED

It is likely that you don't have any clusters that the alignment script considers core i.e. >95% of isolates and dosage/copy number <1.25. You can check that this is the case in the tabular output in the log or summary file. If this is the case and you want to modify the default thresholds to use a more relaxed core definition then you can run the script(s) manually.

create_pangenome_alignment.pl -i /path/to/PIRATE.gene_families.tab -f /path/to/feature_sequences/ -o /path/to/output_fasta.fasta -g /path/to/output_gff.gff -t percentage_threshold (integer) -d dosage_threshold(floating_point)

If dosage/copy number is an issue you may have a problematic dataset. If this is the case let me know and I can suggest some advanced settings that should improve cluster resolution.

from pirate.

Is there somewhere where the options for core threshold and dosage threshold are documented?

I'm running PIRATE v.1.0.3 and the "-t" option seems to set threads, not core threshold, and I can't find the "-d" option in any of the help menus.

Thanks!
Conrad

from pirate.

Hi Conrad,

You will need to run the script create_pangenome_alignment.pl (in /PIRATE/scripts/ directory) separately as detailed above. There is currently no options for changing these defaults from the main PIRATE script. I will consider adding them in a future update.

All the best,
Sion

from pirate.