#BREAKPOINTS DETECTED, ALIGNMENT ACCORDINGLY SPLITTED

################################################
#CALCULATIONS ON FRAGMENT 3

#BUILD PHYLOGENY fragment_3
raxmlHPC-PTHREADS-SSE3 -T 45 -f a -# 100 -x 1234 -p 1234 -s $PROJECT_DIR/fragments/fragment_3/aln/bats_mx1.nt_ali.nogaps.fasta -n nt -m GTRGAMMA -N 1000 -w $PROJECT_DIR/fragments/fragment_3/phylo/

nw_display -v 50 -i 'font-size:11' -l 'font-size:16;font-family:helvetica;font-style:italic' -d 'stroke:black;fill:none;stroke-width:2;stroke-linecap:round' -Il -w 650 -b 'opacity:0' -s $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.nt.corrected > $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.nt.scale.svg

nw_reroot $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.nt.corrected  > $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.nt.corrected.rooted

raxmlHPC-PTHREADS-SSE3 -T 45 -f a -# 100 -x 1234 -p 1234 -s $PROJECT_DIR/fragments/fragment_3/aln/bats_mx1.aa_ali.nogaps.fasta -n aa -m PROTGAMMAWAG -N 1000 -w $PROJECT_DIR/fragments/fragment_3/phylo/

nw_display -v 50 -i 'font-size:11' -l 'font-size:16;font-family:helvetica;font-style:italic' -d 'stroke:black;fill:none;stroke-width:2;stroke-linecap:round' -Il -w 650 -b 'opacity:0' -s $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.aa.corrected > $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.aa.scale.svg

nw_reroot $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.aa.corrected  > $PROJECT_DIR/fragments/fragment_3/phylo/RAxML_bipartitionsBranchLabels.aa.corrected.rooted

#RUN CODEML ON fragment_3
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F61/M0/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F61/M1a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F61/M2a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F61/M7/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F61/M8/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F61/M8a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F1X4/M0/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F1X4/M1a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F1X4/M2a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F1X4/M7/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F1X4/M8/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F1X4/M8a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F3X4/M0/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F3X4/M1a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F3X4/M2a/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F3X4/M7/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F3X4/M8/codeml.ctl
codeml $PROJECT_DIR/fragments/fragment_3/codeml/F3X4/M8a/codeml.ctl
