Comments (8)
-F 0xB00
rather than multiple -F
values, but otherwise yes. I'll make that the default in the next bug fix version.
from methyldackel.
Thanks for reporting the bug!
from methyldackel.
-F 0x400
is exactly the opposite of --keepDupes
. -F
in MethylDackel is the same as -F
in samtools view
, so you're then instructing it to ignore duplicates (but then include secondary and supplemental alignments, since you've overridden the default value of 0xF00
or 3840).
from methyldackel.
Thanks. so, If I need to include the duplicates, I can do:
MethylDackel --keepDupes
# or
MethylDackel -F 0x100 -F 0x200 -F 0x800
what if one uses both --keepDupes
and -F
together? and the order of the flag also matters seems.
e.g. This seems to give different results.
MethylDackel -F 0x400 --keepDupes
MethylDackel --keepDupes -F 0x400
Thanks,
Tommy
from methyldackel.
The -F
option is processed first when reads are filtered, so if you filter it then --keepDupes
is ignored. This looks like a bug, actually, since I'm now not sure that --keepDupes
is properly doing anything since it isn't changing the -F
default. I think I should change the -F
default to 0xB00
, I'll have to test that.
from methyldackel.
Thanks for the detail. I just tested
MethylDackel -F 0x100 -F 0x200 -F0x800
MethylDackel --keepDupes -F 0x100 -F 0x200 -F0x800
I thought these two are the same. -F 0x100 -F 0x200 -F0x800
will keep the duplicates, but one needs to specify --keepDupes
as well for the counting.
For my purpose, I will need to do to keep the duplicates but discarding failed qc, secondary and supplementary reads:
MethylDackel --keepDupes -F 0x100 -F 0x200 -F0x800
right?
Thanks,
Tommy
from methyldackel.
great! thanks for the explaining and user support!
from methyldackel.
Sorry for the delay on this, it's now fixed in the 0.4.0 release.
from methyldackel.
Related Issues (20)
- Genome browsing from MethylDackel bedGraphCpG file
- per-fragment methylation HOT 2
- mbias result is different between bismark and bwameth output HOT 1
- CURL_OPENSSL conflict with samtools HOT 3
- Installation failure: "bigWig.h: No such file or directory" HOT 1
- Clarification on definition of "unmethylated C" HOT 1
- Coverage of C sites HOT 1
- mbias HOT 1
- Does indel effect the methylation calling or C context determination HOT 1
- Positions in cytosine_report did not match the regions in providing bed file
- about M-bias HOT 1
- Could not repeat a CpG extraction with the same reference file and its index HOT 1
- Mixed up reads within bam file
- How to index genome file for MethylDackel? HOT 1
- Confused regarding CTOT, CTOB. Are there suggested values? HOT 3
- Issue running MethylDackel extract in parallel mode using minConversionEfficiency
- Question about minimum coverage
- Different CpG calls when using different regions of inclusion for Methyldackel extract
- Alignment trimming for Soft Clip reads?
- Can MethylDackel extract per-base methylation metrics for TAPS?
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 methyldackel.