How to append two fasta files? Post: 303036022

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

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

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

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

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

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?

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

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

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

LEARN ABOUT DEBIAN

bp_mask_by_search

BP_MASK_BY_SEARCH(1p)					User Contributed Perl Documentation				     BP_MASK_BY_SEARCH(1p)

NAME

       mask_by_search - mask sequence(s) based on its alignment results

SYNOPSIS

	 mask_by_search.pl -f blast genomefile blastfile.bls > maskedgenome.fa

DESCRIPTION

       Mask sequence based on significant alignments of another sequence.  You need to provide the report file and the entire sequence data which
       you want to mask.  By default this will assume you have done a TBLASTN (or TFASTY) and try and mask the hit sequence assuming you've
       provided the sequence file for the hit database.  If you would like to do the reverse and mask the query sequence specify the -t/--type
       query flag.

       This is going to read in the whole sequence file into memory so for large genomes this may fall over.  I'm using DB_File to prevent keeping
       everything in memory, one solution is to split the genome into pieces (BEFORE you run the DB search though, you want to use the exact file
       you BLASTed with as input to this program).

       Below the double dash (--) options are of the form --format=fasta or --format fasta or you can just say -f fasta

       By -f/--format I mean either are acceptable options.  The =s or =n or =c specify these arguments expect a 'string'

       Options:
	   -f/--format=s    Search report format (fasta,blast,axt,hmmer,etc)
	   -sf/--sformat=s  Sequence format (fasta,genbank,embl,swissprot)
	   --hardmask	    (booelean) Hard mask the sequence
			    with the maskchar [default is lowercase mask]
	   --maskchar=c     Character to mask with [default is N], change
			    to 'X' for protein sequences
	   -e/--evalue=n    Evalue cutoff for HSPs and Hits, only
			    mask sequence if alignment has specified evalue
			    or better
	   -o/--out/
	   --outfile=file   Output file to save the masked sequence to.
	   -t/--type=s	    Alignment seq type you want to mask, the
			    'hit' or the 'query' sequence. [default is 'hit']
	   --minlen=n	    Minimum length of an HSP for it to be used
			    in masking [default 0]
	   -h/--help	    See this help information

AUTHOR - Jason Stajich
       Jason Stajich, jason-at-bioperl-dot-org.

perl v5.14.2							    2012-03-02						     BP_MASK_BY_SEARCH(1p)