Comments (4)
Hello,
I think you are right, default alignment parameters for hicstuff are defined here: https://github.com/koszullab/hicstuff/blob/830d04fffbdb14957ba95f196d7bcdd690d84da9/hicstuff/pipeline.py#L89-L101
Generally speaking, we use local
alignment because a global alignment would not allow to partially map chimeric Hi-C reads. However if you want to use different parameters, you can map the reads with whatever software / parameters you prefer and provide the output (name-sorted) BAM files to hicstuff instead of starting from fastq files.
There is a small example starting from an intermediate BAM file in the hicstuff docs: https://hicstuff.readthedocs.io/en/latest/notebooks/demo_cli.html#Starting-from-intermediate-files
I think this is what they did in the instaGRAAL paper, perhaps @baudrly or @nadegeguiglielmoni can confirm ?
Instagraal filters fragments from the Hi-C matrix with low coverage to improve results, this is the --coverage-std
option. One standard dev is indeed the default value.
Hope that helps
from instagraal.
Hello ,
Thanks for your reply, it really helps . In other words , you just run hicstuff pipeline
with the default parameters when you perform your experiment in the paper ? And not setting special parameter ? Am I right ?
Anyway , thanks for you help ,again.
from instagraal.
Hello,
You can indeed run the mapping step directly with hicstuff pipeline
. For 150 bp reads, you should definitely use the parameter --iterative
. You should also put the enzyme used for Hi-C library preparation as the parameter --enzyme
. By default, the mapping tool used is bowtie2, but you can also use bwa instead if you prefer with the parameter --aligner bwa
.
So if you run hicstuff on a library prepared with DpnII and using bowtie2 for mapping, the command would be:
hicstuff pipeline --enzyme DpnII --iterative --outdir hic_folder --genome assembly.fasta R1.fastq R2.fastq
You can also look at the repository with command lines used in the paper https://github.com/koszullab/ectocarpus_scripts
from instagraal.
Hi @nadegeguiglielmoni ,
Thanks for you links, I viewed the scripts and got to know the next step .
Have a nice day and I will close the issue.
from instagraal.
Related Issues (20)
- Should the initial reference assembly and final polished genome have equal sequence lengths? HOT 1
- How to use genetic linkage data?
- What does "--bomb" do? HOT 1
- Description of new_info_frags.txt file HOT 1
- instagraal-polish: problem with -j option HOT 1
- OSError: [Errno 28] No space left on device HOT 2
- module.cpp:1:10: fatal error: boost/python.hpp: No such file or directory HOT 3
- Instagraal docker run gives error about pyramid directory HOT 2
- KeyError on import reference genome HOT 3
- tools option not showing HOT 41
- Instagraal polish - KeyError HOT 2
- instagraal-polish HOT 8
- How to separate "merged" chromosomes? HOT 2
- Module 'numpy' has no attribute 'asscalar'
- pycuda issue with CUDA toolkit 12.0 HOT 8
- Error with --simple HOT 2
- pycuda installation HOT 2
- pycuda version 2022 not supported HOT 1
- usage of __shfl_down is deprecated HOT 2
- FileNotFoundError: [Errno 2] No such file or directory: '/xxx/pyramids' HOT 4
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 instagraal.