Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Bismark - A New Tool for Mapping and Analysis of Bisulfite-Seq Data fkrueger Bioinformatics 649 10-05-2018 02:43 AM
BSMAP: whole genome Bisulfite Sequence MAPping program wei Epigenetics 4 03-20-2014 02:13 PM
which mapping program is good for RainDance Data JCYAO Bioinformatics 5 01-26-2010 12:29 PM
mapping program with Hadoop wykim SOLiD 2 12-20-2009 12:49 AM
whole genome Bisulfite Sequence MAPping program wei Bioinformatics 0 08-07-2009 03:46 PM

Thread Tools
Old 01-12-2012, 02:37 PM   #1
Junior Member
Location: houston

Join Date: Jun 2010
Posts: 6
Default update on bisulfite mapping program: BSMAP-2.4

This is an announcement of the release of BSMAP v2.4, a powerful bisulfite mapping program. This version improved a lot on the run time performance, while maintaining high accuracy and flexibility of previous versions.

Using 8 threads, BSMAP-2.4 can map 28M 76nt pair-end WGBS reads to the human genome in about 2 hours (allow up to 6 mismatches), including the genome indexing time. The memory usage is ~9GB. (CPU: Intel Xeon X5690)

We tested the mapping accuracy using simulated bisulfite reads. BSMAP have significantly higher mapping accuracy than most other bisulfite mapping programs, especially for reads with more than 3 mismatches.

BSMAP is freely available at

Main features of BSMAP-2.4:

1. Reads are directly mapped to original reference genome sequence, no need to preprocess the reads and reference genome to convert C to T.

2. Support both whole genome bisulfite sequencing (WGBS) mode and reduced representation bisulfite sequencing (RRBS) mode. In RRBS mode, reads are guaranteed to be mapped to digestion sites to increase accuracy. The digestion site information can be also changed to support different digestion enzymes.

3. Support both "Lister protocol" (sequence 2 forward strands only) and "Cokus protocol" (sequence all 4 bisulfite converted strands)

4. Support trimming adapter sequences and low quality nucleotides from 3'end of reads

5. Allow trade off between speed/memory usage/mapping sensitivity. For human genome, the RRBS mode uses ~3GB. In WGBS mode, the typical memory usage is ~9GB, but can be as low as 5GB.

6. Allow alignment for other nucleotide transitions, for example, can be set to detect the A=>I(G) transition in RNA editing.

7. Include down stream script to extract methylation ratios from mapping results.

8. Fasta/Fastq/SAM format input, text/SAM output. single/pair-end mapping. read length up to 144nt, max 15 mismatches allowed.

We encourage you try this new version for short bisulfite reads mapping.
Any comments/suggestions/bug reports will be appreciated.

Thank you,

Yuanxin Xi
yxibcm is offline   Reply With Quote
Old 04-13-2013, 12:03 AM   #2
Location: China

Join Date: Jan 2011
Posts: 26

Dear Xi,
here I have a problem in installation of BSMAP, my system is opensuse 12.3 x86_64 GNU/Linux. and when I make the software, the errors:
xiao@localhost:~/Software/bsmap-2.74> make
g++ -DMAXHITS=1000 -DTHREAD -funroll-loops -Lsamtools -Isamtools -Lgzstream -Igzstream -O3 -m64 -c param.cpp -o param.o
param.cpp: In constructor ‘Param::Param()’:
param.cpp:8:20: error: ‘_SC_NPROCESSORS_ONLN’ was not declared in this scope
param.cpp:8:40: error: ‘sysconf’ was not declared in this scope
make: *** [param.o] Error 1

could you tell me what I should do to resolve such a error.

Thank you.

xfh is offline   Reply With Quote
Old 04-29-2013, 07:17 AM   #3
Junior Member
Location: heidelberg

Join Date: Feb 2013
Posts: 1

adding "#include<unistd.h>" in the beginning of param.cpp would solve this problem.
guang920 is offline   Reply With Quote
Old 05-02-2013, 06:38 PM   #4
Location: China

Join Date: Jan 2011
Posts: 26

OKay, thank you!
xfh is offline   Reply With Quote
Old 05-03-2016, 06:57 AM   #5
Junior Member
Location: delhi

Join Date: Apr 2016
Posts: 8

Dear Xi

i used BSmap for mapping and it gave better mapping efficiency as compared to other tools. But there is problem to extract methylation information from .txt file.
plz suggest to solve it
alim123 is offline   Reply With Quote

bisulfite sequencing, bsmap, dna methylation, rna editing, sequence alignment

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 04:11 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO