Analyzing fasta files: BioPython

[medinab]

Hello!!!

I have a file with a list of accession numbers for the organisms of my interest. I also have two other fast files, one with a list of large ribosomal subunits, and one with a list of small ribosomal subunits. What I am trying to do is to get both large and small ribosomal subunits for the organisms of my interest.

Basically, what I have to do is to compare three files and isolate duplicates in all three.

I was hoping you guys can help me figure out a script in BioPython that would help me get the information I want.

Thank you!!!

[TiborNagy]

Create a set from the accession numbers, read the two fasta file and store the headers as a set. Now you can use intersection() method to get common elements.

https://docs.python.org/2/library/sets.html