SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools: bam_pileup.c:112: resolve_cigar2: Assertion `s->k < c->n_cigar' failed SEQquestions Bioinformatics 27 02-15-2013 09:57 AM
Bfast localalign error dg.pooja General 1 01-31-2011 05:35 PM
bfast localalign error nimmi Bioinformatics 6 01-05-2011 08:40 AM
BFAST match/localalign help Esther Bioinformatics 3 08-04-2010 05:19 AM
maq map - Assertion 'matches' failed AnamikaDarwin Bioinformatics 1 12-10-2008 07:56 AM

Reply
 
Thread Tools
Old 09-02-2011, 05:31 AM   #1
seeker
Member
 
Location: Zurich

Join Date: Jan 2011
Posts: 26
Default Assertion failed error in BFAST localalign

Hi all, I'm getting the following error when trying to combine bfast and bwa mapped reads. Any suggestions would be appreciated.

command:
bfast localalign -f rna.fasta -1 p1.f3.reads.bmf -2 p1.f5.reads.bfm -A 1 -t -U -n 4 > p1.reads.baf

error:
Assertion failed: (0 == strcmp(m->readName, r2_name)), function RGMatchesRead_2bmf, file RGMatches.c, line 51.
Abort trap
Assertion failed: (0 == strcmp(m->readName, r2_name)), function RGMatchesRead_2bmf, file RGMatches
seeker is offline   Reply With Quote
Old 09-02-2011, 07:02 AM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Looks like your reads are not in the same order within the two files (also p1.f5.reads.bfm should be p1.f5.read.bmf). Take a look at your input FASTQ files to make sure the read names are in the same order.
nilshomer is offline   Reply With Quote
Old 09-02-2011, 10:02 AM   #3
seeker
Member
 
Location: Zurich

Join Date: Jan 2011
Posts: 26
Default

The first and last 100 lines of reads are in the same order. I've renamed the f5 files to bmf (which was a typo when making the files). Could that be responsible?
seeker is offline   Reply With Quote
Old 09-02-2011, 04:17 PM   #4
seeker
Member
 
Location: Zurich

Join Date: Jan 2011
Posts: 26
Default

Changed the name and still the same error. How could I ensure that all the names are in the same order?
seeker is offline   Reply With Quote
Old 09-02-2011, 06:25 PM   #5
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Where are you getting your FASTQ from?
nilshomer is offline   Reply With Quote
Old 09-02-2011, 07:05 PM   #6
seeker
Member
 
Location: Zurich

Join Date: Jan 2011
Posts: 26
Default

From the solid2fastq script that came with bfast+bwa. I made a seperate fastq for each of the F3 and F5 reads rather than paired.
seeker is offline   Reply With Quote
Old 09-02-2011, 07:07 PM   #7
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Ah, you have to make them together. You cannot convert them unpaired then input them as paired.
nilshomer is offline   Reply With Quote
Old 09-02-2011, 09:33 PM   #8
seeker
Member
 
Location: Zurich

Join Date: Jan 2011
Posts: 26
Default

Oh, I was following the flow from here: http://seqanswers.com/forums/showthr...ight=bfast+bwa

So, you make the fastq paired, then map using the same fastq file for the f3 and f5 seperately?
seeker is offline   Reply With Quote
Reply

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 05:58 AM.


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