RNA SEQ reads assembly for illumina sequenced data
Hi everyone,
I have 10 paired end samples. I have done the quality check and trimming then I merged the trimmed reads and after that I have done the mapping using minimap2 using this command:
minimap2 -ax map-ont ref.fna(reference genome) trimmed.fq(merged trimmed reads) >aln.sam
I want to do the assembly now, but I don't find any tool where I can use this sam file
Please help me with this.
Thank you
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
HLA alignment, approximation and frequencies
written by oumo •
-
Read Classification from ONT
written by maevalefeuvre •
-
No Sequence Segments for secondary alignments in minimap2 while trying to run alignment on ONT data?
written by Adarsh •
-
How to check the version of PacBio HiFi genomic reads?
written by draccident •
-
Polishing PACBIO corrected assembly using ONT and ILLUMINA data (.fastq)
written by K •
-
filtering out nanopore reads by alignment score
written by resug •
-
Question regarding minimap2 presets
written by Rogerio Ribeiro •
-
Mapping to custom reference minimap2/IGV
written by sendhelp •
-
Long-read alignment with Minimap2
written by gokberk
-
aligner for 1D^2 oxford nanopore data
written by KVC_bioinfo
You have paired end reads from Nanopore?
sorry its from illumina
What's your end goal? You should use programs like STAR, Salmon, Kallisto for mapping/generating count matrices.
Then you probably shouldn't use the
map-ontpreset for Nanopore reads? And by merging, you mean you merged overlapping paired-end reads or you merged all 10 samples?If your goal is a reference-guided de-novo assembly of the transcriptome, have a look at StringTie. Here are two tutorials (1,2) or you can use a ready-made pipeline with this capability.
just curious why do you need read assembly, is it known or unknown organism? why minimap2 in case of RNA-seq ?