Sequence similarity networks with python
1
0
Entering edit mode
8 weeks ago
Eduardo • 0

How to build sequence similarity networks within python? Pythoscape seems no longer available and biopython pairwise does not seem the most suitable for thousands of sequences. Any ideas? Thank you

sequence networks similarity • 271 views
ADD COMMENT
1
Entering edit mode
8 weeks ago
Carambakaracho ★ 2.9k

There are lots of algorithms out there for aligning sequences, both pairwise alignments and multiple sequence alignments. These calculations are relatively slow, and you generally wouldn’t want to write such an algorithm in Python. Instead, you can use Biopython to invoke a command line tool on your behalf.

From Biopython-readthedocks.

ADD COMMENT
0
Entering edit mode

Thank you, somehow I didn't find that page before or didn't understand how useful it is to my problem.

Nevertheless, in the meantime, I did a blastp pairwise alignment and I leave the steps here as it might be useful to someone. It is a half baked solution as it seems to duplicate the workload (the same sequence pairs are calculated twice) and it automatically chooses only sequences that produce significant alignments (not all sequences you provide)

1 - prepare a database of all desired sequences

makeblastdb -in All_seqs.fasta -dbtype prot -out my_prot_blast_db

2 - pairwise alignment

blastp -db my_prot_blast_db -query All_seqs.fasta -out all-vs-all.tsv -num_threads 8

I will now try Biopython’s command line wrapper for EMBOSS

ADD REPLY

Login before adding your answer.

Traffic: 1132 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6