Entering edit mode
3.3 years ago
Aaron
•
0
I need to combine two separate reference genomes into one index, so that I can align my sequencing files to one assembly. Is there a way to do this in Hisat2? I need to be able to index them as one index and not multiple indices. Only used Hisat2 for aligning to single genomes, not experienced with creating custom genomes like this!
What did you try? Did you try making one large fasta with both genomes and indexing that?
I haven't, no. This is what I'm trying to do. Do I build the fasta using HIsat2? I would need to do the same for the annotations also.
No, just use cat to combine things, though you may have to see if your annotation files has headers. You will probably want to use sed or something to make sure the chromosome names are different between your two references, if they aren't already.