I have three fasta files and basically all i want is to get is the consensus or union of three fasta files. I mean, i want to pull out all the sequences that are common among three files and put them in a separate fasta file. The three fasta files were generated differently and have different headers but the sequences are identical among three files. Is there a tool or script that does this?
•
written
by
upendrakumar.devisetty
2 answers
Sorry if i didn't explained my question well. Here is what i wanted...
I have three fasta files
fasta1
seq1 ATGATG seq2 GATAGATA seq3 TGGTGG
fasta2
m1 GATAGATA m2 TGGTGG m3 AGGAGG
fasta3
seq1 gene1 ATGATG seq3 gene3 TGGTGG seq4 gene4 AGTGTG
And what i wanted is basically a final fasta file containing the below
final_fasta
seq3 TGGTGG
As you can see this sequence is represented in all fasta files but the header i got is from fasa1
I am currently using a lengthy process to achieve this.
First blast fasta1 to fasta2 and take all the hits and use those hits to blasta fasta3.
Thought it should work it is lengthy and i though if somebody has done something like this before i can use their script/tools
use CDHIT: http://weizhong-lab.ucsd.edu/cd-hit/
Log in to answer this question.
-
Generating a consensus sequence of a large number of amino acid sequences in python
written by unicorn2kang •
-
Compare two protein FASTA files and give a excel that show header with the same sequence
written by robert •
-
Reads Extraction
written by Amit •
-
Choosing random set of 500 fasta files from a folder containing large amount of fasta files in pyth…
written by Fadhalsh441 •
-
Parsing header of FASTA File
written by Shelle •
-
Sequence separation from multi fasta file
written by chetram.nbpgr •
-
Venn Diagram From Cap3 Output
written by upendrakumar.devisetty
-
How To Extract Full Sequence Of A Gene From Bam Or Sam File
written by upendrakumar.devisetty
-
Converting Consensus Fastq To Fasta
written by Love •
-
How To Compare Gene Models
written by Gvj
How large are the fasta files?
The files are not more than 100MB
If there is no existing tool/script, this is a fairly basic bioinformatics programming task which you should learn how to do.
I do understand it is basic bioinformatics programming but i thought if there is something already written then i can just use it rather than trying to write by myself
search for "multiple sequence alignment" in a search engine of your choice or encyclopedia to get a starting point.
Your question title and use of the term "consensus sequence" is misleading, it is reserved for alignments; it seems that you are looking to find exactly identical sequences instead. Please take at least one quarter of the time it will take us to answer your question, to put a proper example, and to explain what you already have tried.