I have a metagenome like a db, and i need to blast a genome with this metagenome, the results are veryy wide, and i need to extract the aligment from this result and put into a fasta file, but i cant, because the result are in a *.html file. i do this..
formatdb
formatdb -i home/biolinux/Desktop/prueba/Metagenoma/l1l2/4499203.3.350.genecalling.coding.faa
Blastp
blastall -p blastp -d /home/biolinux/Desktop/prueba/Metagenoma/l1l2/4499203.3.350.genecalling.coding.faa -T T -i /home/biolinux/Desktop/prueba/Nitrincola.gbk -o /home/biolinux/Desktop/prueba/resulnitri.html
later i try to do a fastacmd, but i cant...
what i can do?