Post #5 is interesting.
I found it can be enhanced, by a ^ anchor (match at the beginning only), and by inserting a ^ character that becomes an anchor for the following grep
Post #7 did not work for me.? Did not further analyze it.
Here is an embedded multi-line awk code that works for me:
I chose $1 not $0 because I found some trailing spaces in the given examples; spaces are stripped if $1 is used.
Hi, I have two files where 1 contains data and the other contains strings eg
file 1
-0.00000 0.00000 0.00000
0.00000 0.00000 0.80000
0.50000 0.50000 0.60000
0.50000 0.50000 0.20000
-0.00000 0.00000 0.40000
file 2
F F F
F F F
T T T
T T T
T T T
How to I append file2 to file 1 to... (1 Reply)
I would like to extract the sequences larger than 10 bases but shorter than 18 along with the identifier from a FASTA file that looks like this:
> Seq I
ACGACTAGACGATAGACGATAGA
> Seq 2
ACGATGACGTAGCAGT
> Seq 3
ACGATACGAT
I know I can extract the IDs alone with the following code
grep... (3 Replies)
I have a fasta file that looks like this:
>Noname
ACCAAAATAATTCATGATATACTCAGATCCATCTGAGGGTTTCACCACTTGTAGAGCTAT
CAGAAGAATGTCAATCAACTGTCCGAGAAAAAAGAATCCCAGG
>Noname
ACTATAAACCCTATTTCTCTTTCTAAAAATTGAAATATTAAAGAAACTAGCACTAGCCTG
ACCTTTAGCCAGACTTCTCACTCTTAATGCTGCGGACAAACAGA
...
I want to... (2 Replies)
I tried to write a script ( not working) to append first value from mylist to a file called my myfirstResult and to another called mysecondResult
awk ' {print $1} >> myfirsResult ' < mylist
awk ' {print $1} >> mysecondResult ' < mylist
$ cat mylist
A 02/16/2012
B 02/19/2012
C... (3 Replies)
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)
Hi All,
I have to append 2 lines at the end of a text file. If those 2 lines are already there then do not append else append the 2 lines to the text file.
Eg: I have a text file, file.txt
This text file might look like this,
/home/kp/make.jsp
/home/pk/model.jsp
I have to append... (1 Reply)
Hi frnds,
My requirement is I have a zip file with name say eg: test_ABC_UH_ccde2a_awdeaea_20150422.zip
within that there are subdirectories on each directory we again have .zip files and in that we have files like mama20150422.gz and so on.
Iam in need of a bash script so that it unzips... (0 Replies)
Hii,
Could someone help me to append string to the starting of all the filenames inside a directory but it should exclude .zip files and subdirectories.
Eg.
file1: test1.log
file2: test2.log
file3 test.zip
After running the script
file1: string_test1.log
file2: string_test2.log
file3:... (4 Replies)
Bio::AlignIO::fasta(3pm) User Contributed Perl Documentation Bio::AlignIO::fasta(3pm)NAME
Bio::AlignIO::fasta - fasta MSA Sequence input/output stream
SYNOPSIS
Do not use this module directly. Use it via the Bio::AlignIO class.
DESCRIPTION
This object can transform Bio::SimpleAlign objects to and from fasta flat file databases. This is for the fasta alignment format, not for
the FastA sequence analysis program. To process the alignments from FastA (FastX, FastN, FastP, tFastA, etc) use the Bio::SearchIO module.
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
Peter Schattner
APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
next_aln
Title : next_aln
Usage : $aln = $stream->next_aln
Function: returns the next alignment in the stream.
Returns : Bio::Align::AlignI object - returns 0 on end of file
or on error
Args : -width => optional argument to specify the width sequence
will be written (60 chars by default)
See Bio::Align::AlignI
write_aln
Title : write_aln
Usage : $stream->write_aln(@aln)
Function: writes the $aln object into the stream in fasta format
Returns : 1 for success and 0 for error
Args : L<Bio::Align::AlignI> object
See Bio::Align::AlignI
_get_len
Title : _get_len
Usage :
Function: determine number of alphabetic chars
Returns : integer
Args : sequence string
width
Title : width
Usage : $obj->width($newwidth)
$width = $obj->width;
Function: Get/set width of alignment
Returns : integer value of width
Args : on set, new value (a scalar or undef, optional)
perl v5.14.2 2012-03-02 Bio::AlignIO::fasta(3pm)