Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
tophat running error figo1019 Bioinformatics 9 05-03-2013 07:18 AM
TopHat error - error running prep_reads ajaffe RNA Sequencing 0 11-29-2012 07:09 AM
velvet running error jjjscuedu Bioinformatics 1 06-05-2012 07:48 AM
tophat Error running running 'prep_reads' victoryhe Bioinformatics 2 10-17-2011 04:53 AM
Error when running TopHat clk RNA Sequencing 3 09-29-2011 12:10 PM

Thread Tools
Old 03-03-2013, 12:57 AM   #1
Location: University of Melbourne

Join Date: Aug 2011
Posts: 10
Default Error running Shrimp aligner

Hi all

I'm trying to map some reads that were found unmapped using bwa with Shrimp agaist chr18.

The command I'm running is as follows:

gmapper-cs -N 8 -E -Q Unmapped-reads.fq chrm-18/chr18-

However I'm getting the following output/error.


SHRiMP 2.1.1b
[ICC Intel(R) C++ g++ 4.1 mode]
Spaced Seeds (weight/span) 111110001111111 (12/15)
111100111001001111 (12/18)
111001001000111001111 (12/21)
1111001000010001001001111 (12/25)
Number of Outputs per Read: 10
Window Generation Mode: 2
Window Length: 140.00%
Window Overlap Length: 20.00%

SW Match Score: 10
SW Mismatch Score: -15
SW Gap Open Score (Ref): -40
SW Gap Open Score (Qry): -40
SW Gap Extend Score (Ref): -7
SW Gap Extend Score (Qry): -7
SW Crossover Score: -14

Window Generation Threshold: 55.00%
SW Vector Hit Threshold: 50.00%
SW Full Hit Threshold: 55.00%

Paired mode: none

Gapless mode: no
Number of threads: 8
Thread chunk size: 1000
Hash Filter Calls: yes
Anchor Width: 8
- Processing genome file [/chrm-18/chr18-homo.fasta]
- Processing contig gi|224589809|ref|NC_000018.9|
Loaded Genome
Automatically trimming genome index lists longer than: 1000
- Processing read file [Unmapped-reads.fq]
gmapper-cs: common/fasta.c:380: bool fasta_get_next_read_with_range(_fasta_t *, read_entry *): Assertion `quality_length <= sequence_length' failed.

Any suggestions is highly appreciated.

Mahtab is offline   Reply With Quote
Old 03-03-2013, 03:11 AM   #2
Senior Member
Location: UK

Join Date: Jan 2010
Posts: 390

It would appear that you have some sequences that have longer qualities than there are bases in the sequence - i.e. the length of the quality is not less than or equal to the sequence length, which is what that error is suggesting.
Bukowski is offline   Reply With Quote
Old 03-03-2013, 04:52 PM   #3
Location: University of Melbourne

Join Date: Aug 2011
Posts: 10

Thanks for you response.

The problem is I'm getting the same error across different sets of input fastq data both synthetic and real data which I had previously mapped using bwa.
Mahtab is offline   Reply With Quote

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:41 AM.

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