Forum:WGBS alignment with BitMapperBS
0
2
Entering edit mode
2.8 years ago
Shicheng Guo ★ 9.0k

Now, it is time to try a faster WGBS alignment tool: BitMapperBS

./bitmapperBS --index human_genome.hg38.fa
trim_galore --paired read_1.fq read_2.fq
./bitmapperBS --search human_genome.hg38.fa --seq1 read1.fq --seq2 read2.fq  -t 6 --bam -o output.bam
MethylDackel extract reference_genome.fa alignments.bam
bismark BitMapperBS Forum • 1.2k views
ADD COMMENT
2
Entering edit mode

I have to admit being surprised by how fast bitmapperBS is. For those wondering you can install it via conda.

ADD REPLY
2
Entering edit mode

I love all these chinese comments spread all over the code ;-)

ADD REPLY
1
Entering edit mode

I'm still adding some new features to bitmapperbs, so these chinese comments have not been removed ...

ADD REPLY
1
Entering edit mode

Thank you a lot : ). Following your suggestions, I'm trying to make BitMapperBS calculate MAPQ like BWA, and then submit it back to PCB ..

ADD REPLY
1
Entering edit mode

Awesome!

ADD REPLY
1
Entering edit mode

Thank you for your poster. I just worry about how to publish BitMapperBS ...

ADD REPLY
1
Entering edit mode

Nature Protocol? Nature Method? Nucleic Acid Research? You can suggest Devon Ryan as one of the preferred reviewers.

ADD REPLY
1
Entering edit mode

It had been rejected by several reviewers. Some of them said they didn't think BitMapperBS is better than existing methods, some other people said it is lack of novelty, and some others said bs-seq aligner is not important ... Now I just want to publish it.

ADD REPLY
1
Entering edit mode

Maybe you didn't suggest the best reviewers. Sure. different people have different comments to same tools. You need to find some experts who appreciate your tools.

ADD REPLY
1
Entering edit mode

Thank you so much! You know in most cases editors do not use the reviewers recommended by the authors, and some editors may directly reject submission without big name. Good luck to me ...

ADD REPLY

Login before adding your answer.

Traffic: 1292 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6