jts / nanopolish

Signal-level algorithms for MinION data
MIT License
566 stars 159 forks source link

Empty vcf #885

Closed madelaam closed 3 years ago

madelaam commented 3 years ago

Hello,

I am trying to get a polished sequence after canu: python3 '/home/grid/nanopolish/scripts/nanopolish_makerange.py' '/data/Varicovorax/Varicovorax/20210204_1720_X2_FAP14911_6902f034/Analysis/Canu/Varicovorax.canu.unitigs.fasta' | parallel --results '/data/Varicovorax/Varicovorax/20210204_1720_X2_FAP14911_6902f034/Analysis/Nanopolish/voricovorazx.results' -P 2 \ '/home/grid/nanopolish/nanopolish' variants --consensus -o '/data/Varicovorax/Varicovorax/20210204_1720_X2_FAP14911_6902f034/Analysis/Nanopolish/voricovorax.polished.{1}.vcf' -w {1} -r '/data/Varicovorax/Varicovorax/20210204_1720_X2_FAP14911_6902f034/Analysis/varicovorax.fasta' -b '/data/Varicovorax/Varicovorax/20210204_1720_X2_FAP14911_6902f034/Analysis/Minimap/voricovorax.sorted.bam' -g '/data/Varicovorax/Varicovorax/20210204_1720_X2_FAP14911_6902f034/Analysis/Canu/Varicovorax.canu.unitigs.fasta' -t 2 --min-candidate-frequency 0.1

And I getting empty vcf files:

fileformat=VCFv4.2

nanopolish_window=tig00000001:0-43935

INFO=

INFO=

INFO=

INFO=

INFO=

INFO=

INFO=

INFO=

INFO=

INFO=

FORMAT=

CHROM POS ID REF ALT QUAL FILTER INFO FORMAT sample

Thank you.

jts commented 3 years ago

In the error logs you should see messages printed that say "post-run summary" or similar. Could you paste some examples of those?

madelaam commented 3 years ago

Hello,

I cannot find the posr-run summary, I have just attached the folder where all the output is.

Thank you.


From: Jared Simpson notifications@github.com Sent: 22 February 2021 10:24 PM To: jts/nanopolish nanopolish@noreply.github.com Cc: MEDINA, Adela (VIAPATH LLP) adela.medina@nhs.net; Author author@noreply.github.com Subject: Re: [jts/nanopolish] Empty vcf (#885)

In the error logs you should see messages printed that say "post-run summary" or similar. Could you paste some examples of those?

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://github.com/jts/nanopolish/issues/885#issuecomment-783719131, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AS6F5XIMSRJZ5RL5YRZKPPTTALKRDANCNFSM4YAL3B6A.

This message originated from outside of NHSmail. Please do not click links or open attachments unless you recognise the sender and know the content is safe.


This message may contain confidential information. If you are not the intended recipient please inform the sender that you have received the message in error before deleting it. Please do not disclose, copy or distribute information in this e-mail or take any action in relation to its contents. To do so is strictly prohibited and may be unlawful. Thank you for your co-operation.

NHSmail is the secure email and directory service available for all NHS staff in England and Scotland. NHSmail is approved for exchanging patient data and other sensitive information with NHSmail and other accredited email services.

For more information and to find out how you can switch, https://portal.nhs.net/help/joiningnhsmail

madelaam commented 3 years ago

Sorry, I though I attached the folder but I did not. Here it is.

Nanopolish.zip

Thank you!

madelaam commented 3 years ago

I done it, Thank you!