Blast results processing lines bigger than 60
0
0
Entering edit mode
2.8 years ago
hdtms ▴ 20

Hi,

I have a doubt processing blast results in linux command line with a column I added, I want the cov result to be bigger than 60. But I can't make this command work. All works except {cov > 60} I also tried (cov > 60) but it doesn't select the results bigger than 60. I must also say that I started by thinking that the cov result would be in decimal like 0.63 but I noticed it gives results in 63.xxx.

Any suggestion would be greatly appreciated.

for file in *.join; do echo $file ; awk '{cov = $4/$13}{cov > 60}{print $1, $2, $3}' $file >> ../MUSTARD_NEW_RESULTS_FILTERED/$file.filtered; done
blast unix commands awk result • 792 views
ADD COMMENT
1
Entering edit mode

What format is your blast output in?

ADD REPLY
0
Entering edit mode

Output format 6

ADD REPLY
0
Entering edit mode

I also rechecked the command and it turns out that the only part that is working is the print command.

ADD REPLY

Login before adding your answer.

Traffic: 2093 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6