-
Notifications
You must be signed in to change notification settings - Fork 9
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
19 changed files
with
23 additions
and
23 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -19,7 +19,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 0a - Download & prepare data * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 0b - Prepare pseudoreference * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -9,7 +9,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 0c - Setup R packages on Metacentrum * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 01 - Raw data processing * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# * based on Weitemier et al. (2014), Applications in Plant Science 2(9): 1400042* | ||
|
@@ -179,12 +179,12 @@ do | |
#Adapter and general quality trimming | ||
echo "Adapter removal & quality trimming..." | ||
if [[ $location == "1" ]]; then | ||
java -jar /software/trimmomatic-0.32/dist/jar/trimmomatic-0.32.jar PE -phred33 $file-noPhiX_1.fq.gz $file-noPhiX_2.fq.gz $file-1P $file-1U $file-2P $file-2U ILLUMINACLIP:../NEBNext-PE.fa:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:36 >Trimmomatic.log 2>&1 | ||
java -jar /software/trimmomatic-0.32/dist/jar/trimmomatic-0.32.jar PE -phred33 $file-noPhiX_1.fq.gz $file-noPhiX_2.fq.gz $file-1P $file-1U $file-2P $file-2U ILLUMINACLIP:../${adapterfile}:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:36 >Trimmomatic.log 2>&1 | ||
#trimmomatic PE -phred33 $file-noPhiX_1.fq.gz $file-noPhiX_2.fq.gz $file-1P $file-1U $file-2P $file-2U ILLUMINACLIP:../NEBNext-PE.fa:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:36 >Trimmomatic.log 2>&1 | ||
elif [[ $location == "2" ]]; then | ||
java -d64 -server -XX:MaxHeapSize=10g -jar ../trimmomatic-0.33.jar PE -phred33 $file-noPhiX_1.fq.gz $file-noPhiX_2.fq.gz $file-1P $file-1U $file-2P $file-2U ILLUMINACLIP:../NEBNext-PE.fa:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:36 >Trimmomatic.log 2>&1 | ||
java -d64 -server -XX:MaxHeapSize=10g -jar ../trimmomatic-0.33.jar PE -phred33 $file-noPhiX_1.fq.gz $file-noPhiX_2.fq.gz $file-1P $file-1U $file-2P $file-2U ILLUMINACLIP:../${adapterfile}:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:36 >Trimmomatic.log 2>&1 | ||
else | ||
java -jar ../trimmomatic-0.33.jar PE -phred33 $file-noPhiX_1.fq.gz $file-noPhiX_2.fq.gz $file-1P $file-1U $file-2P $file-2U ILLUMINACLIP:../NEBNext-PE.fa:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:36 >Trimmomatic.log 2>&1 | ||
java -jar ../trimmomatic-0.33.jar PE -phred33 $file-noPhiX_1.fq.gz $file-noPhiX_2.fq.gz $file-1P $file-1U $file-2P $file-2U ILLUMINACLIP:../${adapterfile}:2:30:10 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:36 >Trimmomatic.log 2>&1 | ||
fi | ||
#Copy trimmed reads to home | ||
if [[ $location == "1" ]]; then | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 03 - Process consensus after mapping, make pslx files * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# * based on Weitemier et al. (2014), Applications in Plant Science 2(9): 1400042* | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -20,7 +20,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 04 - Process pslx files * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# * based on Weitemier et al. (2014), Applications in Plant Science 2(9): 1400042* | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -8,7 +8,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 04b - Put exons to correct reading frame and translate * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -19,7 +19,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 05 - Missing data handling * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
@@ -243,7 +243,7 @@ do | |
#Calculate percentage of missing data as integer (for evaluation in bash) | ||
percentage=$(echo -e "scale=0;($count * 100) / $y" | bc) | ||
echo -e "$header\t$percentageprint" >> ${file}_${MISSINGPERCENT}percN.fas | ||
if [ "$percentage" -lt "$MISSINGPERCENT" ]; then | ||
if [ "$percentage" -le "$MISSINGPERCENT" ]; then | ||
echo "$line" >> ${file}_modif${MISSINGPERCENT}.fas | ||
else | ||
echo -e "Deleting $header with $percentageprint of missing data" | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -19,7 +19,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 06a2 - summary of RAxML gene trees * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -17,7 +17,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 06a - RAxML gene tree building * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 06b - FastTree gene tree building * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -17,7 +17,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 07 - Root gene trees * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 08a - Astral species tree * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -19,7 +19,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 08b - Astrid species tree * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 08c - MRL species tree * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 08e - concatenated species tree * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -8,7 +8,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 08f - ExaML concatenated tree * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -18,7 +18,7 @@ | |
# ******************************************************************************** | ||
# * HybPhyloMaker - Pipeline for Hyb-Seq data processing and tree building * | ||
# * Script 09 - Update trees * | ||
# * v.1.4.1 * | ||
# * v.1.4.2 * | ||
# * Tomas Fer, Dept. of Botany, Charles University, Prague, Czech Republic, 2017 * | ||
# * [email protected] * | ||
# ******************************************************************************** | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -8,7 +8,7 @@ | |
# Tomas Fer, 2017 # | ||
# [email protected] # | ||
# https://github.com/tomas-fer/HybPhyloMaker # | ||
# v.1.4.0 # | ||
# v.1.4.2 # | ||
########################################################################################################################## | ||
|
||
#Carefully set your distribution | ||
|