Sponsored Content
Top Forums UNIX for Advanced & Expert Users Map snps into a ref gene file Post 302989924 by marwah on Thursday 19th of January 2017 09:06:21 PM
Old 01-19-2017
I'm expecting that a gene might have more than snpID,
and there might be genes that don't have snps it will be NA
and there might be one snpID for pre one gene

Code:
awk 'FNR==1 {next} FILENAME=="pre_snpinfo_tumor.txt" {k++; POS[k]=$2; ID[k]=$2;} \  
                   FILENAME=="refFlat.txt" {i++; \
                                     if(POS[i]>=$5 && POS[i]<=$6) \
                                          print $1, ID[i], $3, POS[i], $5, $6} \
    ' pre_snpinfo_tumor.txt  refFlat.txt

but there is an error can you help please
 

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
Bio::Map::Prediction(3pm)				User Contributed Perl Documentation				 Bio::Map::Prediction(3pm)

NAME
Bio::Map::Prediction - An object representing the predictions of something that can have multiple locations in several maps. SYNOPSIS
use Bio::Map::Prediction; use Bio::Map::Position; # normally you would get predictions from a run wrapper like # Bio::Tools::Run::Meme, but here we create some manually: my $pred1 = Bio::Map::Prediction->new(-source => 'meme'); Bio::Map::Position->new(-element => $prediction1, -map => Bio::Map::GeneMap->get(-gene => 'gene1', -species => 'species1'), -start => 950, -end => 960); Bio::Map::Position->new(-element => $prediction1, -map => Bio::Map::GeneMap->get(-gene => 'gene1', -species => 'species2'), -start => 1950, -end => 1960); Bio::Map::Position->new(-element => $prediction1, -map => Bio::Map::GeneMap->get(-gene => 'gene2', -species => 'species1'), -start => 955, -end => 965); Bio::Map::Position->new(-element => $prediction1, -map => Bio::Map::GeneMap->get(-gene => 'gene2', -species => 'species2'), -start => 1955, -end => 1965); my $pred2 = Bio::Map::Prediction->new(-source => 'gerp'); Bio::Map::Position->new(-element => $prediction2, -map => Bio::Map::GeneMap->get(-gene => 'gene1', -species => 'species1'), -start => 950, -end => 960); # etc. # find the places where predictions agree use Bio::Map::GeneRelative; my $rel = Bio::Map::GeneRelative->new(-gene => 0); my $di = Bio::Map::Mappable->disconnected_intersections([$pred1, $pred2], -min_mappables_percent => 100, -min_map_percent => 100, -relative => $rel); my @positions = $di->get_positions; DESCRIPTION
For example, used to model transcription factor binding site predictions, which can have multiple locations in several maps. FEEDBACK
Mailing Lists User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the Bioperl mailing list. Your participation is much appreciated. bioperl-l@bioperl.org - General discussion http://bioperl.org/wiki/Mailing_lists - About the mailing lists Support Please direct usage questions or support issues to the mailing list: bioperl-l@bioperl.org rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible. Reporting Bugs Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the web: https://redmine.open-bio.org/projects/bioperl/ AUTHOR - Sendu Bala Email bix@sendu.me.uk APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ new Title : new Usage : my $prediction = Bio::Map::Prediction->new(); Function: Builds a new Bio::Map::Prediction object Returns : Bio::Map::Prediction Args : -name => string : name of the mappable element -id => string : id of the mappable element -source => string : name of the prediction program source Title : name Usage : $mappable->name($new_name); my $name = $mappable->name(); Function: Get/Set the name for this Mappable Returns : A scalar representing the current name of this Mappable Args : none to get string to set perl v5.14.2 2012-03-02 Bio::Map::Prediction(3pm)
All times are GMT -4. The time now is 05:04 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy