Hello I struggled with alignment step in rna-seq pipeline (my data is whole transcriptome of prokaryote by using direct RNA sequencing method with nanopore). The manuals and tutorials of alignment to reference transcriptome that I found were short read only. I have tried using minimap2 for alignment the RNA from nanopore FASTQ file to reference genome (FASTA) and gene annotation file (GTF) (this tools need convert the GTF,GFF3 to BED first). However, the SAM result I got is exactly the same to alignment to reference genome (FASTA file only). So my question is what is difference of output that they generate in SAM file between mapping with mapping to reference genome (FASTA) only and mapping with reference genome (FASTA) with gene annotation file (GTF or GFF). Thanks in advance!
•
written
by
sonsunjirachote •
Thank you. Do you have any recommend tools for long read alignment that generate output file as GTF/GFF file? I'm not sure that hisat2 and stringtie2 are work for long read data. Thanks in advance :)
Do you have any recommend tools for long read alignment that generate output file as GTF/GFF file
Output of an alignment is a SAM/BAM format file and not GTF/GFF. GTF/GFF files contain known information about gene models that is used along with a SAM/BAM alignment file to do read counting.
Log in to answer this question.
can you share the command for alignment with minimap? by default the output is a SAM file
1.FASTA+GTF
paftools.js gff2bed [my.gtf] > [my.bed]
minimap2 -a -x map-ont --junc-bed [my.bed] [my.fasta] [my.fastq] > [my.sam]
2.FASTA only
minimap2 -a -x map-ont [my.fasta] [my.fastq] > [my.sam]
There is a published pipeline from ONT that you should consider: https://github.com/nanoporetech/pipeline-transcriptome-de
Leaving this link here as a reference for future visitors.