03-08-2013
Well, this gives me exactly all the entries common with original and allele freq file without the MAF values
chr1 979748
chr1 1120377
chr1 1178925
chr1 1222958
10 More Discussions You Might Find Interesting
1. UNIX for Dummies Questions & Answers
Hi, i was looking for unix command(s) for :
find the first occurrence of a given pattern with in a file and print the remaining part.
below is an example of what i am looking for:
lets say, a file named myfile.txt
now, the command i am looking for will do the following (4 Replies)
Discussion started by: nurulamin862
4 Replies
2. Shell Programming and Scripting
Hi,
I have two files. 1st file has 1 column (huge file containing ~19200000 lines) and 2nd file has 2 columns (small file containing ~6000 lines).
#################################
huge_file.txt
a
a
ab
b
##################################
small_file.txt
a 1.5
b 2.5
ab ... (4 Replies)
Discussion started by: AshwaniSharma09
4 Replies
3. Shell Programming and Scripting
Hello all, please help. There are two file like this:
file1:
1197510.0 294777.7 9666973.0 21.6 1839.8
1197510.0 294777.7 9666973.0 413.2 2075.9
1197510.0 294777.7 9666973.0 689.3 2260.0
... (1 Reply)
Discussion started by: attila
1 Replies
4. UNIX for Dummies Questions & Answers
Hello,
I am trying to modify 2 files, to yield results in a 3rd file.
File-1 is a 8-columned file, separted with tab.
1234:1 xyz1234 blah blah blah blah blah blah
1234:1 xyz1233 blah blah blah blah blah blah
1234:1 abc1234 blah blah blah blah blah blah
n/a RRR0000 blah blah blah... (1 Reply)
Discussion started by: ad23
1 Replies
5. Shell Programming and Scripting
I have a file with very specific column spacing formatting,
I wish to do the following:
awk '{print $1, $2, $3, $4, $5, $6, $19-$7, $20-$8, $21-$9, $10, $11, $12}' merge.pdb > vector.pdb
but the format gets ruined.
I have tried with print -f but to no avail.... (7 Replies)
Discussion started by: chrisjorg
7 Replies
6. Shell Programming and Scripting
Hi All,
I have to match each row in file 1 with 1st row in file 2 and print the corresponding column from file2. I am trying to use an awk script to do this. For example
cat File1
X1
X3
X4
cat File2
ID X1 X2 X3 X4
A 1 6 2 1
B 2 7 3 3
C 3 8 4 1
D 4 9 1 1 (3 Replies)
Discussion started by: newpro
3 Replies
7. Shell Programming and Scripting
Hello,
I have two files
File 1 with 10 columns
rsid position ........
xx 1:10000
File 2
position
1:10000
2:2000
....
I need to extract the IDs given in file 2(column1) from file 1 (column2) and print all columns from file1. I am trying this command (1 Reply)
Discussion started by: nans
1 Replies
8. UNIX for Beginners Questions & Answers
Hello, I have two tab files with headers
File1: with 4 columns
header1 header2 header3 header4
44 a bb 1
57 c ab 4
64 d d 5
File2: with 26 columns
header1.. header5 header6 header7 ... header 22...header26
id1 44 a bb
id2 57 ... (6 Replies)
Discussion started by: nans
6 Replies
9. UNIX for Beginners Questions & Answers
I have two text files. File 1 has 150 ids but all the ids exists in duplicates so it has 300 ids in total. File 2 has 1500 ids but all exists in duplicates so file 2 has 300 ids in total. i want to match the first occurance of every id in file 1 with first occurance of thet id in file 2 and 2nd... (2 Replies)
Discussion started by: limd
2 Replies
10. Shell Programming and Scripting
In the awk below I am trying to output those lines that Match between file1 and file2, those Missing in file1, and those missing in file2. Using each $1,$2,$4,$5 value as a key to match on, that is if those 4 fields are found in both files the match, but if those 4 fields are not found then missing... (0 Replies)
Discussion started by: cmccabe
0 Replies
LEARN ABOUT DEBIAN
bio::alignio::maf
Bio::AlignIO::maf(3pm) User Contributed Perl Documentation Bio::AlignIO::maf(3pm)
NAME
Bio::AlignIO::maf - Multiple Alignment Format sequence input stream
SYNOPSIS
Do not use this module directly. Use it via the Bio::AlignIO class.
use Bio::AlignIO;
my $alignio = Bio::AlignIO->new(-fh => *STDIN, -format => 'maf');
while(my $aln = $alignio->next_aln()){
my $match_line = $aln->match_line;
print $aln, "
";
print $aln->length, "
";
print $aln->num_residues, "
";
print $aln->is_flush, "
";
print $aln->num_sequences, "
";
$aln->splice_by_seq_pos(1);
print $aln->consensus_string(60), "
";
print $aln->get_seq_by_pos(1)->seq, "
";
print $aln->match_line(), "
";
print "
";
}
DESCRIPTION
This class constructs Bio::SimpleAlign objects from an MAF-format multiple alignment file.
Writing in MAF format is currently unimplemented.
Spec of MAF format is here:
http://genome.ucsc.edu/FAQ/FAQformat
FEEDBACK
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 the bugs and their resolution. Bug reports can be submitted via the
web:
https://redmine.open-bio.org/projects/bioperl/
AUTHORS - Allen Day
Email: allenday@ucla.edu
APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
new
Title : new
Usage : my $alignio = Bio::AlignIO->new(-format => 'maf'
-file => '>file',
-idlength => 10,
-idlinebreak => 1);
Function: Initialize a new L<Bio::AlignIO::maf> reader
Returns : L<Bio::AlignIO> object
Args :
next_aln
Title : next_aln
Usage : $aln = $stream->next_aln()
Function: returns the next alignment in the stream.
Throws an exception if trying to read in PHYLIP
sequential format.
Returns : L<Bio::SimpleAlign> object
Args :
perl v5.14.2 2012-03-02 Bio::AlignIO::maf(3pm)