SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools "is recognized as '*'" "truncated file" error axiom7 Bioinformatics 3 11-26-2014 03:53 AM
Relatively large proportion of "LOWDATA", "FAIL" of FPKM_status running cufflink ruben6um Bioinformatics 3 10-12-2011 01:39 AM
The position file formats ".clocs" and "_pos.txt"? Ist there any difference? elgor Illumina/Solexa 0 06-27-2011 08:55 AM
"Systems biology and administration" & "Genome generation: no engineering allowed" seb567 Bioinformatics 0 05-25-2010 01:19 PM
SEQanswers second "publication": "How to map billions of short reads onto genomes" ECO Literature Watch 0 06-30-2009 12:49 AM

Reply
 
Thread Tools
Old 09-12-2011, 08:30 AM   #1
sdvie
Member
 
Location: Spain

Join Date: Jul 2010
Posts: 68
Default bfast postprocess "AlignedEntryGetAlignment" error

Dear all,

unexpectedly, I am getting the following error when trying to run bfast postprocess (alignment of 95 bp Illumina paired-end reads):

Code:
bfast+bwa-0.6.4e/bin/bfast postprocess -i t_8.baf -f hg19_indices/hg19.fa -a 3 -z -A 0 -R -O 1 -n 4 -U > t_8_alternative.sam
************************************************************
Checking input parameters supplied by the user ...
Validating fastaFileName /projects/bg/genomed/data/hg19_indices_ucsc/hg19.fa. 
Validating alignFileName t_8.baf. 
Input arguments look good!
************************************************************
************************************************************
Printing Program Parameters:
programMode:			[ExecuteProgram]
fastaFileName:			/projects/bg/genomed/data/hg19_indices_ucsc/hg19.fa
alignFileName:			t_8.baf
algorithm:			[Best Score]
space:				[NT Space]
unpaired:			[Using]
reversePaired:			[Not Using]
avgMismatchQuality:		10
scoringMatrixFileName:		[Not Using]
randomBest:			[Not Using]
minMappingQuality:		-2147483648
minNormalizedScore:		-2147483648
pairingStandardDeviation:	2.000000
gappedPairingRescue		[Not Using]
numThreads:			4
queueLength:			50000
outputFormat:			[SAM]
outputID:			[Not Using]
RGFileName:			[Not Using]
timing:				[Not Using]
************************************************************
************************************************************
Reading in reference genome from /projects/bg/genomed/data/hg19_indices_ucsc/hg19.fa.nt.brg.
In total read 93 contigs for a total of 3137161264 bases
************************************************************
Postprocessing...
************************************************************
In function "AlignedEntryGetAlignment": Fatal Error[OutOfRange]. Message: Could not get reference sequence.
 ***** Exiting due to errors *****
************************************************************
I do not know what to make of this error - I ran a lot of analysis with the same reference sequence (in the same folder).

Suggestions are very much appreciated.

thanks and cheers,
Sophia
sdvie is offline   Reply With Quote
Old 09-12-2011, 09:23 AM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by sdvie View Post
Dear all,

unexpectedly, I am getting the following error when trying to run bfast postprocess (alignment of 95 bp Illumina paired-end reads):

Code:
bfast+bwa-0.6.4e/bin/bfast postprocess -i t_8.baf -f hg19_indices/hg19.fa -a 3 -z -A 0 -R -O 1 -n 4 -U > t_8_alternative.sam
************************************************************
Checking input parameters supplied by the user ...
Validating fastaFileName /projects/bg/genomed/data/hg19_indices_ucsc/hg19.fa. 
Validating alignFileName t_8.baf. 
Input arguments look good!
************************************************************
************************************************************
Printing Program Parameters:
programMode:			[ExecuteProgram]
fastaFileName:			/projects/bg/genomed/data/hg19_indices_ucsc/hg19.fa
alignFileName:			t_8.baf
algorithm:			[Best Score]
space:				[NT Space]
unpaired:			[Using]
reversePaired:			[Not Using]
avgMismatchQuality:		10
scoringMatrixFileName:		[Not Using]
randomBest:			[Not Using]
minMappingQuality:		-2147483648
minNormalizedScore:		-2147483648
pairingStandardDeviation:	2.000000
gappedPairingRescue		[Not Using]
numThreads:			4
queueLength:			50000
outputFormat:			[SAM]
outputID:			[Not Using]
RGFileName:			[Not Using]
timing:				[Not Using]
************************************************************
************************************************************
Reading in reference genome from /projects/bg/genomed/data/hg19_indices_ucsc/hg19.fa.nt.brg.
In total read 93 contigs for a total of 3137161264 bases
************************************************************
Postprocessing...
************************************************************
In function "AlignedEntryGetAlignment": Fatal Error[OutOfRange]. Message: Could not get reference sequence.
 ***** Exiting due to errors *****
************************************************************
I do not know what to make of this error - I ran a lot of analysis with the same reference sequence (in the same folder).

Suggestions are very much appreciated.

thanks and cheers,
Sophia
Please submit a test case to the bfast mailing list.
nilshomer is offline   Reply With Quote
Old 09-13-2011, 01:47 AM   #3
sdvie
Member
 
Location: Spain

Join Date: Jul 2010
Posts: 68
Default

Dear all and Nils,

when preparing the test set, I realized that the reference genome used for bfast postprocess was not the same I had used for alignment after all. (Too many editions of the hg19 around!) That is what caused the error.
However, as so often, to post the problem helped me to solve it .

Sorry for the confusion!

cheers,
Sophia
sdvie is offline   Reply With Quote
Old 09-13-2011, 06:25 PM   #4
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Excellent, thanks for the great insight.
nilshomer is offline   Reply With Quote
Old 03-02-2016, 10:13 AM   #5
xiangwulu
Member
 
Location: ireland

Join Date: Apr 2014
Posts: 18
Default

what does "avgMismatchQuality" mean?
does it refer to something with the base quality of mismach base?
xiangwulu is offline   Reply With Quote
Reply

Tags
bfast, error, outofrange, postprocess, reference

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 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO