Question: Nanopolish makerange for smaller contigs
gravatar for SUDOsundu
12 months ago by
SUDOsundu0 wrote:

My aim is to characterize plant viral genome from nanopore reads. I am doing reference based assembly. First mapping the reads to the reference genome with minimap2 and polishing with nanopolish. In nanopolish variants, the reference file which I have to use contains 3 contigs

For example

    > contig1 
    > contig2 
    > contig3 

From what I understood, I can do

nanopolish variants -r reads.fasta  -b reads.sorted.bam -g draft.fa -w "contig1:0-4" -o polished.vcf

If it is larger genome can be used to split. By default it is splitting a single contig into 50kb of segments. But my reference file is a viral genome which has approx 1kb of 3 contigs

If I want to to call variants for all three contig in parallel and to result in polished.vcf for all three contigs how to change the parameters?

Help please. Thank you.

assembly genome • 312 views
ADD COMMENTlink written 12 months ago by SUDOsundu0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2262 users visited in the last hour