Align two sets of amino acid sequences

[zerhacker]

Hi all, I am a microbiology student with little knowledge to bioinformatics and programming outside of Illumina reads alignment and denovo assembly. Recently I was tasked to compare two closely related strains of the same species and to identify unique pathways amongst the two that allow each to specialize in their niche based on the complete sequences present on genbank prior to any wet lab procedures.

To approach this, I plan to extract all open reading frames from both strains, pull out shared/highly similar ORF and unique ORF, then find the pathway the unique ORFs are involved in to draw conclusion. I've extracted all open reading frames from both organisms using prodigal, based solely on the in frame non-interupted sequence between start and stop codon. Each of strain have about 2000 real and hypothetical AA sequences.

now I'm stuck trying to extract shared and unique sequences from both organisms. Are there any programs that is suitable for this task? All replies are appreciated!

[GenoMax]

If these two strains are relatively closely related then you can identify the similarities using BLAT (https://genome.ucsc.edu/FAQ/FAQblat.html). Post-alignment processing will have to be done to extract the information you need from the results.

You could learn to do some of this but if you are working against a deadline then it may be better to find a programmer friend or your local bioinformatics support facility. They should be able to this for you.

[GenoMax]

CD-HIT-2D may be useful: http://weizhong-lab.ucsd.edu/cdhit_s...?cmd=cd-hit-2d

Best of all you can try it yourself without waiting for someone's help. You may still need to do some parsing afterwards.

[zerhacker]

If these two strains are relatively closely related then you can identify the similarities using BLAT (https://genome.ucsc.edu/FAQ/FAQblat.html). Post-alignment processing will have to be done to extract the information you need from the results.

You could learn to do some of this but if you are working against a deadline then it may be better to find a programmer friend or your local bioinformatics support facility. They should be able to this for you.

thank you! I checkout out the programs that you suggested, but I ended up generating a fake sets of illumina reads out of both sequences using Simseq https://github.com/jstjohn/SimSe,
then I used bowtie2 to align them to each other and pulled out reads that dont align, then denovo assemble them into short contigs and extracted their ORF which codes for unique proteins.
I'm book marking BLAT as it seem like a fairly useful program.

Edited: bolded out my procedure to make it easier to read

[GenoMax]

Long as you were able to get what you needed :-)

What program did you use to generate the "illumina" reads. Just for the record. For someone running across this thread later-on via a search.

[zerhacker]

Long as you were able to get what you needed :-)

What program did you use to generate the "illumina" reads. Just for the record. For someone running across this thread later-on via a search.

https://github.com/jstjohn/SimSe

I think Simseq works great. but I used a python script wrote by the departments programmer that works similarly.