Hello, I could use help with the following problem -
I've assembled a transcriptome and I can BLAST against Drosophila to identify some of the transcripts. Now I have two files: the BLAST output and transcript fasta. I would like to append the names (subject IDs) from the BLAST file to the fasta >headers.
The closest I can get is to format my BLAST output to a CSV with the query and subject IDs in the first two columns, but I'm not sure how to put together a script that will search for transcripts in my fasta and add the IDs.
Any advice or other solutions?
Thanks
I've assembled a transcriptome and I can BLAST against Drosophila to identify some of the transcripts. Now I have two files: the BLAST output and transcript fasta. I would like to append the names (subject IDs) from the BLAST file to the fasta >headers.
The closest I can get is to format my BLAST output to a CSV with the query and subject IDs in the first two columns, but I'm not sure how to put together a script that will search for transcripts in my fasta and add the IDs.
Any advice or other solutions?
Thanks
Comment