Comments (3)
- As far as I can tell, there is only one graph involved here. What do you mean by "merging graphs"?
- If you have a PCR step in your library prep, any standard duplicate removal should work for the BAM output. MarkDuplicates is fine.
from vg.
Thanks for response.
- If I have technical replicates and wish to add more depth before variant calling, I would want to combine the aligned BAM files. Iām asking if there are any accuracy gains by merging graphs from technical replicates versus merging BAMs after conversion of graph to BAM.
- Yes - thanks for confirming.
from vg.
It sounds like you are actually describing merging the reads before mapping. Is that correct? There is no second graph to be merged here. If I'm correct that you're referring to reads, then it shouldn't matter much. There might be a very small improvement in accuracy if you map the replicates separately and then combine the BAMs, because that will allow each replicate to estimate its own fragment length distribution during read mapping. I would expect the difference to be tiny though.
from vg.
Related Issues (20)
- Feature request: Accept jellyfish kmer counts for vg haplotypes HOT 2
- Using vg to analyze transcriptomic data from a mixed background HOT 4
- `vg giraffe` for long reads HOT 1
- Missing '--alt-prefix' option in recent versions of vg deconstruct.
- Problem with vg autoindex with phased VCF HOT 2
- No sample & read group field in BAM output with vg map HOT 5
- I cant get .vg file using vg construct HOT 2
- `vg autoindex` fails on GFA file containing 20k SARS-CoV-2 sequences. HOT 1
- vg giraffe mapping HiFi HOT 1
- error vg autoindex: Input is not sufficient to create indexes HOT 2
- Oddly high mapping quality when aligning reads with Ns HOT 6
- Release vg v1.56.0
- I have a question about gam file. HOT 2
- vg surject: incorrect behavior near polymorphic {repetitive} insertions HOT 4
- Can vg construct a trans-chromosomal graph ? HOT 5
- ERROR: Signal 6 occurred. VG has crashed. HOT 8
- vg construct will let you create a graph without reference paths, with different path names than Cactus
- VG convert to gfa and back will drop path data. HOT 7
- ERROR: Signal 6 occurred. VG has crashed. Terminate called after throwing an instance of 'std::bad_alloc' what(): std::bad_alloc HOT 1
- VG has crashed
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 vg.