minimap2, before or after assembly?
1
0
Entering edit mode
9 weeks ago

Dear all, I feel confused because I saw someone uses the minimap2 after demultiplexing, but before proceeding with the assembly (CANU) [case 1], and someone using the minimap2/samfile/BFCTools/medaka after assembly (always CANU) [case 2]. In case 1, the reference file used to align against the fastq files was a public available sequence from NCBI (NC_003310), while in the case 2 the file used as reference were the contigs from CANU to obtain the polished consensus. To be honest, I don't understand when is it convenient to use the 1st or the 2nd approach to get the best result. Does it only depend on having or not the public reference genome? Please, can someone give me some help? Thank you very much. Emilio

canu nanopore minimap2 • 550 views
ADD COMMENT
0
Entering edit mode

It entirely depends what you want to do.

ADD REPLY
0
Entering edit mode

minimap2 is a sequence alignment tool and is blind to whether or not you are aligning from or to the reference genome. In your examples - Case 1 is two genome sequences alignment, Case 2 is aligning contig sequences to a reference genome

ADD REPLY
0
Entering edit mode

Why don't you start by listing

  • your dataset
  • your research question
  • your goal (de novo assembly, find SNPs, something else)

Then people can help you more to get to your goal.

ADD REPLY
0
Entering edit mode

You are right. I am sorry to have been not so clear. Dataset: we have a collection of fast5 files from Nanopore MinION, in total 10 barcodes Research question: metagenomic analysis, in detail: identify viral sequences as much precise as possible My goal:to improve the assembled contigs obtained with MegaHit/CANU to get the whole genome

Thank you very much

ADD REPLY
0
Entering edit mode
8 weeks ago

Option 1

  • With viruses, you can get a long way with a single corrected read. You generated these corrected reads as part of running Canu.
  • You can align these reads with eg minimap2 against a custom FASTA database of the viruses you intended to get. Does it agree with the contig level de novo ? Do you get more information ?

Option 2

Option 3

  • Further improve current metagenomic assembly ( I think option 2 is better).
  • polish with Racon, then Medaka
  • Possibly then scaffold further if possible with Ragtag https://github.com/malonge/RagTag
ADD COMMENT

Login before adding your answer.

Traffic: 3309 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6