Question: how to make short custom reference sequence to use in RNA-seq alignment
0
gravatar for genya35
23 months ago by
genya3520
genya3520 wrote:

Hello,

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?

Thanks

reference sequence • 808 views
ADD COMMENTlink written 23 months ago by genya3520

Try it, see what happens.

ADD REPLYlink written 23 months ago by st.ph.n2.5k

it didn't work with BBmap but I will try a different aligner, perhaps?

ADD REPLYlink written 23 months ago by genya3520

Post your command, and the error (if any), and perhaps an example of your FASTA.

ADD REPLYlink written 23 months ago by st.ph.n2.5k

I got it to work with HiSAT2 and it's awesome!

ADD REPLYlink written 23 months ago by genya3520

Glad to hear you got it working.

ADD REPLYlink written 23 months ago by st.ph.n2.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1752 users visited in the last hour