Alignment of 2 fasta files already aligned
Hi guys,
I'm currently working on something, which is quite hard to me. Indeed, I have 2 fasta files already aligned, and I want to align the 2 fasta files which contain protein sequences already aligned on BioPython.
Thank's for your help.
sequence
alignment
• 2,129 views
•
link
updated
by
lieven.sterck
•
written
by
ptitfab94 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Sorting and writing multifasta entries to new fasta files
written by lachiemck •Hi, first post here. So I'm trying take the CDS out of various species' orthologous sequences. I'm running on a Linux server, and am mainly …
-
Fishing nucleotide sequences from genome files given aligned protein FASTA files using biopython
written by mirali1025 •Hello everyone: I have written a script using biopython to try and fish out nucleotide sequences of aligned proteins (within FASTA files). My task is …
-
Multiple Sequence Alignment Python
written by Manu Madhavan •Hi, Currently I am working on classifying coding/non coding RNAs based on sequence features. I want to try MSA score as a feature for sequence …
-
Replace the sequence
written by 22.venkat •Hi, I have two multiple sequence fasta files <p>file-1 I have 4000 fasta sequences with different IDs <p>file-2 I have 100 fasta sequences which IDs …
-
pairwise alignning of multiple protein sequences
written by vahideh •Hi! I have 4 protein sequences in FASTA format that its 2 sequences is for Human and others are for Mouse! I want to align …
-
Obtaining alignment from psiblast
written by nilus1432 •Hi everyone, I am using stand alone version of Psi-blast. I want to get just the query sequence followed by matched-hit sequences (both aligned) in …
-
how to compare protein sequence in Python
written by Jason Lin •Hi there, So I have lists which contain the PDBID for protein sequence and protein sequence such as: ``` >102L MNIFEMLRIDEGLRLKIYKDTEGYYTIGIGHLLTKSPSLNAAAKSELDKAIGRNTNGVITKDEAEKLFNQDVDAAVRGILRNAKLKPVYDSLDAVRRAALINMVFQMGETGVAGFTNSLRMLQQKRWDEAAVNLAKSRWYNQTPNRAKRVITTFRTGTWDAYKNL >103L MNIFEMLRIDEGLRLKIYKDTEGYYTIGIGHLLTKSPSLNSLDAAKSELDKAIGRNTNGVITKDEAEKLFNQDVDAAVRGILRNAKLKPVYDSLDAVRRAALINMVFQMGETGVAGFTNSLRMLQQKRWDEAAVNLAKSRWYNQTPNRAKRVITTFRTGTWDAYKNL >104L MNIFEMLRIDEGLRLKIYKDTEGYYTIGIGHLLTKSPSLNAAKSAAELDKAIGRNTNGVITKDEAEKLFNQDVDAAVRGILRNAKLKPVYDSLDAVRRAALINMVFQMGETGVAGFTNSLRMLQQKRWDEAAVNLAKSRWYNQTPNRAKRVITTFRTGTWDAYKNL …
-
blast output for FASTA (aligned sequences) to get organism name
written by s.m.karim •<p>Hi,</p> <p>I want to get the organism name from blast output (FASTA (aligned sequences)) for the specific protein.</p> <p>I am a beginner in python. Please …
-
How To Use Coordinates In Order To Extract Sequences In Fasta File?
written by srpatel144 •<p>I have obtained a .coords file from MUMmer after aligning my assembly to the reference genome.</p> <p>I have been asked to write a script that …
-
Sub-Record Indexing Of Fasta Files In Python?
written by interfect •<p>Hello,</p> <p>I have several genomes, each spread across multiple FASTA files. Some of these files contain multiple short sequences, and some of them contain large, …
How To Ask Good Questions On Technical And Scientific Forums
Apart from Pierre Lindenbaum comment, does it have to be with biopython? Several MSA software will do this as well (eg. muscle, t-coffee, ...)
Yes, it has to be down on BioPython
Do you want to dealign and realign all the sequences, or do a profile-profile alignment?
To realign what what previously aligned
Please make this easier on us. It's still not clear (maybe it's just me) what you want? As a general rule, if your question fits into a sentence or two, there's not enough information in it. Please at least provide some sample input data we can use, and ideally, some exemplar desired output.
Do you want to:
or
Yes sorry if i'm not clear, it is exactly like "Take all the sequences from both alignments, de-align them all to get back to 'normal' sequences, and then create one 'super alignment'?"
Ok, that's quite straight forward.
Can you explain your requirement for BioPython more? Do the alignments have to be done with pure python, or are you allowed to use, for example
MafftCommandlinefrom within BioPython?I don't know, we didn't have any indications about this. But I have never hear anything about the "MafftCommandline". So if you have both solution, with and without this this command, I would like to see !
is this a homework / course assignment?
It is a kind of "summer project", which is very long, and I am blocked to this specific part of the work.