match, mismatch, gap parameters for short DNA motifs
Hello all,
I have millions of short DNA sequences (20-150 nt) and I suspect some of these show similarity. I would like to compare them using pairwise alignment method, however, I am having trouble in determining which match score, mismatch penalty and gap penalty (linear method) I should use. My research didnt result any guideline for such short sequences, so any help is appreciated.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Using sequence percent similarity from alignments to compare duplication events
written by Zeng Hao •
-
How to find if sequences perfectly overlap
written by gero.knittel •
-
Finding short, highly similar matches to long sequences using BLASTN
written by CC •
-
Pairwise alignment of multi-FASTA file sequences
written by aurora •
-
How to apply BioPython function to every row of CSV
written by ishackm
-
How to force blastn to align full subject
written by Bastien Hervé
-
multilabel clustering/classification based on sequence silimarity
written by odoluca •
-
Gap Continuation Penalty With Dynamic Programming ?
written by User 5037 •
-
How Do I Calculate The Overall Gap Penalty
written by User 5037 •
-
Similarity Score Of Multiple Sequence Alignment
written by Ananth •
Are these in fasta format? You may want to use CD-HIT to reduce redundancy instead of trying to do pair-wise alignments.
I checked the CD-HIT. Its webservers don't allow any file bigger than 50 MBs. So I have to compile myself. However, I already coded a pairwise comparison algorithm, so I want to stick with it, in which otherwise would take longer. Thanks anyways.