SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Building bowtie index with mirBase hairpin.fa file Gators RNA Sequencing 6 05-07-2015 11:43 AM
bowtie index problem (bowtie-build and then bowtie-inspect) tgenahmet Bioinformatics 4 09-10-2013 11:51 AM
Building bfast index with btestindexes bre Bioinformatics 35 09-29-2011 04:34 AM
Bowtie mapping problem harrike Bioinformatics 1 05-09-2011 06:29 AM
tophat-bowtie building index repinementer Bioinformatics 1 07-17-2010 10:53 PM

Reply
 
Thread Tools
Old 08-04-2010, 05:02 PM   #1
Gangcai
Member
 
Location: Shanghai, China

Join Date: Nov 2009
Posts: 30
Question strange bowtie index building and mapping problem

Hello everyone,
I have met a quite strange problem when running bowtie for my SOLID sequencing result.
I have succeeded to build bowtie colorspace index for chimpanzee whole genome, but when I run bowtie , it showed the error like:
"run_bowtie.sh: line 32: 1898 Segmentation fault $bowtie -C --phred33-quals -n 3 -a -m 3 --best --strata -p 4 -S $index -q $infile $out_dir/$out_name
"
Then I checked the output, the output just contain sam header and a few line of alignment result. The strange thing is that the sam header is like this(chromosome name is incorrect):
@HD VN:1.0 SO:unsorted
@SQ SN:IB<8D><EA>^RP<F3>{<<FD><FA><93>ދ<FC>^X%^M<C3>^R^V<DA>A LN:135001995
@SQ SN:1 LN:8402541
@SQ SN:2 LN:134204764
@SQ SN:3 LN:8412303
@SQ SN:4 LN:135371336
@SQ SN:5 LN:2259969
@SQ SN:6 LN:115868456
@SQ SN:7 LN:9231045
@SQ SN:8 LN:107349158
@SQ SN:9 LN:2108736
@SQ SN:10 LN:100063422
@SQ SN:11 LN:3087076
@SQ SN:12 LN:90682376
@SQ SN:13 LN:6370548
@SQ SN:14 LN:83384210
@SQ SN:15 LN:5078517
@SQ SN:16 LN:77261746
@SQ SN:17 LN:1841920
@SQ SN:18 LN:64473437
@SQ SN:19 LN:2407237
......

I have done same procedure for human, the sam header is :
@HD VN:1.0 SO:unsorted
@SQ SN:chr10 LN:135534747
@SQ SN:chr11 LN:135006516
@SQ SN:chr11_gl000202_random LN:40103
@SQ SN:chr12 LN:133851895
@SQ SN:chr13 LN:115169878
@SQ SN:chr14 LN:107349540
@SQ SN:chr15 LN:102531392
@SQ SN:chr16 LN:90354753
@SQ SN:chr17_ctg5_hap1 LN:1680828
.....

So it seems bowtie either didn't build a correct colorspace index or read the index incorrectly. By the way, I have checked the chimpanzee genome I used, it seems OK(the chromosome name is chr1,chr2a, .....chrX, ...).

Have anybody met the same problem before? I am quite precious to hear your suggestion.
Gangcai is offline   Reply With Quote
Reply

Tags
bowtie, software, solid

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 02:19 AM.


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