awk script


 
Thread Tools Search this Thread
Top Forums Shell Programming and Scripting awk script
# 1  
Old 02-27-2015
awk script

Hi,

I wrote an awk script to compare intervals between two files
Code:
awk 'FNR==NR {a[$1]; next} $1 in a'

My input files I am comparing for overlaps are arranged like this:
Code:
chr1    23110    23140 
chr1    23120    23371    MIR    SINE/MIR

But my script only gives the overlap of one file over the other.
What do I have to do to get an output like this:
Code:
chr1    23110    23140  chr1    23120    23371    MIR  SINE/MIR

regards.

Moderator's Comments:
Mod Comment Please use code tags for your code and data next time

Last edited by vbe; 02-27-2015 at 11:16 AM..
# 2  
Old 02-27-2015
Hi,
The input files are attached
# 3  
Old 02-27-2015
Quote:
Originally Posted by genebuster
Hi,
The input files are attached
---------- Post updated at 10:43 AM ---------- Previous update was at 10:42 AM ----------

Hi,

I am trying to find out overlapping regions in two given files, so I wrote a awk script
Code:
 awk 'FNR==NR {a[$1]; next} $1 in a'

My first input file is arranged like this:
Code:
 
chr7 751143 766356
chr7 825291 826114
chr7 826118 857092
chr7 912969 916281
chr7 935060 938556
chr7 995044 1009012
chr7 1015146 1022846
chr7 1028459 1037259

and the second input file like this:
Code:
 
chr1 10001 10468 (CCCTAA)n Simple_repeat
chr1 10469 11447 TAR1 Satellite/telo
chr1 11505 11675 L1MC5a LINE/L1
chr1 11678 11780 MER5B DNA/hAT-Charlie
chr1 15265 15355 MIR3 SINE/MIR
chr1 16713 16749 (TGG)n Simple_repeat
chr1 18907 19048 L2a LINE/L2
chr1 19972 20405 L3 LINE/CR1
chr1 20531 20679 Plat_L3 LINE/CR1

If I want to find out the overlapping regions between the two files having an output like this which contains elements from the first and the second file :
Code:
chr7 751143 766356 chr7 751143 766356 MIR3 SINE/MIR

How should I write the script ?
The output is just an example I have given to illustrate the overlaps between two files

Last edited by Scrutinizer; 02-27-2015 at 11:47 AM.. Reason: code tags
# 4  
Old 02-27-2015
Hello genebuster,

Kindly use code tags while using commands/codes in your posts. Not clear still, could you please let me know what is the common attribute between 2 files which we will take and start making our code? Also let us know on what basis you want output too?


Thanks,
R. Singh
# 5  
Old 02-27-2015
This is already implemented in the bedtools package, see tutorials below...

Bedtools intersect
# 6  
Old 02-27-2015
I'm sorry, I can't see any overlap between the numbers in the two files nor can I see how the desired output is obtained...
Login or Register to Ask a Question

Previous Thread | Next Thread

10 More Discussions You Might Find Interesting

1. UNIX for Beginners Questions & Answers

Shell script to call and sort awk script and output

I'm trying to create a shell script that takes a awk script that I wrote and a filename as an argument. I was able to get that done but I'm having trouble figuring out how to keep the header of the output at the top but sort the rest of the rows alphabetically. This is what I have now but it is... (1 Reply)
Discussion started by: Eric7giants
1 Replies

2. Shell Programming and Scripting

awk script to call another script based on second column entry

Hi I have a text file (Input.txt) with two column entries separated by tab as given below: aaa str1 bbb str2 cccccc str3 dddd str4 eee str3 ssss str2 sdf str3 hhh str1 fff str2 ccc str3 ..... ..... ..... (1 Reply)
Discussion started by: my_Perl
1 Replies

3. UNIX for Dummies Questions & Answers

Passing shell script parameter value to awk command in side the script

I have a shell script (.sh) and I want to pass a parameter value to the awk command but I am getting exception, please assist. diff=$1$2.diff id=$2 new=new_$diff echo "My id is $1" echo "I want to sync for user account $id" ##awk command I am using is as below cat $diff |... (1 Reply)
Discussion started by: Sarita Behera
1 Replies

4. Post Here to Contact Site Administrators and Moderators

Unable to pass shell script parameter value to awk command in side the same script

Variable I have in my shell script diff=$1$2.diff id=$2 new=new_$diff echo "My id is $1" echo "I want to sync for user account $id" ##awk command I am using is as below cat $diff | awk -F'~' ''$2 == "$id"' {print $0}' > $new I could see value of $id is not passing to the awk... (0 Replies)
Discussion started by: Ashunayak
0 Replies

5. Shell Programming and Scripting

Calling shell script within awk script throws error

I am getting the following error while passing parameter to a shell script called within awk script. Any idea what's causing this issue and how to ix it ? Thanks sh: -c: line 0: syntax error near unexpected token `newline' sh: -c: line 0: `./billdatecalc.sh ... (10 Replies)
Discussion started by: Sudhakar333
10 Replies

6. Shell Programming and Scripting

Passing awk variable argument to a script which is being called inside awk

consider the script below sh /opt/hqe/hqapi1-client-5.0.0/bin/hqapi.sh alert list --host=localhost --port=7443 --user=hqadmin --password=hqadmin --secure=true >/tmp/alerts.xml awk -F'' '{for(i=1;i<=NF;i++){ if($i=="Alert id") { if(id!="") if(dt!=""){ cmd="sh someScript.sh... (2 Replies)
Discussion started by: vivek d r
2 Replies

7. Shell Programming and Scripting

Help: How to convert this bash+awk script in awk script only?

This is the final first release of the dynamic menu generator for pekwm (WM). #!/bin/bash function param_val { awk "/^${1}=/{gsub(/^${1}="'/,""); print; exit}' $2 } echo "Dynamic {" for CF in `ls -c1 /usr/share/applications/*.desktop` do name=$(param_val Name $CF) ... (3 Replies)
Discussion started by: alexscript
3 Replies

8. Shell Programming and Scripting

Call shell script function from awk script

hi everyone i am trying to do this bash> cat abc.sh deepak() { echo Deepak } deepak bash>./abc.sh Deepak so it is giving me write simply i created a func and it worked now i modified it like this way bash> cat abc.sh (2 Replies)
Discussion started by: aishsimplesweet
2 Replies

9. Shell Programming and Scripting

want to pass parameters to awk script from shell script

Hello, I have this awk script that I want to execute by passing parameters through a shell script. I'm a little confused. This awk script removes duplicates from an input file. Ok, so I have a .sh file called rem_dups.sh #!/usr/bin/sh... (4 Replies)
Discussion started by: script_op2a
4 Replies

10. Shell Programming and Scripting

create a shell script that calls another script and and an awk script

Hi guys I have a shell script that executes sql statemets and sends the output to a file.the script takes in parameters executes sql and sends the result to an output file. #!/bin/sh echo " $2 $3 $4 $5 $6 $7 isql -w400 -U$2 -S$5 -P$3 << xxx use $4 go print"**Changes to the table... (0 Replies)
Discussion started by: magikminox
0 Replies
Login or Register to Ask a Question