Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

change charset to charpartition for genomes #45

Open
roblanf opened this issue Aug 30, 2023 · 0 comments
Open

change charset to charpartition for genomes #45

roblanf opened this issue Aug 30, 2023 · 0 comments

Comments

@roblanf
Copy link
Owner

roblanf commented Aug 30, 2023

Right now when I run one of the alignments in IQ-TREE it loads all the charsets:

iqtree2 -s alignment.nex -p alignment.nex --redo

Loading 5 partitions...
Subset  Type    Seqs    Sites   Infor   Invar   Model   Name
1               145     531     54      425     1       COI_1stpos
2               145     531     3       502     2       COI_2ndpos
3               145     530     303     109     3       COI_3rdpos
4               145     1445    283     924     4       16S
5               145     3037    643     1960    1       mitochondrial_genome
Degree of missing data: 0.000
Info: multi-threading strategy over partitions

So instead of this:

begin SETS;

	[loci]
	CHARSET	COI_1stpos = 1-1592\3;
	CHARSET	COI_2ndpos = 2-1592\3;
	CHARSET	COI_3rdpos = 3-1592\3;
	CHARSET	16S = 1593-3037;

    CHARPARTITION loci = 1:COI_1stpos, 2:COI_2ndpos, 3:COI_3rdpos, 4:16S;

	[genomes]
	CHARSET	mitochondrial_genome = 1-3037;

    CHARPARTITION genomes = 1:mitochondrial_genome;

end;

we should have this:

begin SETS;

	[loci]
	CHARSET	COI_1stpos = 1-1592\3;
	CHARSET	COI_2ndpos = 2-1592\3;
	CHARSET	COI_3rdpos = 3-1592\3;
	CHARSET	16S = 1593-3037;

    CHARPARTITION loci = 1:COI_1stpos, 2:COI_2ndpos, 3:COI_3rdpos, 4:16S;

	[genomes]
	CHARPARTITION	mitochondrial_genome = 1-3037;

    CHARPARTITION genomes = 1:mitochondrial_genome;

end;
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant