Calculate pairwise distance for 3000 small sequences
I have a .fasta with 3000 ascensions, each of which is 10 bases long. I need to calculate the alignment score between each pair of sequences - that's 3000^2 comparisons. Using R or python, what's the best way to go about doing this? Technical details like code snippets are particularly helpful.
Thanks!!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
What exactly is Biopython's alignment.score variable?
written by Tom •
-
Pairwise alignment and Ka/Ks computation
written by maxime.policarpo
-
ML tree inferred from multiple genes
written by Lada •
-
Computing pairwise AA global sequence alignment between all pairs in a vector of sequences
written by rubic
-
weighted mean pairwise distance among samples containing multiple sequences?
written by Thomas Moore •
-
Pairwise alignment of multi-FASTA file sequences
written by aurora •
-
Adding sequences to preexisting MAF alignment
written by Carla •
-
Analyzing pair by pair DNA sequences with for loops in python : could it be faster?
written by nlehmann
-
Split a concatenated alignment in multiple files
written by dago
-
limit size of alignment generated by clustalw? Or related aligners?
written by user
This post has some tips and example python code: Massive Pairwise Comparison Using Biopython
This sounds like an interesting challenge for a 100 lines C program. ;)
You can use biostrings package with R. It is samplest methode I think.see