Daniel:Notebook/HiResChrPaint/2014-6-19
Jump to navigation
Jump to search
Probe Design (Started 06/17/2014)
Probe Mining
The OligoArray script might actually need the blast database broken up into 1kb chunks for some reason. So, let's try that.
- Create BLAST database from fragmented genome
- formatdb -i in_hg38.fas -p T -o F
[formatdb] WARNING: Cannot add sequence number 1953021 (lcl|range=hg38:1953020000-1953020999) because it has zero-length. [formatdb] FATAL ERROR: Fatal error when adding sequence to BLAST database.
So that's a fatal error.
The instructions call for concatenating all the fasta files into one. I did this step using cat, but they suggest using their java script, so I'll try that.
- Concatenate fasta files
- java ConcatenateFiles hg38_cat.fas *.fa