From 41d57ed01da5a967ac25d80aee368c8995226e42 Mon Sep 17 00:00:00 2001 From: Yutian Feng Date: Thu, 18 Oct 2018 17:47:38 -0400 Subject: [PATCH] Delete mlsa.sh --- mlsa.sh | 47 ----------------------------------------------- 1 file changed, 47 deletions(-) delete mode 100644 mlsa.sh diff --git a/mlsa.sh b/mlsa.sh deleted file mode 100644 index fa8fb94..0000000 --- a/mlsa.sh +++ /dev/null @@ -1,47 +0,0 @@ -#assumes blast+ and iqtree, and muscle. also requires supermatrix.py and wrapper_supermatrix.py in the a usable path or PATH - -for f in * (names of genomes in fasta).fasta; do makeblastdb -i $f -dbtype nucl; done -#need names of genomes in nucleotide fasta format - -for x in *.blast; do while IFS=$'\t' read -r -a myArray; do NAME="${myArray[1]}"; START="${myArray[8]}"; STOP="${myArray[9]}"; done < "$x" ; if [ "$START" -gt "$STOP" ]; then STRAND=2; RANGE="$STOP"-"$START"; else STRAND=1; RANGE="$START"-"$STOP"; fi; - g=${x%.blast}; fastacmd -d /home/yuf17006/Frankia_project/Small_Data_tests/Genomes/"$g" -p F -s "$NAME" -L "$RANGE" -S "$STRAND" -o "$x".grab; done; -#grabs top scoring hit for each gene outputs them to the same gene file - -#align them all with muscle -for g in *.grab; do muscle -in $g -out $g.align -maxiters 4; done - -#all taxa must have the same names in every gene file - - - -#instructions for using supermatrix concatenater -#Attached are two py scripts. On the cluster, create a directory in your home folder called scripts. copy the two files there. -#In the texteditor of your choice, edit the file .bash_profile (dot files are usually not listed ...) -#In that file should be a description of your PATH. Something like -#PATH=$PATH:$HOME/bin:/opt/454/bin -#You need to add the script folder to the path. in the example above, after editing the line would read -#​P​ATH=$PATH:$HOME/bin:/opt/454/bin:$HOME/scripts -​#Then, -#create a folder that contains the alignments you want to concatenate -#Create a text file that lists all the names of the individual alignments. -#save that file (possible Alist.txt) -#Invoke the wrapper by typing -#wrapper_supermatrix.py Alist.txt -#​The contents of alist would look something like this: -#list1_aligned.fst -#list2_aligned.fst -#list3_aligned.fst -#list4_aligned.fst -#list5_aligned.fst -#list7_aligned.fst -#​The program creates a lot of output on the screen and in the directory where your alignments were. -#Essentially if you ave 31 files, the program runs 30 time, creating a concatenated alignment every time. Outfile30 would be the concatenation you want to work with. But check the logs, the missingseqs and partition lists files, in case something went wrong ... And, obviously, look at the final alignment in seaview or similar. -​ - -wrapper_supermatrix.py Alist.txt - -#take the final concatenated alignments and turn them into codons, proteins or leave them -#can use as input into iqtree -iqtree -st (codon, AA, or NT) -bb 1000+ -nt AUTO - -#view tree \ No newline at end of file