Sponsored Content
Top Forums UNIX for Beginners Questions & Answers How to append two fasta files? Post 303036069 by Scrutinizer on Thursday 13th of June 2019 12:30:46 PM
Old 06-13-2019
awk version you could try:
Code:
awk 'BEGIN{FS=RS; RS=">"; ORS=""} FNR>1{A[$1]=RS $0} END{for(i in A) print A[i]}'  file1 file2

note: ensure that there are no excess trailing spaces as is the case with the file samples.

Last edited by Scrutinizer; 06-13-2019 at 01:50 PM..
These 3 Users Gave Thanks to Scrutinizer For This Post:
 

10 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

append two files

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)
Discussion started by: princessotes
1 Replies

2. UNIX for Dummies Questions & Answers

grep FASTA files

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)
Discussion started by: Xterra
3 Replies

3. UNIX for Dummies Questions & Answers

renaming (renumbering) fasta files

I have a fasta file that looks like this: >Noname ACCAAAATAATTCATGATATACTCAGATCCATCTGAGGGTTTCACCACTTGTAGAGCTAT CAGAAGAATGTCAATCAACTGTCCGAGAAAAAAGAATCCCAGG >Noname ACTATAAACCCTATTTCTCTTTCTAAAAATTGAAATATTAAAGAAACTAGCACTAGCCTG ACCTTTAGCCAGACTTCTCACTCTTAATGCTGCGGACAAACAGA ... I want to... (2 Replies)
Discussion started by: Oyster
2 Replies

4. Shell Programming and Scripting

append to two files

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)
Discussion started by: Sara_84
3 Replies

5. 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

6. UNIX for Dummies Questions & Answers

Append Files

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)
Discussion started by: pavan_test
1 Replies

7. UNIX for Dummies Questions & Answers

Append file name to fasta file headers in Linux

How do we append the file name to fasta file headers in multiple fasta-files in Linux? (10 Replies)
Discussion started by: Mauve
10 Replies

8. Shell Programming and Scripting

Unzip all the files with subdirectories present and append a part of string from the main .zip files

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)
Discussion started by: Ravi Kishore
0 Replies

9. Shell Programming and Scripting

Append string to all the files inside a directory excluding subdirectories and .zip files

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)
Discussion started by: Ravi Kishore
4 Replies

10. Shell Programming and Scripting

Help with reformat single-line multi-fasta into multi-line multi-fasta

Input File: >Seq1 ASDADAFASFASFADGSDGFSDFSDFSDFSDFSDFSDFSDFSDFSDFSDFSD >Seq2 SDASDAQEQWEQeqAdfaasd >Seq3 ASDSALGHIUDFJANCAGPATHLACJHPAUTYNJKG ...... Desired Output File >Seq1 ASDADAFASF ASFADGSDGF SDFSDFSDFS DFSDFSDFSD FSDFSDFSDF SD >Seq2 (4 Replies)
Discussion started by: patrick87
4 Replies
awk(1)							      General Commands Manual							    awk(1)

Name
       awk - pattern scanning and processing language

Syntax
       awk [-Fc] [-f prog] [-] [file...]

Description
       The  command scans each input file for lines that match any of a set of patterns specified in prog.  With each pattern in prog there can be
       an associated action that will be performed when a line of a file matches the pattern.  The set of patterns may appear literally  as  prog,
       or in a file specified as -f prog.

       Files  are  read  in  order;  if there are no files, the standard input is read.  The file name `-' means the standard input.  Each line is
       matched against the pattern portion of every pattern-action statement; the associated action is performed for each matched pattern.

       An input line is made up of fields separated by white space.  (This default can be changed by using FS, as described  below.)   The  fields
       are denoted $1, $2, ... ; $0 refers to the entire line.

       A pattern-action statement has the form

	    pattern { action }

       A missing { action } means print the line; a missing pattern always matches.

       An action is a sequence of statements.  A statement can be one of the following:

	    if ( conditional ) statement [ else statement ]
	    while ( conditional ) statement
	    for ( expression ; conditional ; expression ) statement
	    break
	    continue
	    { [ statement ] ... }
	    variable = expression
	    print [ expression-list ] [ >expression ]
	    printf format [ , expression-list ] [ >expression ]
	    next # skip remaining patterns on this input line
	    exit # skip the rest of the input

       Statements  are terminated by semicolons, new lines or right braces.  An empty expression-list stands for the whole line.  Expressions take
       on string or numeric values as appropriate, and are built using the operators +, -, *, /, %,  and concatenation	(indicated  by	a  blank).
       The  C operators ++, --, +=, -=, *=, /=, and %= are also available in expressions.  Variables may be scalars, array elements (denoted x[i])
       or fields.  Variables are initialized to the null string.  Array subscripts may be any string, not necessarily numeric; this allows  for  a
       form of associative memory.  String constants are quoted "...".

       The  print  statement prints its arguments on the standard output (or on a file if >file is present), separated by the current output field
       separator, and terminated by the output record separator.  The statement formats its expression list according to the format.  For  further
       information, see

       The  built-in  function	length	returns the length of its argument taken as a string, or of the whole line if no argument.  There are also
       built-in functions exp, log, sqrt, and int.  The last truncates its argument to an integer.  substr(s, m, n) returns the  n-character  sub-
       string  of  s that begins at position m.  The function sprintf(fmt, expr, expr, ...)  formats the expressions according to the format given
       by fmt and returns the resulting string.

       Patterns are arbitrary Boolean combinations (!, ||, &&, and parentheses)  of  regular  expressions  and	relational  expressions.   Regular
       expressions  must be surrounded by slashes and are as in egrep.	Isolated regular expressions in a pattern apply to the entire line.  Regu-
       lar expressions may also occur in relational expressions.

       A pattern may consist of two patterns separated by a comma; in this case, the action is performed for all lines between	an  occurrence	of
       the first pattern and the next occurrence of the second.

       A relational expression is one of the following:

	    expression matchop regular-expression
	    expression relop expression

       where a relop is any of the six relational operators in C, and a matchop is either ~ (for contains) or !~ (for does not contain).  A condi-
       tional is an arithmetic expression, a relational expression, or a Boolean combination of these.

       The special patterns BEGIN and END may be used to capture control before the first input line is read and after the last.   BEGIN  must	be
       the first pattern, END the last.

       A single character c may be used to separate the fields by starting the program with

	    BEGIN { FS = "c" }

       or by using the -Fc option.

       Other  variable	names  with special meanings include NF, the number of fields in the current record; NR, the ordinal number of the current
       record; FILENAME, the name of the current input file; OFS, the output field separator (default blank); ORS,  the  output  record  separator
       (default new line); and OFMT, the output format for numbers (default "%.6g").

Options
       -	 Used for standard input file.

       -Fc	 Sets interfield separator to named character.

       -fprog	 Uses prog file for patterns and actions.

Examples
       Print lines longer than 72 characters:
	    length > 72

       Print first two fields in opposite order:
	    { print $2, $1 }

       Add up first column, print sum and average:
		 { s += $1 }
	    END  { print "sum is", s, " average is", s/NR }

       Print fields in reverse order:
	    { for (i = NF; i > 0; --i) print $i }

       Print all lines between start/stop pairs:
	    /start/, /stop/

       Print all lines whose first field is different from previous one:
	    $1 != prev { print; prev = $1 }

Restrictions
       There  are  no explicit conversions between numbers and strings.  To force an expression to be treated as a number add 0 to it; to force it
       to be treated as a string concatenate "" to it.

See Also
       lex(1), sed(1)
       "Awk - A Pattern Scanning and Processing Language" ULTRIX Supplementary Documents Vol. II: Programmer

																	    awk(1)
All times are GMT -4. The time now is 01:15 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy