Sponsored Content
Top Forums UNIX for Dummies Questions & Answers Probably some stupid mistake... Post 302491722 by sluvah on Friday 28th of January 2011 06:09:37 AM
Old 01-28-2011
[solved]Probably some stupid mistake...

Hi everyone !

I have a file wich look like this :

>Sis01
[lots of ATCG stuff]

> Sis02
[lots of ATCG stuff]

...

>Sis44
[lots of ATCG stuff]


I want to separe each paragraphe in a different file, so I decide to use the "FOR" loop + sed.

Code:
for f in {01..44}
do
(
   sed -n '/^>Sis `echo $f'/,/^>/p' Sis_110125.fasta | sed '$d' > Sis_110125_contig`echo $f`
)
done

I get empty files, i think due to the space between ^> and the echoed $f


I try a lot of thing, but I was unable to find a solution !

Thanks for your help !

Last edited by sluvah; 01-28-2011 at 02:10 PM.. Reason: Solved
 

8 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

Crontab Mistake!!!

Hi. I hope someone can help me with this problem. Being a novice to Unix, I editted my crontab directly by typing " crontab -e ". Well, I needed to make some changes so, I typed " crontab -r ". Now I have no crontab, and I can't seem to get crontab to write a new file. I' ve tried: vi... (4 Replies)
Discussion started by: cstovall
4 Replies

2. AIX

Did a Mistake with HACMP

Hi, I needed space on a FS, and when I've added the space on the filesystem, I did it trough the regular smitty fs inteface and not with smitty cl_lvm. Can someone help me to repair the situat before a faileover happen ? Thanks for your help,:mad: (13 Replies)
Discussion started by: azzed27
13 Replies

3. Shell Programming and Scripting

Is there any mistake in this code:

cat $1 | sort -n | uniq | $1 in other words, I sort the content of the file and put the ouput in the same file, is there any mistakes in this cshell code ??? (4 Replies)
Discussion started by: Takeeshe
4 Replies

4. Solaris

Renamed lib directory by mistake

Let's say someone accidentally renamed the lib directory in Solaris 8, and now they cannot get into the terminal or even rename the folder via file manager.What would one do? (37 Replies)
Discussion started by: jetjaguar
37 Replies

5. Red Hat

Changing Desktop Environment by mistake

Hi, I'm using fedora 15, my defualt DE is XFCE I once saw that there is a option in startup DE's that I can select and it was OPENBOX. I just wanted to test it. but after choosing it as my DE here I am. I have just access to terminal, firefox.. what's in my desktop. there is no panel here. I... (0 Replies)
Discussion started by: hoseinit
0 Replies

6. Shell Programming and Scripting

Can anyone find the mistake in this script file

#!/bin/ksh db_user=`echo $DB_USER_NAME` db_pwd=`echo $DB_PASSWORD` db_sid=`echo $TWO_TASK` sqlplus -s $db_user/$db_pwd@$db_sid << EOF a = select ACK_PARTY_NAME,bus_event_seq_nbr from bus_event where ack_party_name like 'MOVE_USAGE_DAEMON%' and bus_event_seq_nbr='3969094' set -- echo $a |... (17 Replies)
Discussion started by: rkrish
17 Replies

7. UNIX for Dummies Questions & Answers

Can anyone help me to spot my mistake?

Hi there can anyone help me to spot my mistake and please explain why it appears My code : #!/usr/bin/gawk -f BEGIN { bytes =0} { temp=$(grep "datafeed\.php" | cut -d" " -f8) bytes += temp} END { printf "Number of bytes: %d\n", bytes } when I am running ./q411 an411 an411: ... (6 Replies)
Discussion started by: FUTURE_EINSTEIN
6 Replies

8. Shell Programming and Scripting

Strange type mistake?!

Hi, I want to start MY_PROGRAM in a bash script with additional parameters given in the CONFIGURATION_ARRAY. IFS="'" CONFIGURATION_ARRAY=( '-N 0 -m 0' '-N 0 -m 1' ) for configuration in ${CONFIGURATION_ARRAY} do //DEBUG N=${configuration%-*} //-N 0 M=-${configuration##*-} //-m 0... (5 Replies)
Discussion started by: xraystorm
5 Replies
Bio::SeqIO::tab(3pm)					User Contributed Perl Documentation				      Bio::SeqIO::tab(3pm)

NAME
Bio::SeqIO::tab - nearly raw sequence file input/output stream. Reads/writes id" "sequence" " SYNOPSIS
Do not use this module directly. Use it via the Bio::SeqIO class. DESCRIPTION
This object can transform Bio::Seq objects to and from tabbed flat file databases. It is very useful when doing large scale stuff using the Unix command line utilities (grep, sort, awk, sed, split, you name it). Imagine that you have a format converter 'seqconvert' along the following lines: my $in = Bio::SeqIO->newFh(-fh => *STDIN , '-format' => $from); my $out = Bio::SeqIO->newFh(-fh=> *STDOUT, '-format' => $to); print $out $_ while <$in>; then you can very easily filter sequence files for duplicates as: $ seqconvert < foo.fa -from fasta -to tab | sort -u | seqconvert -from tab -to fasta > foo-unique.fa Or grep [-v] for certain sequences with: $ seqconvert < foo.fa -from fasta -to tab | grep -v '^S[a-z]*control' | seqconvert -from tab -to fasta > foo-without-controls.fa Or chop up a huge file with sequences into smaller chunks with: $ seqconvert < all.fa -from fasta -to tab | split -l 10 - chunk- $ for i in chunk-*; do seqconvert -from tab -to fasta < $i > $i.fa; done # (this creates files chunk-aa.fa, chunk-ab.fa, ..., each containing 10 # sequences) 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 one of the Bioperl mailing lists. 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 the bugs and their resolution. Bug reports can be submitted via the web: https://redmine.open-bio.org/projects/bioperl/ AUTHORS
Philip Lijnzaad, p.lijnzaad@med.uu.nl APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ next_seq Title : next_seq Usage : $seq = $stream->next_seq() Function: returns the next sequence in the stream Returns : Bio::Seq object Args : write_seq Title : write_seq Usage : $stream->write_seq($seq) Function: writes the $seq object into the stream Returns : 1 for success and 0 for error Args : Bio::Seq object perl v5.14.2 2012-03-02 Bio::SeqIO::tab(3pm)
All times are GMT -4. The time now is 08:25 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy