nanopolish: error: could not parse model string:

Running nanopolish version 0.7.1

After assembly with canu, I’m trying to run nanopolish on a dataset, but I keep getting the error: could not parse model string:

I’m working with Albacore basecalled 2D data I ran on a r9.4 flow cell. I used nanopolish to extract the fastq/fasta reads. I looked at the fast5 files and I noticed that the the models are:

UniqueGlobalKey:context_tags: local_bc_comp_model = complement_r9.4_250bps_5mer local_bc_temp_model = template_r9.4_250bps_5mer

Do I need to import new models for nanopolish to use? I noticed there are r9 250bps 5mer models, but not for r9.4. Just a guess on my part.

INPUT command:

 python ~/nanopolish/scripts/nanopolish_makerange.py draft.fa | parallel --results nanopolish.results -P 8 \
> nanopolish variants --consensus polished.{1}.fa -w {1} -r reads.fa -b reads.sorted.bam -g draft.fa -t 4 --min-candidate-frequency 0.1

Full output:

##fileformat=VCFv4.2
##INFO=<ID=TotalReads,Number=1,Type=Integer,Description="The number of event-space reads used to call the variant">
##INFO=<ID=SupportFraction,Number=1,Type=Float,Description="The fraction of event-space reads that support the variant">
##INFO=<ID=BaseCalledReadsWithVariant,Number=1,Type=Integer,Description="The number of base-space reads that support the variant">
##INFO=<ID=BaseCalledFraction,Number=1,Type=Float,Description="The fraction of base-space reads that support the variant">
##INFO=<ID=AlleleCount,Number=1,Type=Integer,Description="The inferred number of copies of the allele">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  sample
[fai_load] build FASTA index.
error: could not parse model string:

About this issue

  • Original URL
  • State: closed
  • Created 7 years ago
  • Comments: 15 (14 by maintainers)

Most upvoted comments

Thanks, that is helpful. I think I know what the problem is now and I can provide a fix later.