Entering edit mode
5.4 years ago
genya35 ▴ 40
I would like to align reads from RNA fastq to a custom short reference sequence. I made a small fasta file with >Name in the header followed by DNA sequence. Do I need to do anything else to the file for tools lie BBmap or HISAT2 to be able to use it as a reference?
Try it, see what happens.
it didn't work with BBmap but I will try a different aligner, perhaps?
Post your command, and the error (if any), and perhaps an example of your FASTA.
I got it to work with HiSAT2 and it's awesome!
Glad to hear you got it working.