#NEXUS

BEGIN TAXA;
	TITLE Taxa;
	DIMENSIONS NTAX=10;
	TAXLABELS
		Anolis_carolinensis Anser_erythropus Gallus_gallus Pitta_guajana Megalaima_virens Urocolius_indicus Phalocrocorax_carbo Eudromia_elegans Dromaius_novahollandiae Struthio_camelus 
	;

END;

BEGIN TREES;
	Title 'deep-bird-each-part';
	LINK Taxa = Taxa;
	TRANSLATE
		1 Anolis_carolinensis,
		2 Anser_erythropus,
		3 Gallus_gallus,
		4 Pitta_guajana,
		5 Megalaima_virens,
		6 Urocolius_indicus,
		7 Phalocrocorax_carbo,
		8 Eudromia_elegans,
		9 Dromaius_novahollandiae,
		10 Struthio_camelus;
	TREE con_50_majrule = (((((4:0.018805,(5:0.021407,6:0.014289)1.00:6.82E-4)1.00:6.21E-4,7:0.007659)1.00:0.002343,(2:0.007618,3:0.012711)1.00:0.002677)1.00:0.002551,((8:0.015589,9:0.006108)1.00:0.001447,10:0.008537)1.00:0.003334):0.0523295,1:0.0523295);

END;

begin mrbayes;
	charset HKYI = 1 - 81008;
	charset GTRG = 81009 - 147480;
	charset JC69I = 147481 - 147858;
	charset HKYG = 147859 - 242253;
	charset GTR = 242254 - 254430;
	charset GTRI = 254431 - 298276;
	charset GTRIG = 298277 - 303680;
	charset SYM = 303681 - 303988;
	charset HKY = 303989 - 338215;
	charset SYMG = 338216 - 338831;
	charset HKYIG = 338832 - 344773;
	charset F81 = 344774 - 346521;
	charset F81I = 346522 - 347208;
	charset F81G = 347209 - 348085;
	charset K2PI = 348086 - 349068;
	charset K2PG = 349069 - 350617;
	charset K2P = 350618 - 352120;
	partition partial = 17 : HKYI , GTRG , JC69I , HKYG , GTR , GTRI , GTRIG , SYM , HKY , SYMG , HKYIG , F81 , F81I , F81G , K2PI , K2PG , K2P;
	set partition = partial;
	lset applyto = ( 1 ) nst = 2 rates = propinv;
	lset applyto = ( 2 ) nst = 6 rates = gamma;
	lset applyto = ( 3 ) nst = 1 rates = propinv;
	prset applyto = ( 3 ) statefreqpr = fixed ( equal );
	lset applyto = ( 4 ) nst = 2 rates = gamma;
	lset applyto = ( 5 ) nst = 6;
	lset applyto = ( 6 ) nst = 6 rates = propinv;
	lset applyto = ( 7 ) nst = 6 rates = invgamma;
	lset applyto = ( 8 ) nst = 6;
	prset applyto = ( 8 ) statefreqpr = fixed ( equal );
	lset applyto = ( 9 ) nst = 2;
	lset applyto = ( 10 ) nst = 6 rates = gamma;
	prset applyto = ( 10 ) statefreqpr = fixed ( equal );
	lset applyto = ( 11 ) nst = 2 rates = invgamma;
	lset applyto = ( 12 ) nst = 1;
	lset applyto = ( 13 ) nst = 1 rates = propinv;
	lset applyto = ( 14 ) nst = 1 rates = gamma;
	lset applyto = ( 15 ) nst = 2 rates = propinv;
	prset applyto = ( 15 ) statefreqpr = fixed ( equal );
	lset applyto = ( 16 ) nst = 2 rates = gamma;
	prset applyto = ( 16 ) statefreqpr = fixed ( equal );
	lset applyto = ( 17 ) nst = 2;
	prset applyto = ( 17 ) statefreqpr = fixed ( equal );
	unlink statefreq = ( all ) revmat = ( all ) shape = ( all ) pinvar = ( all ) tratio = ( all );

END;