Comments (10)
Going to make idxdepth ignore ALT contigs.
from paragraph.
Binned coverage for bam files is disabled in idxdepth. That's why you cannot get binnedcov.txt for your bam file. I have no idea why @pkrusche applies this. While waiting for his reply, I'll try to figure it out.
from paragraph.
I am getting the same error too
[2019-07-15 16:42:23.776] [idxdepth] [12150] [info] Thread 46971021932288 done estimating depth for chrUn_JTFH01001987v1_decoy ; DP = 579.61 after 5230 us
[2019-07-15 16:42:23.776] [idxdepth] [12150] [info] Thread 46971021932288 estimating depth for chrUn_JTFH01001996v1_decoy
[2019-07-15 16:42:25.712] [idxdepth] [12152] [info] Thread 46971442169600 done estimating depth for chrUn_JTFH01001989v1_decoy ; DP = 0 after 4112.74 us
[2019-07-15 16:42:25.712] [idxdepth] [12152] [info] Thread 46971442169600 estimating depth for chrUn_JTFH01001997v1_decoy
[2019-07-15 16:42:25.748] [idxdepth] [12147] [info] Thread 46971015628544 done estimating depth for chrUn_JTFH01001990v1_decoy ; DP = 0 after 4109.9 us
[2019-07-15 16:42:25.748] [idxdepth] [12147] [info] Thread 46971015628544 estimating depth for chrUn_JTFH01001998v1_decoy
[2019-07-15 16:42:25.948] [idxdepth] [12153] [info] Thread 46971444270848 done estimating depth for chrUn_JTFH01001992v1_decoy ; DP = 0.44 after 3825.24 us
[2019-07-15 16:42:25.948] [idxdepth] [12153] [info] Thread 46971444270848 estimating depth for HLA-A*01:01:01:01
[2019-07-15 16:42:25.968] [idxdepth] [12149] [info] Thread 46971019831040 done estimating depth for chrUn_JTFH01001991v1_decoy ; DP = 1.87 after 3982.76 us
[2019-07-15 16:42:25.968] [idxdepth] [12149] [info] Thread 46971019831040 estimating depth for HLA-A*01:01:01:02N
[2019-07-15 16:42:26.218] [idxdepth] [12154] [info] Thread 46971643496192 done estimating depth for chrUn_JTFH01001993v1_decoy ; DP = 0.03 after 2784.95 us
[2019-07-15 16:42:26.218] [idxdepth] [12154] [info] Thread 46971643496192 estimating depth for HLA-A*01:01:38L
[2019-07-15 16:42:26.244] [idxdepth] [12148] [info] Thread 46971017729792 done estimating depth for chrUn_JTFH01001994v1_decoy ; DP = 2.34 after 2709.57 us
[2019-07-15 16:42:26.244] [idxdepth] [12148] [info] Thread 46971017729792 estimating depth for HLA-A*01:02
[2019-07-15 16:42:26.482] [idxdepth] [12151] [info] Thread 46971240843008 done estimating depth for chrUn_JTFH01001995v1_decoy ; DP = 2.66 after 2946.74 us
[2019-07-15 16:42:26.482] [idxdepth] [12151] [info] Thread 46971240843008 estimating depth for HLA-A*01:03
[2019-07-15 16:42:26.488] [idxdepth] [12150] [info] Thread 46971021932288 done estimating depth for chrUn_JTFH01001996v1_decoy ; DP = 1.13 after 2711.84 us
[2019-07-15 16:42:26.488] [idxdepth] [12150] [info] Thread 46971021932288 estimating depth for HLA-A*01:04N
terminate called after throwing an instance of 'std::invalid_argument'
what(): stoll
Here's my argument
idxdepth -b NA19240.alt_bwamem_GRCh38DH.20150715.YRI.high_coverage_markedDupsRemoved.bam \
-o NA19240.paragraph_cov.txt \
-O NA19240.paragraph_binnedcov.txt \
-r GRCh38_full_analysis_set_plus_decoy_hla.fa \
--threads 8
Is there a way to not process the decoy contigs?
Thank you.
from paragraph.
I tried to use these options
--autosome-regex '(chr)?[1-9][0-9]?' --sex-chromosome-reg '(chr)?[XY]?'
But it started to process the ALT contigs in hg38.
[2019-07-16 07:00:39.498] [idxdepth] [34604] [info] Thread 47330071938816 done estimating depth for chr22 ; DP = 52.7 after 52444 us
[2019-07-16 07:00:39.498] [idxdepth] [34604] [info] Thread 47330071938816 estimating depth for chr1_KI270706v1_random
[2019-07-16 07:00:44.000] [idxdepth] [34604] [info] Thread 47330071938816 done estimating depth for chr1_KI270706v1_random ; DP = 86.97 after 4501.77 us
[2019-07-16 07:00:44.000] [idxdepth] [34604] [info] Thread 47330071938816 estimating depth for chr1_KI270707v1_random
[2019-07-16 07:00:45.477] [idxdepth] [34600] [info] Thread 47329517483776 done estimating depth for chrM ; DP = 15677.6 after 17895.2 us
[2019-07-16 07:00:45.477] [idxdepth] [34600] [info] Thread 47329517483776 estimating depth for chr1_KI270708v1_random
[2019-07-16 07:00:46.694] [idxdepth] [34599] [info] Thread 47329515382528 done estimating depth for chr18 ; DP = 80.98 after 70156.3 us
[2019-07-16 07:00:46.694] [idxdepth] [34599] [info] Thread 47329515382528 estimating depth for chr1_KI270709v1_random
[2019-07-16 07:00:48.232] [idxdepth] [34601] [info] Thread 47329602176768 done estimating depth for chr17 ; DP = 65.13 after 74335.1 us
[2019-07-16 07:00:48.233] [idxdepth] [34601] [info] Thread 47329602176768 estimating depth for chr1_KI270710v1_random
[2019-07-16 07:00:51.749] [idxdepth] [34604] [info] Thread 47330071938816 done estimating depth for chr1_KI270707v1_random ; DP = 89.22 after 7748.92 us
[2019-07-16 07:00:51.749] [idxdepth] [34604] [info] Thread 47330071938816 estimating depth for chr1_KI270711v1_random
from paragraph.
Thanks please update this thread. For now I am using bedtools genomecov
from paragraph.
v2.3 released with idxdepth excluding ALT contig as default. Let me know if this problem still exists with v2.3.
from paragraph.
Seems to work now, except I do not see the -O NA19240.paragraph_binnedcov.txt
output.
It skips the ALT contigs but I get this error
[2019-07-29 09:54:42.694] [idxdepth] [27225] [critical] BAM index reading not available yet.
The NA19240.paragraph_cov.txt
is present and looks good.
Thanks for quickly addressing this issue!
from paragraph.
Estimating coverage from a CRAM index is done by proxy (number of bytes in chunk), which is easy to implement (the method was shown here: https://github.com/brentp/goleft/tree/master/indexcov). On a BAM index this should be possible also, but will have to be done differently -- it wasn't implemented because there was no need when I wrote this (all relevant files were CRAMs).
from paragraph.
Ok I thought the "binned" coverage was something similar to CNVnator, which calculates read depth in bins of a specific (user-defined) length.
So it seems to not be relevant to this concept, the output for idxdepth is fine now. Thanks for all your help!
from paragraph.
Thank you all for your discussion. v2.4 dded another warning message in idxdepth for BAM file for this issue.
from paragraph.
Related Issues (20)
- Exception: Different padding base for REF and ALT HOT 6
- Can paragraph be used for indel from 2 bp to 30 bp? HOT 1
- ValueError: Invalid VariantRecord. Number of samples does not match header HOT 2
- Error with idxdepth: "Assertion failed: _impl->header_contig_map.count(chr) != 0" HOT 1
- Missing key SEQ for <INS> HOT 6
- --vcf-split option with no description
- subprocess.CalledProcessError HOT 2
- grmpy error: [E::cram_itr_query]
- index file
- How to merge multi-samples SVs and obtain breakpoints for genotyping a population
- no BGZF EOF marker
- Install paragraph
- Stop using Werror and Wall
- Genotyping for SNP
- idxdepth regex option not working
- Problem starting the script multigrmpy.py HOT 2
- Error when working with `--ins-info-key` HOT 1
- Add support for VCFv4.2 breakend notation
- Issue's with Native Build and Boost
- Format error in vcf line: HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from paragraph.