Hi I am trying to BLAST one large genome file with a specific set of sequences and I cannot figure out what is wrong with my script (I am obviously new to scripting). Any help would be greatly appreciated.
This is my script:
-in P.reformat.CDS.transcripts.fa -out p.full dbtype=nucl
blastn -query=NGF.fasta -db=p.full -outfmt="6" -out=blast
I want to search the "P.reformat.CDS.transcripts.fa" file for anything similar to the sequences in the "NGF.fasta" file.
This is my script:
-in P.reformat.CDS.transcripts.fa -out p.full dbtype=nucl
blastn -query=NGF.fasta -db=p.full -outfmt="6" -out=blast
I want to search the "P.reformat.CDS.transcripts.fa" file for anything similar to the sequences in the "NGF.fasta" file.
Comment