pairwise comparison local BLAST
Hi,
I have to perform multiple pairwise blast comparisons, via the command line interface, in order to check how well overlap sequence match (typically 1000 bases at one end should match 1000 bases at the start of the other). Can I do this without having to compile a database for each piece?
• 4,164 views
•
link
updated
by
Ram
•
written
by
michiel.karrenbelt •
•
written
by
michiel.karrenbelt •
2 answers
Use the -subject option. But it will be faster if you generate database for each sequences.
Furthermore I recommend use as low word size (-word_size) as you can. Default word size of blastn (megablast) is 28 which may be too high for your purpose.
• 0 views
•
link
written
by
Csaba Kerepesi
Or maybe you can define -task dc-megablast if you want to find sequences with not a very high homology
• 0 views
•
link
written
by
Antonio R. Franco
Log in to answer this question.
More posts like this
-
Pathogenicity prediction experiment
written by Newbie •
-
Use of --detect_adapter_for_pe and order of adapter removal
written by SergeW •
-
Check for deletions at particular neighboring positions - Pysam pileup or fetch
written by z.klinovska •
-
How to check quality of multiple sequence alignment (MSA)?
written by Chirag Parsania
-
RNAseq differential expression analysis : learning design matrix
written by guillaume.rbt
-
Trying to reproduce mirbase results locally with BLAST
written by avino •
-
Dealing with CNV overlap
written by Max
-
A Tool For Flagging Blast Searches?
written by ishengomae
-
Visualizing Pairwise Distance Between A Group Of Sequences
written by David M
-
The "Read Bases" Column In Samtools Mpileup
written by dustar1986
Is there a reason you want to use blast for this?
I want to check the specificity of the overlaps