Skip to content
Snippets Groups Projects
Commit 8c35a32f authored by peguerin's avatar peguerin
Browse files

first contact with this script lol

parent 70794ce9
Branches
No related tags found
No related merge requests found
...@@ -77,9 +77,9 @@ for sample in `ls $main_dir/"$pref"_sample_*.fasta`; ...@@ -77,9 +77,9 @@ for sample in `ls $main_dir/"$pref"_sample_*.fasta`;
do do
sample_sh="${sample/.fasta/_cmd.sh}" sample_sh="${sample/.fasta/_cmd.sh}"
echo "bash "$sample_sh >> $all_samples_parallel_cmd_sh echo "bash "$sample_sh >> $all_samples_parallel_cmd_sh
# Drplication des reads en squences uniques # Dereplicate sequences with tally
dereplicated_sample="${sample/.fasta/.uniq.fasta}" dereplicated_sample="${sample/.fasta/.uniq.fasta}"
echo "/usr/bin/time $tally -i "$sample" --fasta-in --fasta-out -o "$dereplicated_sample > $sample_sh; echo "$tally -i "$sample" --fasta-in --fasta-out -o "$dereplicated_sample > $sample_sh;
# On garde les squences de plus de 20pb sans bases ambigues # On garde les squences de plus de 20pb sans bases ambigues
good_sequence_sample="${dereplicated_sample/.fasta/.l20.fasta}" good_sequence_sample="${dereplicated_sample/.fasta/.l20.fasta}"
echo "$obigrep -s '^[ACGT]+$' -l 20 "$dereplicated_sample" > "$good_sequence_sample >> $sample_sh echo "$obigrep -s '^[ACGT]+$' -l 20 "$dereplicated_sample" > "$good_sequence_sample >> $sample_sh
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment