minimap2 setting for RNAseq of 150PE
Hello I need to map Illumina's 150PE RNAseq reads to a reference genome. I wonder if minimap is suitable for this task or I should use different software. The goal of the mapping is variant calling. If minimap is a good choice, which parameters should I use? Thank you
1 answer
Minimap2 is commonly used for long reads (Nanopore and company...) while for Illumina short reads one typically uses tools such as most commonly Hisat2 or STAR. They obviously need to be splice-aware which these two are (so is minimap2 but as said it is usually used for long reads). Any of the two should be a good choice.
Log in to answer this question.
More posts like this
-
Transcriptome vs genome reference in gene expression analysis
written by bioinfo1994 •
-
Alignment of RNASeq with STAR for variant calling
written by Gilles •
-
quantifying splicing transcripts for nanopore RNAseq
written by jshen •
-
RNASeq analysis without reference genome using Trinity
written by Dhananjai •
-
Best practices for variant calling on multiple sequencing runs of the same sample
written by Damian Kao
-
How to add sequence information to vcf from MPileup
written by helge.fluess •
-
choice of human genome before RNASeq mapping
written by Omics data mining
-
What are the applications of read-reference sequence mapping?
written by manetsus •
-
what is the difference between reference genome sequence from UCSC and ensembl ?
written by Pei
-
Map quality difference in Varscan somatic
written by andxiwang •
