Question: Nanopolish makerange for smaller contigs
0
gravatar for SUDOsundu
12 months ago by
SUDOsundu0
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 
    ATGC
    > contig2 
    ATGCA
    > contig3 
    ATGCAT

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 nanopolish_makerange.py 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.

Help
Access

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
_