Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- for X in $(ls);do N=$(grep -c "." $X);
- for W in $(seq 1 $N);do Y=$(sed -n "$W"p $X);
- AN=$(echo $Y | cut -f1 -d "+");
- AN=$(echo $AN | tr '[:upper:]' '[:lower:]');
- Seq=$(echo $Y | cut -f2 -d "+");
- extractseq -sequence /Users/remironzano/Desktop/Bioinfo/BioInfoProjet/Unisequence.fasta_avec_nom_accessionnumber/$AN.fasta -regions $Seq -outseq /Users/remironzano/Desktop/Bioinfo/BioInfoProjet/Resultats_seq_signal.txt; done; done
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement