Sponsored Content
Top Forums UNIX for Advanced & Expert Users Map snps into a ref gene file Post 302989923 by Don Cragun on Thursday 19th of January 2017 08:47:30 PM
Old 01-19-2017
And what output are you trying to get from the two sample input files you provided?

What happens if there is no ID in the 1st file that appears in a range specified by the 2nd file?

What happens if there is more than one ID in the 1st file that fits in a range specified by a single line in the 2nd file?

What happens if there is no range in the 2nd file for a position specified in the 1st file?

What have you tried to solve this problem on your own?
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Reading a path (including ref to shell variable) from file

Hi! 1. I have a parameter file containing path to log files. For this example both paths are the same, one is stated directly and the second using env variables. /oracle/admin/orcl/bdump/:atlas:trc:N ${ORACLE_BASE}/admin/${ORACLE_SID}/bdump/:${ORACLE_SID}:trc:N 2. I try to parse the path... (1 Reply)
Discussion started by: lojzev
1 Replies

2. Shell Programming and Scripting

File merging using first column as the ref

I had two files 1.txt 2.txt. I want a 3rd file(o/p) 3.txt like below (using awk) 1.txt 11 a1 12 a2 13 a3 14 a4 15 a5 16 a6 17 a7 18 a8 19 a9 20 a10 2.txt 14 b1 15 b2 16 b3 (8 Replies)
Discussion started by: p_sai_ias
8 Replies

3. Shell Programming and Scripting

Append file from ref file AWK

FILE1 abc:xxx:abc:123:wer:AAA:12 csf:xxx:123:aeg:sar:BBB:13 asq:yer:321:wsa:qqq:CCC:14 FILE2 AAA:12:SET1:R1 AAA:12:SSS1:RR1 AAA:11:SET4:R3 BBB:13:SET2:R2 OUTPUT abc:xxx:abc:123:wer:AAA:12:SET1:R1:SSS1:RR1 csf:xxx:123:aeg:sar:BBB:13:SET2:R2::... (4 Replies)
Discussion started by: greycells
4 Replies

4. Shell Programming and Scripting

Script to search and extract the gene sub-location from gff file.

Hi, my problem is that I have two files. File no. 1 is a gff text file (say gi1) that has gene information like : ******************** gene 39389788..39395643 /gene="RPSA" /note="Derived by automated computational analysis using ... (2 Replies)
Discussion started by: reena2305
2 Replies

5. UNIX for Advanced & Expert Users

cannot find map file

Hi, all: My writed network device driver works fine when the transmitted file is under several MegaBytes, but above this size, especially dozens of or hundreds of MegaBytes, the kernel panic ocurres! I check the kern.log and find this error : 522 Nov 14 19:35:32 liklstar-server kerneNov 14... (2 Replies)
Discussion started by: liklstar
2 Replies

6. Shell Programming and Scripting

File merging using first column as the ref

I had two files 1.txt 2.txt. I want a 3rd file(o/p) 3.txt like below based on the common elements from the first coulmns of 1.txt and 2.txt. 1.txt 11 12 13 14 15 16 17 18 19 20 21 2.txt (6 Replies)
Discussion started by: p_sai_ias
6 Replies

7. UNIX for Dummies Questions & Answers

Breaking a fasta formatted file into multiple files containing each gene separately

Hey, I've been trying to break a massive fasta formatted file into files containing each gene separately. Could anyone help me? I've tried to use the following code but i've recieved errors every time: for i in *.rtf.out do awk '/^>/{f=++d".fasta"} {print > $i.out}' $i done (1 Reply)
Discussion started by: Ann Mc Cartney
1 Replies

8. Shell Programming and Scripting

Extract a string between 2 ref string from a file

Hi, May i ask if someone share some command for extracting a string between 2 ref string in a txt file My objective: i had a file with multiple lines and wants only to extract the string "watch?v=IbkAXOmEHpY" or "watch?v=<11 random character>", when i used "grep 'watch?=*' i got a results per... (4 Replies)
Discussion started by: jao_madn
4 Replies

9. Shell Programming and Scripting

awk to average target and gene

I am trying to modify the awk below to include the gene name ($5) for each target and can not seem to do so. Also, I'm not sure the calculation is right (average of all targets that are the same is $4 using the values in $7)? Thank you :). awk '{if((NR>1)&&($4!=last)){printf("%s\t%f\t%s\n",... (1 Reply)
Discussion started by: cmccabe
1 Replies

10. UNIX for Beginners Questions & Answers

Snps annotation

I have the following Snps data CHROM POS ID chr7 78599583 rs987435 chr15 33395779 rs987436 chr1 189807684 rs987437 chr20 33907909 rs987438 chr12 75664046 rs987439 and the following gene data genename name chrom strand txstart txend... (8 Replies)
Discussion started by: marwah
8 Replies
mlib_VideoColorBlendABGR(3MLIB) 			    mediaLib Library Functions				   mlib_VideoColorBlendABGR(3MLIB)

NAME
mlib_VideoColorBlendABGR, mlib_VideoColorBlendABGR_ResetAlpha - image blend SYNOPSIS
cc [ flag... ] file... -lmlib [ library... ] #include <mlib.h> void mlib_VideoColorBlendABGR(mlib_u32 *dst, const mlib_u32 *src1,const mlib_u32 *src2, mlib_s32 src1_w, mlib_s32 src1_h, mlib_s32 src2_w,mlib_s32 src2_h, mlib_s32 src2_x, mlib_s32 src2_y, mlib_s32 dst_lb,mlib_s32 src1_lb, mlib_s32 src2_lb, mlib_blend src1_blend, mlib_blend src2_blend); void mlib_VideoColorBlendABGR_ResetAlpha(mlib_u32 *dst, const mlib_u32 *src1,const mlib_u32 *src2, mlib_s32 src1_w, mlib_s32 src1_h, mlib_s32 src2_w,mlib_s32 src2_h, mlib_s32 src2_x, mlib_s32 src2_y, mlib_s32 dst_lb, mlib_s32 src1_lb, mlib_s32 src2_lb, mlib_blend src1_blend, mlib_blend src2_blend); DESCRIPTION
The functions use the following equation for blending images: dst = (src1 * src1_blend) + (src2 * src2_blend) The two multi-banded source images (src1 and src2) are blended together and stored in the destination image (dst). The image buffers pointed to by dst, src1, and src2 contain 4-banded ABGR images, 8 bits per component. src1_w and src1_h are the dimensions of the src1 input buffer. src2_w and src2_h are the dimensions of the src2 input buffer. The output buffer must be at least as large as the src1 input buffer. src2_x and src2_y are the offset of the src2 input buffer relative to src1. Where pixels in src2 overlap pixels in src1, the pixels are blended. Pixels in src1 which are outside of src2 are copied into dst. Pixels in the dst image outside of src1 are left unchanged. src1_blend specifies the blend function to be applied to the pixels of src1 image and src2_blend specifies the blend function to be applied to the pixels of src2. Possible blend functions are: MLIB_BLEND_ZERO MLIB_BLEND_ONE MLIB_BLEND_SRC_ALPHA MLIB_BLEND_ONE_MINUS_SRC_ALPHA MLIB_BLEND_DST_ALPHA MLIB_BLEND_ONE_MINUS_DST_ALPHA MLIB_BLEND_SRC_ALPHA is the alpha component of image src2 scaled to the range 0.0 to 1.0. MLIB_BLEND_DST_ALPHA is the alpha component of image src1 scaled to the range 0.0 to 1.0. All pixel components are treated as unsigned 8-bit quantities and the output pixel component values are clamped to the range 0 to 255. For the mlib_VideoColorBlendABGR_ResetAlpha() function, the alpha value of every pixel in destination image is set to 0 after blending is complete. PARAMETERS
Each of the functions takes the following arguments: dst Pointer to output image. src1 Pointer to 1st input image. src2 Pointer to 2nd input image. src1_w src1 image width in pixels. src1_h src1 image height in rows. src2_w src2 image width in pixels. src2_h src2 image height in rows. src2_x src2 horizontal displacement (in pixels), relative to the upper-left corner of src1. src2_y src2 vertical displacement (in rows), relative to the upper-left corner of src1. dst_lb Linebytes for output image. src1_lb Linebytes for 1st input image. src2_lb Linebytes for 2nd input image. src1_blend Blend function for src1 image. src2_blend Blend function for src2 image. RETURN VALUES
None. ATTRIBUTES
See attributes(5) for descriptions of the following attributes: +-----------------------------+-----------------------------+ | ATTRIBUTE TYPE | ATTRIBUTE VALUE | +-----------------------------+-----------------------------+ |Interface Stability |Committed | +-----------------------------+-----------------------------+ |MT-Level |MT-Safe | +-----------------------------+-----------------------------+ SEE ALSO
mlib_VideoColorBlendABGR_Inp(3MLIB), attributes(5) SunOS 5.11 2 Mar 2007 mlib_VideoColorBlendABGR(3MLIB)
All times are GMT -4. The time now is 10:13 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy