column 2 describes start, column 3 end and column 4 a value.
I want to print al lines between the lines that have a zero value PLUS have a 4 in their column 4. And only print the first start and last end.
So for this example, I would print out:
Code:
chr1 569791 570061
Because this is the only region which is between two zeros and does contain a 4.
Last edited by Scott; 10-17-2012 at 04:27 PM..
Reason: Code tags, please...
Hi, Thanks so much! That works exactly how it should. Could you explain your code?
awk '$4==0{$0=""}{$1=$1}1' file
replace lines with zero in fourth field with empty lines
awk -vRS=
process the file using paragraphs (lines enclosed between empty lines)
'{for (i=4;i<=NF;i+=4) if ($i==4) print $1,$2,$(NF-1)}'
check every fourth field of each paragraph (fourth field of each line) and if it equals "4", then print first, second and second to last field from that paragraph.
---------- Post updated at 02:55 PM ---------- Previous update was at 02:54 PM ----------
Hi,
I would like to know how can I get a mean score value by ID over a defined genomic region. Here it is an example:
file1
12 100 103 id1
12 110 112 id1
12 200 203 id2
file2
12 100 101 1
12 101 102 0.8
12 102 103 0.7
12 110 111 2.5
12 111 112 2.8
12 200 201 10.1
12 201 202... (7 Replies)
Hello, I am trying to concatenate two strings by merging the overlapped region. E.g.
Seq1=ACGTGCCC
Seq2=CCCCCGTGTGTGT
Seq_merged=ACGTGCCCCCGTGTGTGTFunction strcat(char *dest, char *src) appends the src string to the dest string, ignoring the overlapped parts (prefix of src and suffix of dest).... (30 Replies)
Input file
2 5 ASFGEWTEWRQWEQ
10 20 QEWIORUEIOUEWORUQWEQWRQRQWGQWGFQ
1 6 WRQTQWTQTQWTQT
Desired output file
2 5 ASFGEWTEWRQWEQ
10 20 QEWIORUEIOUEWORUQWEQWRQRQWGQWGFQ
1 6 WRQTQWTQTQWTQT
Column 1 is the start region of underline the text in column 3;
Column 2 is the end region of... (13 Replies)
Dears my rootvg is missed up i can not extend the /opt
as soon as i try to extend the Filesystem its give me that there is not enough space .
as there any way to change the REGION of the LVs in HDISK1 ?
lspv -p hdisk0
hdisk0:
PP RANGE STATE REGION LV NAME TYPE ... (8 Replies)
Hi guys,
I hope everybody is doing fine. I have written this small program which solves the critical region problem. Only on of the two threads can make changes to a common variable called counter. I am using two semaphores, is it possible to write the same program using only one semaphore? Here... (0 Replies)
Hello all,
I have a question about what you think the best practice is to determine what region you are running on when you have a system setup with a DEV/TEST, QA, and PROD regions running the same scripts in all.
So, when you run in DEV, you have a different directory structure, and you... (4 Replies)
Hi:
Let's suppose I want to replace all the | by > ONLY when | is between . Usually (and it works) I would do something like
sed -e 's/\(\*\)|\(*\]\)/\1>\2/g'
where I have to "save" some portions of the matched region and use them with the \n metacharacter. I was wondering if I could... (2 Replies)
10-17-2012
