Create a single bash script that does the following:
a. Print out the number of occurrences for each motif that is found in the bacterial genome and output to a file called motif_count.txt
b. Create a fasta file for each motif (so 3 in total) which contains all of the genes and their corresponding sequences that have that motif. Each file should be named after the motif (ie ATG.txt) and outputted to a new directory called motifs
motif file is a txt file with these motifs : ATG, GGGGG, ATTTT
the bacterial genome file is fasta file with the following lines
Moderator's Comments:
edit by bakunin: please use CODE-tags to let data, code and terminal output stand out. Thank you.
I have a very basic bash shell script, which has many "while... done; for .... done" loop clauses, like the following
~~
#!/bin/bash
while blablalba; do
....
done < /tmp/file
for line in `cat blablabla`; do grep $line /tmp/raw ; done > /tmp/1;
while blablalba2; do
....
done <... (2 Replies)
SEND_MESSAGE=test
echo $SEND_MESSAGE
if
then
echo `date` > update_dt_ccaps.lst
echo "The file transfer failed" >> update_dt_ccaps.lst
SEND_MESSAGE=false
fi
The above code is showing error in bash shell as :
./test: line 5: [: test: integer expression expected
... (2 Replies)
I'm trying to search all .odt files in a directory for a string in the text of the file.
I've found a bash script that works, except that it can't handle whitespace in the filenames.
#!/bin/bash
if ; then
echo "Usage: searchodt searchterm"
exit 1
fi
for file in $(ls *.odt); do
... (4 Replies)
Hi,
I am working on bash script after a long time. I am getting error near done statement while running a for loop snippet. The error says "Syntax error near unexpcted token 'done'"
please suggest what could be wrong. here is the snippet
elements=${#option_arr} //an array of values... (1 Reply)
Having issues with an expect script. I've been scripting bash, python, etc... for a couple years now, but just started to try and use Expect. Trying to create a script that takes in some arguments, and then for now, just runs a pwd command(for testing, final will be command I pass).
Here is... (0 Replies)
Hi Folks,
I have a loop that goes through an array and the output is funky.
sample:
array=( 19.239.211.30 )
for i in "${array}"
do
echo $i
iperf -c $i -P 10 -x CSV -f b -t 50 | awk 'END{print '$i',$6}' >> $file
done
Output:
19.239.211.30
19.2390.2110.3 8746886
seems that when... (2 Replies)
I have the following while loop that I put in a script, demo.sh:
while read rna; do
aawork=$(echo "${rna}" | sed -n -e 's/\(...\)\1 /gp' | sed -f rna.sed)
echo "$aawork" | sed 's/ //g'
echo "$aawork" | tr ' ' '\012' | sort | sed '/^$/d' | uniq -c | sed 's/*\(*\) \(.*\)/\2: \... (3 Replies)
Hi Everybody,
I'm a newbie to shell scripting, and I'd appreciate some help. I have a bunch of .txt files that have some unwanted content. I want to remove lines 1-3 and 1028-1098.
#!/bin/bash
for '*.txt' in <path to folder>
do
sed '1,3 d' "$f";
sed '1028,1098 d' "$f";
done
I... (2 Replies)
Dear all,
Please help with the following.
I have a file, let's call it data.txt, that has 3 columns and approx 700,000 lines, and looks like this:
rs1234 A C
rs1236 T G
rs2345 G T
Please use code tags as required by forum rules!
I have a second file, called reference.txt,... (1 Reply)
Hi everyone
I need some help
I want to create an script which does some processing
it takes the two arguments 201901010000 and 201901020200 - so YYYMMDDHHMM
I want to split processing into hours from start until end,
I dont get why this works but when I add to a future variable... (1 Reply)
Discussion started by: kl1ngac1k
1 Replies
LEARN ABOUT DEBIAN
bp_revtrans-motif
BP_REVTRANS-MOTIF(1p) User Contributed Perl Documentation BP_REVTRANS-MOTIF(1p)NAME
revtrans-motif - Reverse translate a Profam-like protein motif
VERSION
Version 0.01
SYNOPSIS
From a file:
revtrans-motif.pl -i motifs.txt
Using pipes:
revtrans-motif.pl < motifs.txt > output.txt
Using interactively at the command prompt:
$ revtrans-motif.pl
MAAEEL[VIKP]
1. ATGGCNGCNGARGARYTNVHN
[^P]H(IW){2,3}
2. NDNCAY(ATHTGG){2,3}
DESCRIPTION
This script takes a protein motif as input and returns a degenerate oligonucleotide sequence corresponding to it. The main reason for doing
this is to design degenerate primers that amplify a given sequence pattern.
The input motif consists of a string of one-letter residues, with any of the following syntactic elements:
[...] : Redundant position.
A position in which more than one residue is allowed. Example:
[TS]YW[RKSD]
^^ ^^^^
[^...] : Negated position.
A position in which any residue is allowed, saved for those between brackets. Example:
[^PW]MK[LAE]
^^
(...){n,m,...} : Repeated motif.
A motif that is repeated n or m times. It can have any of the previous syntactic elements. Example:
A[SN]C(TXX){2,4,8}
^^^
The allowed letters are those that correspond to the 20 natural aminoacids, plus:
B = N + D
Z = Q + E
X = All
OPTIONS -i input-file:
A file with a list of motifs to reverse translate.
-h
Display this help message.
AUTHOR
Bruno Vecchi, "vecchi.b at gmail.com"
BUGS
Please report any bugs or feature requests to "vecchi.b at gmail.com"
COPYRIGHT & LICENSE
Copyright 2009 Bruno Vecchi, all rights reserved.
This program is free software; you can redistribute it and/or modify it under the same terms as Perl itself.
perl v5.14.2 2012-03-02 BP_REVTRANS-MOTIF(1p)