This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to find identity percentage of two genomes sequences stored in two variable by carrying using BLAST in biopython

Suppose I have the genomes of two organisms stored in variable 'a' and 'b'. I want to BLAST them against each other and store the identity percentage in another variable 'c'. I'm new to Biopython. Can anyone please provide me the code my which I can solve this problem.

alignment gene python biopython
• 1,357 views • link
updated by ATpoint ATpoint • written by prittam32 • prittam32

You won't get a single percentage identity with BLAST - what do you want to do with multiple numbers?. BioPython itself is also probably not the best way to go about this. Can you not make use of a commandline tool?

• 0 views • link
written by Joe Joe

Thank You for your reply. I have a list containing the genome sequences of 10 different organisms. I will make another list of antibiotic resistance gene sequences and check their presence in each genome using BLASTn. I tried pairwise2 but it goes out of memory. So i tried splitting genomic sequence into various overlapping sequence. Although it works, It is very time-consuming. So I wish to use BLASTn to reduce the time.

• 0 views • link
written by prittam32 • prittam32

It should be feasible to blast each ARG against each genome sequence, one genome at a time.

Please share the code you’re currently use?

• 0 views • link
written by Joe Joe

0 answers

No answers yet.

Log in to answer this question.

More posts like this