Nanopolish variant calling. Need to specify input region of the reference genome
1
0
Entering edit mode
5.4 years ago
SUDOsundu ▴ 80

I mapped the reads with the reference fasta with minimap2. Out put: BAM file

Then I used the BAM, reference fasta, reads to nanopolish variant calling. And I got this error.

Error: genome has multiple contigs, please use -w to specify input region

My reference fasta contains 6 fasta sequences. It has no chromosome locations in it. Without knowing the chromosome locations (format: <chromsome_name>:<start>-<end>)) how do I specify the input region?

genome assembly • 1.4k views
ADD COMMENT
2
Entering edit mode
5.4 years ago

Chromosomes are irrelevant. The tool looks at contigs regardless of how well assembled these are, and considers every sequence in your reference fasta a contig.

Say that you have the following fasta file:

>my_fabulous_contig  
ATGCGCGCGC    
>and_another_contig  
GCTTACTTGTGTC

Then -w my_fabulous_contig:2-5 would be perfectly valid.

ADD COMMENT

Login before adding your answer.

Traffic: 2309 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