Here's what I'm trying to do-
I'm trying to assemble a couple a candidate genes from a pool of illumina reads. I have a species close to drosophila and would like to align 60 million reads to a transcript from flybase and pull only those reads out and assemble them.
If I use standalone BLAST I can get csv output but the program only returns aligned portions and I want the whole read. Is there a way to format the output or another way to get the sequence if all I have is a list of IDs?
I'm trying to assemble a couple a candidate genes from a pool of illumina reads. I have a species close to drosophila and would like to align 60 million reads to a transcript from flybase and pull only those reads out and assemble them.
If I use standalone BLAST I can get csv output but the program only returns aligned portions and I want the whole read. Is there a way to format the output or another way to get the sequence if all I have is a list of IDs?
Comment