Hi! I'm trying to blast a sequence where I have a TSS signal (e.g. using CapSeq) that comes from the minus strand, against a chosen genome. I already have the reverse compliment of the sequence i.e. the reverse compliment of the position of the genome, so I'm trying to force BLASTN to only use the minus strand of the chosen subject. Sadly I only found the command -strand, which takes the rev. compliment of the query and not the subject.
Did I miss the parameter in the documentation or how do you do this?
Thanks in advance!
1 answer
After blastn, you can only look at the Plus/Minus strand results (that is your query vs. the reverse complement of the subject). Alternatively, you can do reverse complement for your subject sequence first, and then blastn and choose only the plus/plus strand results. You can use the default setting ("-strand both") or "-strand plus". It doesn't really matter in terms of the running speed using a local blast, I think.
Log in to answer this question.
-
How to determine the coordinate of the transposon insertion on the reverse reading?
written by Елизавета •
-
Reverse complement only sequences resulted from BLASTN minus
written by Oscar •
-
Distinguish between forward strand, plus strand, sense strand
written by tien •
-
Extracting Promoter Sequence from reverse Strand using NCBI Entrez and Eukaryotic Promoter Database?
written by Lys •
-
Best way for finding annotated counterparts of unknown transcript after BLAST
written by mwittep •
-
homer-annotation of the minus strand genes
written by c2c12 •
-
Megabast multistrand alignment
written by rangerunseen •
-
position specific sequence retrieval from whole genome sequence
written by psiwach29 •
-
Can I compare the e-value results across different Blast program?
written by qingxiangg •
-
Blastn(Strand= Plus/Minus and Plus/Plus) explication
written by midox
