Hi, I am doing reciprocal blast to find ortholog. for that i run i create database and run blastn and retrieve hit sequence fasta through awk 'BEGIN { OFS = "\n" } { print ">"$2, $3 }' blast.out > hit.fa but when i am trying to align hit.fa to reference genome database, i am getting this issue.
I used this command.--
blastn -query ~/Blast/hit.fa -db ~/Blast/omykiss/genome/omykiss -outfmt 6 -out blast_out -num_alignments 1
please help to solve it.
Thanks in advance
•
written
by
baurumon •
1 answer
Hyphens (coding for indels relative to another sequence, from BLAST in this case) are not in {A, C, G, T, N}. You should remove them or replace them with Ns, depending on what you're trying to do.
A quick Google search reveals this was an issue in some other applications, as well.
Log in to answer this question.
-
Question on how to include TAXIDs when running blast on database
written by Mani •
-
Puzzling NCBI BLAST+ output after transcriptome assembly
written by ella •
-
Local blast with multiple databases and multiple files
written by akang
-
Blastn .XML to FASTA
written by baurumon •
-
Reciprocal blast commandline
written by baurumon •
-
blastn: how to show all query list in my blastn results, even if they have no match.
written by nuria.mauri •
-
BLAST Multiple Subject sequences?
written by zgayk •
-
Missing Blast fields in outfmt 6
written by zgayk •
-
parameter in blastn to limit the output to top 20 best hits
written by J.F.Jiang
-
data from SRA against blast databases method.
written by arronar
