I have this shell script that I wrote to check an input file to see if it is empty or not, and then clean the file from any line that starts with the sign "<" (without quotation marks" and then spell the number of line of the file, and the empty lines, too. The script then will create two output files, DNA.out and RNA.out.
First I get an error message at that says:
Code:
./script.sh: line 3: [: dna_input.txt: integer expression expected
but it gives me the results I want.
Here is the code:
Code:
#!/bin/bash
#check to see if there is an input file:
if [ $1 -lt 1 ]
then
echo "Usage: $0 file ..."
exit 1
fi
#Check if the file is empty or not
file=$1
if [[ -s $1 ]]
then
echo ""
echo "**** $file has data."
echo "Number of non-empty lines:"
grep -cve '^\s*$' $file
echo "Mumber of empty lines:"
grep -ce '^\s*$' $file
#grep -cvP '^\s*$' $file -- above line originally was like this one
echo ""
cat $1 | sed 's/>/\n>/g' > temp1.txt
cat temp1.txt | sed '/^>/ d' > temp2.txt
#Remove duplicate empty lines:
awk '!NF{if(++n <=1) print; next}; {n=0; print}' < temp2.txt > DNA.out
echo ""
#convert to mRNA and remove temp1, temp2 to avoid confusion
tr ACGT UGCA < DNA.out > RNA.out
rm temp1.txt temp2.txt
else
echo "**** $1 has no data, or file does not exist."
echo "**** done!"
echo ""
fi;
After I get the two output files (DNA.out and RNA.out) I use another script to convert the contents of these two files into Amino Acids. The conversion script is:
Code:
#!/bin/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/[ ]*\([0-9]*\) \(.*\)/\2: \1/'
done
Can my scripts be combined together in just one script? My script gives the frequencies of each Amino Acid. Also, I ask to learn how to add a while loop in a script if it takes an input of its own.
I like my code to run in /tmp file when I am using a static files.
Code:
#!/bin/bash
for i in /home/*/Desktop/test.txt; do cp $i /tmp/test.txt && echo "Usage: $i file ... " || echo "File is not where you think it should be";
sleep 3
done
exit 0
Consider trapping ( think that is what it called. Do a google) the file. So after checking that's the file exists with code above, follow up with ...
Code:
cd /tmp
sudo ./test.txt
And then within test.txt. Have this line of code. Which also checks the file was placed in /tmp. Then runs your test.txt file.
This is called trapping and insure the program will stop running when it is completed. You can also run sub routines under the trap before running main code body.
Okay I want to try my luck at fixing bugs for the Fedora OS, but I guess this question deals with any Linux distro or any open source OS for that matter.
I want to know how I can start fixing bugs on the OS level. For example the particular bug that I want to target is this logout bug
I mean... (6 Replies)
can someone please help me spot and fix the issue with the following code:
awk -F, -v SEARCHPATT="(Wed|Tue)" -v ADDISTR="Mon|Tue|Wed|Thu|Fri|Sat|Sun" -vVF="$VALFOUND"
"BEGIN{ {D = D = 1
D = D = 2
}
$0 ~ "," VF "," {L = 1 ... (9 Replies)
Hey guys,
I am trying to make print a pattern with * on a 10*10 two dimensional array in a for loop and I want the incoming 10*10 to overlap the previous 10*10 so that the * look like it is moving. is there a way to fix the screen position?
ever time it prints a 10*10 the screen moves.
... (3 Replies)
Hi
I use Rapid 7 to check some servers ( AIX 5.3 ) for security problems. There are 2 problems I don't know to deal with
1.
Problem : TCP Sequence Number Approximation Vulnerability
Solution :
_Enable TCP MD5 Signature
2.
Problem : HTTP Basic Authentication Enable
Solution :
_ Use... (5 Replies)
1. The problem statement, all variables and given/known data:
I need help I get a variant of syntax errors when compiling my script to maintain a database. It's a simple database meant to create/view/maintain vehicles.
2. Relevant commands, code, scripts, algorithms:
my if statements have... (5 Replies)
Hi Experts
I've problem in a my office server (solaris 10 - x86) version. x4600 M2 hardware
This system is getting rebooted because of zfs bug
I've applied patch using live upgrade with live new environment created and
applied the patch which oracle suggested(
144501-19), it asks for... (3 Replies)
Is there a way to fix the width of the word being printed to a file? I am trying to create an output to a file with columns , like a spread sheet. I have used "\t" to adjust the columns but still it does not show well in the file, mainly due to the variable length values in the column so \t does... (1 Reply)
I am trying to replace the value of $f3 but its not working . I don't know what I am missing here .
cat dim_copy.20080516.sql | grep -i "create view" | grep -v OPSDM002 | while read f1 f2 f3 f4 f5 f6 f7 f8 f9
do
echo " $f3 "
sed -e... (13 Replies)