Extract Big and continuous regions Post: 302969222

Sponsored Content

Top Forums Shell Programming and Scripting Extract Big and continuous regions Post 302969222 by amrutha_sastry on Sunday 20th of March 2016 04:13:57 AM

03-20-2016

Registered User

Extract Big and continuous regions

Hi all,
I have a file like this I want to extract only those regions which are big and continous

Code:

chr1    3280000 3440000
chr1    3440000 3920000
chr1    3600000 3920000 # region coming within the 3440000 3920000. so i don't want it to be printed in output
chr1    3920000 4800000 
chr1    4080000 4800000 # region coming within the 3920000 4800000 . so i don't want it to be printed in output
chr1    4160000 4360000 # region coming within the 3920000 4800000 . so i don't want it to be printed in output
chr1    4160000 4760000 # region coming within the 3920000 4800000 . so i don't want it to be printed in output
chr1    4360000 4760000 # region coming within the 3920000 4800000 . so i don't want it to be printed in output
chr1    4800000 4920000 # region coming within the 4800000 6160000 . so i don't want it to be printed in output
chr1    4800000 5200000 # region coming within the 4800000 6160000 . so i don't want it to be printed in output
chr1    4800000 5280000 # region coming within the 4800000 6160000 . so i don't want it to be printed in output
chr1    4800000 6160000 
chr1    4920000 5200000 # region coming within the 4800000 6160000 . so i don't want it to be printed in output
chr1    5280000 6160000 # region coming within the 4800000 6160000 . so i don't want it to be printed in output

Finally the below regions could be printed in output

Code:

chr1    3280000 3440000
chr1    3440000 3920000
chr1    3920000 4800000
chr1    4800000 6160000

So can anyone help me to get that

amrutha_sastry

View Public Profile for amrutha_sastry

Find all posts by amrutha_sastry

8 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

How to view a big file(143M big)

1 . Thanks everyone who read the post first. 2 . I have a log file which size is 143M , I can not use vi open it .I can not use xedit open it too. How to view it ? If I want to view 200-300 ,how can I implement it 3 . Thanks

2. UNIX for Dummies Questions & Answers

How big is too big a config.log file?

I have a 5000 line config.log file with several "maybe" errors. Any reccomendations on finding solvable problems?

3. Shell Programming and Scripting

awk: union regions

Hi all, I have difficulty to solve the followign problem. mydata: StartPoint EndPoint 22 55 2222 2230 33 66 44 58 222 240 11 25 22 60 33 45 The union of above...

4. Shell Programming and Scripting

Extract certain entries from big file:Request to check

Hi all I have a big file which I have attached here. And, I have to fetch certain entries and arrange in 5 columns Name Drug DAP ID disease approved or notIn the attached file data is arranged with tab separated columns in this way: and other data is...

5. UNIX for Dummies Questions & Answers

extract regions of file based on start and end position

Hi, I have a file1 of many long sequences, each preceded by a unique header line. file2 is 3-columns list: headers name, start position, end position. I'd like to extract the sequence region of file1 specified in file2. Based on a post elsewhere, I found the code: awk...

6. Shell Programming and Scripting

Obtain the names of the flanking regions

Hi I have 2 files; usually the end position in the file1 is the start position in the file2 and the end position in file2 will be the start position in file1 (flanks) file1 Id start end aaa1 0 3000070 aaa1 3095270 3095341 aaa1 3100822 3100894 aaa1 ...

7. Shell Programming and Scripting

Extract certain columns from big data

The dataset I'm working on is about 450G, with about 7000 colums and 30,000,000 rows. I want to extract about 2000 columns from the original file to form a new file. I have the list of number of the columns I need, but don't know how to extract them. Thanks!

8. Shell Programming and Scripting

Want to extract certain lines from big file

Hi All, I am trying to get some lines from a file i did it with while-do-loop. since the files are huge it is taking much time. now i want to make it faster. The requirement is the file will be having 1 million lines. The format is like below. ##transaction, , , ,blah, blah...

LEARN ABOUT ULTRIX

xintersectregion

XIntersectRegion(3X11)						     MIT X11R4						    XIntersectRegion(3X11)

Name
       XIntersectRegion, XUnionRegion, XUnionRectWithRegion, XSubtractRegion, XXorRegion, XOffsetRegion, XShrinkRegion - region arthmetic

Syntax
       XIntersectRegion(sra, srb, dr_return)
	  Region sra, srb, dr_return;

       XUnionRegion(sra, srb, dr_return)
	  Region sra, srb, dr_return;

       XUnionRectWithRegion(rectangle, src_region, dest_region_return)
	  XRectangle *rectangle;
	  Region src_region;
	  Region dest_region_return;

       XSubtractRegion(sra, srb, dr_return)
	  Region sra, srb, dr_return;

       XXorRegion(sra, srb, dr_return)
	  Region sra, srb, dr_return;

       XOffsetRegion(r, dx, dy)
	  Region r;
	  int dx, dy;

       XShrinkRegion(r, dx, dy)
	  Region r;
	  int dx, dy;

Arguments
       dest_region_return
		 Returns the destination region.

       dr_return Returns the result of the computation.  ds Dy move or shrink

       dx
       dy	 Specify the x and y coordinates, which define the amount you want to  the specified region.

       r	 Specifies the region.

       rectangle Specifies the rectangle.

       sra
       srb	 Specify the two regions with which you want to perform the computation.

       src_region
		 Specifies the source region to be used.

Description
       The function computes the intersection of two regions.

       The function computes the union of two regions.

       The function updates the destination region from a union of the specified rectangle and the specified source region.

       The function subtracts srb from sra and stores the results in dr_return.

       The function calculates the difference between the union and intersection of two regions.

       The function moves the specified region by a specified amount.

       The function reduces the specified region by a specified amount.  Positive values shrink the size of the region, and negative values expand
       the region.

See Also
       XCreateRegion(3X11), XDrawRectangle(3X11), XEmptyRegion(3X11)
       X Window System: The Complete Reference, Second Edition, Robert W. Scheifler and James Gettys

															    XIntersectRegion(3X11)