question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

READ TAG: n_elements is zero

See original GitHub issue

Hi,

I am trying to run DeepVariant 1.2.0 on a few human samples PacBio HiFi data (about 30x coverage per sample). I first ran my samples through the PEPPER-Margin pipeline r0.4 to get a haplotagged BAM file. Then I ran DeepVariant as follows:

singularity exec -B ${SOME_PATHS} deepvariant_1.2.0.sif bash /opt/deepvariant/bin/run_deepvariant --model_type PACBIO --ref ${PATH_TO_REF} --reads MARGIN_PHASED.PEPPER_SNP_MARGIN.happlotagged.bam  --output_vcf sample.vcf.gz --output_gvcf sample.g.vcf.gz --num_shards 24 --make_examples_extra_args="realign_reads=false,min_mapping_quality=5" --sample_name MYSAMPLE --use-hp-information;

I have two problems:

  1. Right from the beginning (CALL VARIANT MODULE SELECTED), for each interval processed. I get thousands of READ TAG: n_elements is zero messages in the console. What does it mean and is it a problem or just a warning?
  2. I allocate 200GB of RAM for per job and they all seem to systematically fail on memory. I do not recall DeepVariant using that much memory in the past but I might be wrong. Is 200GB too light for a human genome PacBio Hifi 30x coverage dataset?

Thank you for your help, Guillaume

Issue Analytics

  • State:closed
  • Created 2 years ago
  • Comments:9

github_iconTop GitHub Comments

2reactions
GuillaumeHolleycommented, Oct 25, 2021

Just to confirm: removing the tags worked. Thanks!

0reactions
GuillaumeHolleycommented, Oct 22, 2021

Wow, thank you @MariaNattestad @pichuan @williamrowell @kishwarshafin, this all went really quickly. I am currently in the process of generating {fi,fp,ri,rp}-tagless BAM files and will rerun DeepVariant on those. In the meantime, I will close this issue as I am fairly confident you found the solution to the problem. Thank you again!

Read more comments on GitHub >

github_iconTop Results From Across the Web

How to tell XMLReader not to read empty Elements
You can use reader.IsStartElement to check if element is empty or not. while (reader.Read ...
Read more >
XmlTextReader.IsEmptyElement Property (System.Xml)
Gets a value indicating whether the current node is an empty element (for example, <MyElement/>). ... Read(); //Read the start tag. if (reader....
Read more >
XML Elements - W3Schools
Empty XML Elements​​ An element with no content is said to be empty. The two forms produce identical results in XML software (Readers,...
Read more >
The Noscript element - HTML: HyperText Markup Language
The HTML element defines a section of HTML to be inserted if a script ... zero or more <style> elements, and zero or...
Read more >
Reading XML with XmlReader - Yazan Diranieh
public enum XmlNodeType { None = 0, // Read method has not been Element = 1, ... EndElement = 15, // An end...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found