I'm having trouble aligning 3 sequences together - the total gene (~2500bp) and the forward and reverse sequences that are ~1000bp in length. Now I'm using Clustal X2, but the problem is that the program does complete alignment, which is useless for me right now. Could anyone recommend or suggest any parameters or methods on how to align more that 2 sequences together, but that each sequence is aligned to the first one (the total gene)? Thanks!
Please clarify your specific problem or add additional details to highlight exactly what you need. As it's currently written, it’s hard to tell exactly what you're asking. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.
Your question is not very clear, but I think what you are trying to do is align sequencing results to a reference gene sequence.
Your total gene sequence spans ~2500 bp. It sounds like you have a forward sequence at the start of gene, reading only about 1000 bp and the reverse strand reading in 1000 bp from the end (if this is Sanger sequencing, you likely only have 800-900 bp of accurate base reads by the way). Therefore, you have a gap of about 500 bp in the middle of your gene of interest.
There are (at least) two ways to do the alignment, and I would suggest the manual method if you are trying to verify gene sequencing.
In a word processing document, paste in your gene-of-interest (or the recombinant gene you are trying to make), as well as the base reads from the forward strand, and the reverse complementof your reverse strand. Make sure you remove any whitespace like line breaks and spaces as this will disrupt the rest of the process.
Manually select a small portion of your forward read (20-30 bp for example) and search for the same "phrase" in your gene-of-interest. If you find it, highlight it in both the sequence it came from and the gene-of-interest so that you know you have already matched them.
Carry on with adjacent bases for the forward strand. Repeat for the reverse complement of the reverse strand.
Go to NCBI Blast, and under the heading Basic Blast, click "nucleotide blast".
Check "Align 2 or more sequences."
Enter your Query Sequence as your 1000 bp sequences, entered in FASTA format. FASTA format will allow you to enter multiple sequences, but if you just enter one at a time, that will work perfectly well. Results will only display one at a time anyway.
> sequence name here ATCGATCGATCGATCGATCGATCG....
Enter your gene of interest as your Subject Sequence. Format does not matter here, but for consistency, I use FASTA format.
Keep "Highly similar sequences (megablast)" selected. You shouldn't need to play with these parameters, they just change how the alignment method penalizes alignments for gaps and mismatches.
Interpret your results. I've included an example of BLAST results below.
At the top of this window is a schematic, showing exactly how much of the the query sequence aligns with the subject sequence. The longer the line, the more bases align. Similarly, the closer the colour of the line is to right of the scale (e.g., red), then the more continuous the segment is indicating a longer, continuous alignment.
Below the graphic is a table of results. Briefly, the table highlights all possible alignments that the BLAST algorithm has found, including how much coverage there was, confidence of the alignment and length of alignment coverage.
Below the results table, are the detailed results which were summarized in the table. This detailed report shows you exactly which bases align or mismatch, or any gaps. The top line is the subject sequence; the bottom is the query sequence. The other thing to keep in mind here is the Strand line. In this case, it says "Plus/Plus, meaning that the aligned sequences were the forward sense, or plus strand, of each of subject and query sequence, respectively. If you entered a reverse complement sequence instead (such as the reverse sequencing results) then this would change to "Plus/Minus".
From here, you will still need to manually line these results up to your gene of interest to determine exactly which portion of the gene is missing sequence coverage, as well as any mutations that may have cropped up in the cloning process.
Multiple sequence alignment is a common bioinformatics problem and there are many tools available to solve it.
This paper provides a good comparison of the following tools:
The paper addresses the following methods:
If you want to align sequencing results over the reference sequence and check the sequences for presence of spontaneous mutations, the easiest way would be to use Sequlator.
Sequlator is a convenient because it has a pairwise alignment, reverse complement, editing and navigation functions but lacks irrelevant functions. With Sequlator one can manipulate any way existing alignments, or manually assemble sequencing reads into an alignment with the reference sequence. After elimination of poorly sequenced fragments, discrepancies with the reference sequence are revealed in the consensus of the resulting alignment. And it is free and Java-based meaning that it works on Windows, Linux and Mac computers