Biostar Beta. Not for public use.
Question: Long-read alignment with Minimap2
0
Entering edit mode

Hi all,

I want to align the coding sequence of a gene (which I cropped out from hg38) to genomic read contigs of a non-human ape. All contigs are in a single fasta file (2.8 GB) and the length of the coding sequence that I want to align to it is ~400 bp. I've been using Minimap2, here is my command:

minimap2 -ax map-ont ape_genome_contigs.fasta human_gene_coding_seq.fa > aln.sam

And here is the output of this command:

goekberk@bonobo:/disk1/goekberk$ minimap2/minimap2 -ax map-ont ape_genome_contigs.fasta human_gene_coding_seq.fa.fa > aln.sam
[M::main::11.782*1.00] loaded/built the index for 2771 target sequence(s)
[M::mm_mapopt_update::13.971*1.00] mid_occ = 630
[M::mm_idx_stat] kmer size: 15; skip: 10; is_hpc: 0; #seq: 2771
[M::mm_idx_stat::15.343*1.00] distinct minimizers: 99452876 (39.18% are singletons); average occurrences: 5.363; average spacing: 5.304
[M::main] Version: 2.17-r943-dirty
[M::main] CMD: minimap2/minimap2 -ax map-ont ape_genome_contigs.fasta human_gene_coding_seq.fa.fa
[M::main] Real time: 16.021 sec; CPU: 16.022 sec; Peak RSS: 7.437 GB

However, for some reason, sam files that I generate have 0 aligments.

$ samtools view -c aln.sam
0

I tried map-pb and asm10 options as well as minimap2 -a [-x preset] target.mmi query.fa > output.sam command but could not manage to have a successful alignment. I was wondering if I'm missing something crucial about the long-read alignments to a set of genomic contigs. Any help is much appreciated.

Cheers, Gökberk

ADD COMMENTlink 9 months ago gokberk • 30
Entering edit mode
2

Have you tried using blast+?

ADD REPLYlink 9 months ago
genomax
68k
Entering edit mode
0

I tried that but since the genome is about 2.8 GB, I could not manage to upload it and received errors even though I tried multiple times.

ADD REPLYlink 9 months ago
gokberk
• 30
Entering edit mode
1

You should try this with a local install of blast+. Even blat may be worth a try. No indexing needed.

ADD REPLYlink 9 months ago
genomax
68k
2
Entering edit mode

I guess my query sequence was not formatted as a proper fasta file, I fixed it and problem is solved now, sorry for bothering :).

ADD COMMENTlink 9 months ago gokberk • 30

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0