mkatari-bioinformatics-august-2013-blastnotes
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
mkatari-bioinformatics-august-2013-blastnotes [2013/08/16 18:14] – mkatari | mkatari-bioinformatics-august-2013-blastnotes [2014/07/02 13:58] – mkatari | ||
---|---|---|---|
Line 1: | Line 1: | ||
[[mkatari-bioinformatics-august-2013|Back to Manny' | [[mkatari-bioinformatics-august-2013|Back to Manny' | ||
+ | |||
+ | ====== Getting specific fasta sequences from reference file ====== | ||
+ | |||
+ | In case you need to retrieve a specific sequence from a larger fasta file, use one of my perl scripts. Simply provide a pattern that it needs to match in the definition field and it will retrieve the sequence. The **-f** is the reference file, **-o** is the output file name, and **-query** is the pattern it will try to match. | ||
+ | |||
+ | < | ||
+ | perl / | ||
+ | -o scaffold12498.fa \ | ||
+ | -query scaffold12498 | ||
+ | </ | ||
+ | |||
+ | ====== Creating the blast database ====== | ||
+ | |||
+ | Before we can actually perform the blast, we need to prepare the database using **makeblastdb**. | ||
+ | |||
+ | < | ||
+ | makeblastdb -in cassavaV5 -input_type " | ||
+ | </ | ||
====== Run Blast using sbatch ====== | ====== Run Blast using sbatch ====== | ||
Line 96: | Line 114: | ||
sbatch / | sbatch / | ||
</ | </ | ||
- |
mkatari-bioinformatics-august-2013-blastnotes.txt · Last modified: 2015/06/04 12:38 by mkatari