how to make short custom reference sequence to use in RNA-seq alignment
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
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Align unmapped reads to reference vector sequence
written by Tonkatsu •
-
RNA-seq alignment in UNIX
written by Hien •
-
filtering alignned seqeunces
written by priscilladraj •
-
Return number of mismatches for multiple sequences
written by genya35 •
-
HISAT2 index generation
written by deshpande.neha2 •
-
How to achieve the same alignment with HiSAT2 command line vs. HiSAT2 Galaxy
written by genya35 •
-
hisat2 command line fastq to bam
written by genya35 •
-
how to make a .bam file from fastQ for RNA-seq ion torrent with BBmap
written by genya35 •
-
How to align reads using bowtie against multiple customized reference files
written by Jackie
-
BWA mapping with selected mismatch allowed?
written by lim24m •
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.