Exploring DNA Motifs with Custom Bash Script

gene_x 0 like s 658 view s

Tags: bash, Motif Discovery

#!/bin/bash

#./search_motif4.sh test1.fasta GRG 5

if [ $# -ne 3 ]; then
  echo "Usage: $0 <fasta_file> <motif> <context>"
  exit 1
fi

fasta_file=$1
motif=$2
context=$3

motif_regex=$(echo $motif | sed 's/R/[AG]/g' | sed 's/Y/[CT]/g' | sed 's/S/[GC]/g' | sed 's/W/[AT]/g' | sed 's/K/[GT]/g' | sed 's/M/[AC]/g' | sed 's/B/[CGT]/g' | sed 's/D/[AGT]/g' | sed 's/H/[ACT]/g' | sed 's/V/[ACG]/g')

grep -B1 -A1 -i -E -o ".{0,$context}${motif_regex}.{0,$context}" $fasta_file

like unlike

点赞本文的读者

还没有人对此文章表态

本文有评论

没有评论

Exploring DNA Motifs with Custom Bash Script

本文有评论

看文章，发评论，不要沉默

最受欢迎文章

最新文章

最多评论文章

推荐相似文章