Originally posted by NGS group
View Post
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I'm always getting the error below. Tried bwa and bowtie, same error. Also different input files. Can anybody help?
Thanks,
S.
=>Thu Jun 26 13:58:40 2014: Building Bowtie index for contigs
Warning: Empty input file
Reference file does not seem to be a FASTA file
Command: /x/bin/SSPACE-STANDARD-3.0_linux-x86_64/bowtie/bowtie-build --quiet --noref 13.18a/tmp.13.18a/subset_contigs.fasta 13.18a/alignoutput/13.18a.13.18.bowtieIndex
Bowtie-build error; 256 at /x/bin/SSPACE-STANDARD-3.0_linux-x86_64/bin/PairingAndScaffolding.pl line 80.
**************************************************
Process 'mapping reads' failed on Thu Jun 26 13:58:40 2014
#the 'summary' file looks like this:
READING READS 13.18:
------------------------------------------------------------
Total inserted pairs = 2200000
------------------------------------------------------------
CONTIG EXTENSION:
------------------------------------------------------------
Number of contig sequences =139
Number of unmapped reads used for contig extension = 1009817
Number of contigs extended = 110
Number of bases extended = 25878
------------------------------------------------------------
LIBRARY 13.18 STATS:
################################################################################Last edited by susanklein; 06-25-2014, 09:32 PM.
Comment
-
Hello guys,
I'm trying to extend contigs / join them into scaffolds but at some point I get this error message:
Code:Process 'extend/format contigs' failed
I'm running sspace with -v option but even so I'm not getting extra error output. Has any of you got through this?
This is the command I'm using to call sspace:
Code:perl /home/george/SSPACE-STANDARD-3.0_linux-x86_64/SSPACE_Standard_v3.0.pl -x 1 -k 5 -z 100 -a 0.7 -m 35 -o 20 -b snake_extension -l libraries.txt -T 80 -s /work/snake/montagem/abyss/20140421_k91/ snake-6.fa -v 1>output.log 2>output.err
Thanks in advance fellows!
Condomitti.
Comment
-
Does this error mean that no reads were unmaped an therefore the file was empty??
Thanks, anyone??
Originally posted by susanklein View PostI'm always getting the error below. Tried bwa and bowtie, same error. Also different input files. Can anybody help?
Thanks,
S.
=>Thu Jun 26 13:58:40 2014: Building Bowtie index for contigs
Warning: Empty input file
Reference file does not seem to be a FASTA file
Command: /x/bin/SSPACE-STANDARD-3.0_linux-x86_64/bowtie/bowtie-build --quiet --noref 13.18a/tmp.13.18a/subset_contigs.fasta 13.18a/alignoutput/13.18a.13.18.bowtieIndex
Bowtie-build error; 256 at /x/bin/SSPACE-STANDARD-3.0_linux-x86_64/bin/PairingAndScaffolding.pl line 80.
**************************************************
Process 'mapping reads' failed on Thu Jun 26 13:58:40 2014
#the 'summary' file looks like this:
READING READS 13.18:
------------------------------------------------------------
Total inserted pairs = 2200000
------------------------------------------------------------
CONTIG EXTENSION:
------------------------------------------------------------
Number of contig sequences =139
Number of unmapped reads used for contig extension = 1009817
Number of contigs extended = 110
Number of bases extended = 25878
------------------------------------------------------------
LIBRARY 13.18 STATS:
################################################################################
Comment
-
Originally posted by condomitti View PostHello guys,
I'm trying to extend contigs / join them into scaffolds but at some point I get this error message:
Code:Process 'extend/format contigs' failed
I'm running sspace with -v option but even so I'm not getting extra error output. Has any of you got through this?
This is the command I'm using to call sspace:
Code:perl /home/george/SSPACE-STANDARD-3.0_linux-x86_64/SSPACE_Standard_v3.0.pl -x 1 -k 5 -z 100 -a 0.7 -m 35 -o 20 -b snake_extension -l libraries.txt -T 80 -s /work/snake/montagem/abyss/20140421_k91/ snake-6.fa -v 1>output.log 2>output.err
Thanks in advance fellows!
Condomitti.
Any clue here? Anyone who could have faced something like that?
Thanks,
Condomitti.
Comment
-
Originally posted by boetsie View PostHmmm, It seems that in the newest version of perl they removed the getopts.pl library (see http://search.cpan.org/~rjbs/perl-5....s_and_Pragmata). At this site they explain how to solve this issue:
Hope this helps.
Boetsie
I looked at the nasa site and I am confused about how to solve the getopts.pl problem.
I found this on anothe page at the nasa site. The link to the fix it script is broken though....
Perl scripts with newer Perl:
"Many of the older Perl scripts in HEASoft (e.g. fhelp) require the getopt.pl or getopts.pl libraries which are no longer available in the newest Perl distributions, e.g. Perl 5.16.2 which ships with Fedora 18. This can result in the error "Can't locate getopt.pl in @INC..." when trying to run a script. To get around this, users may either download and install the Perl4::CoreLibs module from CPAN which includes getopt.pl, or download and run (in $HEADAS/../ftools/<architecture>/bin) a script which repairs our older Perl scripts to use Getopt::Std instead."
Nathan
Comment
-
Originally posted by maubp View PostNathan's post just above outlines two ways to fix this, the simpler one sounded like installing Perl4::CoreLibs module from CPAN.
Comment
-
I now have another problem with the ExtendOrFormatContigs.pl script. Attempts to use bowtie failed due to problems with the output of the aligner and so I tried bwa as the aligner. However now after alignement I get this error call from the ExtendOrFormatContigs.pl script.
Bwa error; 0 at ExtendOrFormatContigs.pl line 288.
this is the line
system("$aligninput") == 0 || die "\nBwa error; $?" if($aligninput ne "");
Not sure what this means. My alignmentment folder looks fine.
I'm trying again without the extend contigs option but I sure would like to extend my contigs before scaffolding...
NathanLast edited by fireant; 08-07-2014, 10:56 AM.
Comment
-
Hi,
Has anyone seen this problem, and know how to fix it. We were going to remove the offending reads, but the we don't know how to convert the read names in the Reads dir, back to the input read names for removal. Or is there something else we can do? These are BAC end sequences - 3730 reads.
Thanks,
Pat
Error: Read (read1ad27Error: Read (Error: Read (rError: Read (read891/1)
is less than 4 characters long
203/Error: Read (read74read815/2) is less than 3 characters long
2) is less than 4 characters long
8/2) is less than 3 characters long
ead1170/2) is less than 3 characters long
73/2) is less than 3 characters long07/2
) is less than 3 characters long
Error: Read (read1576/2) is less than 3 characters long
terminate called recursively
terminate called recursively
Error: Read (Error: Read (ead273/2) is less than 4 characters long
Error: Read (eadError: Read (read748/2) is less than 3 characters long
read207/2) is less than 3 characters long
295/2) is less than 3 characters long
Error: Read (read815/2) is less than 3 characters long
Error: Read (read891/1) is less than 4 characters long
Error: Read (read1073/2) is less than 3 characters long
Error: Read (read1170/2) is less than 3 characters long
terminate called after throwing an instance of 'terminate called
recursively
int'
Comment
-
Originally posted by fireant View PostBwa error; 0 at ExtendOrFormatContigs.pl line 288.
this is the line
system("$aligninput") == 0 || die "\nBwa error; $?" if($aligninput ne "");
Not sure what this means. My alignmentment folder looks fine.
I'm trying again without the extend contigs option but I sure would like to extend my contigs before scaffolding...
Nathan
Cheers
Comment
-
Greetings all! Another plea for wisdom from SSPACE users...
So, I am working on a 1.7Gb diploid de novo genome assembly (http://seqanswers.com/forums/showthread.php?t=42555) and wanted to try and employ some enhanced scaffolding to improve continuity of some of my existing draft assemblies. However, I get an error message upon completion of the sspace run:
"Process 'extend/format contigs' failed on Sat Dec 13 07:49:06 2014"
I have accepted default values to start with on the majority of the settings:
perl /home/jpummil/SSPACE-STANDARD-3.0_linux-x86_64/SSPACE_Standard_v3.0.pl -l libraries.txt -s Ray-53-contigs.fasta -x 1 -T 32 -b Ray_standard_out_v2_extend
Though sizeable (400GB memory and 18 hour runtime on 32 compute cores), not hitting any memory barriers, etc...job does complete and exit properly.
Think my libraries.txt is OK...
[jpummil@razor-l3 Snake_v2]$ more libraries.txt
lib1 bowtie /scratch/jpummil/Snake_v2/L001_R1_MiSeq.fastq /scratch/jpummil/Snake_v2/L001_R2_MiSeq.fastq 479 0.25 FR
lib2 bowtie /scratch/jpummil/Snake_v2/s1_R1_PE_Phred33.fastq /scratch/jpummil/Snake_v2/s1_R2_PE_Phred33.fastq 150 0.25 FR
lib3 bowtie /scratch/jpummil/Snake_v2/s2_R1_PE_Phred33.fastq /scratch/jpummil/Snake_v2/s2_R2_PE_Phred33.fastq 150 0.25 FR
Anything obvious come to mind? I had understood that with -x 1 flag, the assembly would be physically modified during the scaffolding process?
Thanks for any tips, and HAPPY HOLIDAYS!
Comment
-
Originally posted by jpummil View PostThink my libraries.txt is OK...
[jpummil@razor-l3 Snake_v2]$ more libraries.txt
lib1 bowtie /scratch/jpummil/Snake_v2/L001_R1_MiSeq.fastq /scratch/jpummil/Snake_v2/L001_R2_MiSeq.fastq 479 0.25 FR
lib2 bowtie /scratch/jpummil/Snake_v2/s1_R1_PE_Phred33.fastq /scratch/jpummil/Snake_v2/s1_R2_PE_Phred33.fastq 150 0.25 FR
lib3 bowtie /scratch/jpummil/Snake_v2/s2_R1_PE_Phred33.fastq /scratch/jpummil/Snake_v2/s2_R2_PE_Phred33.fastq 150 0.25 FR
So a library file would have been something like:
Code:Org1 bwa Truseq_1_LRTrimmed.fastq Truseq_2_LRTrimmed.fastq 430 0.50 FR
Last edited by rhinoceros; 12-17-2014, 11:46 PM.savetherhino.org
Comment
-
Originally posted by rhinoceros View PostIs your insert size really just 150 ± 38 bp? 150 bp sounds more like your read length. As far as I recall, I never got SSPACE working with bowtie. However, it worked just fine with BWA..
So a library file would have been something like:
Code:Org1 bwa Truseq_1_LRTrimmed.fastq Truseq_2_LRTrimmed.fastq 430 0.50 FR
Current Filename File Size Encoding Length Insert Size Type
L001_R1_MiSeq.fastq 16GB Illumina 1.9 300bp 479 Pair Ended
L001_R2_MiSeq.fastq 16GB Illumina 1.9 300bp 479 Pair Ended
L003_R1_MP.fastq 21GB Illumina 1.9 101bp 6000 Mate Pairs
L003_R2_MP.fastq 21GB Illumina 1.9 101bp 6000 Mate Pairs
s1_R1_PE.fastq 27GB Illumina 1.5 101bp 150 Pair Ended
s1_R2_PE.fastq 27GB Illumina 1.5 101bp 150 Pair Ended
s2_R1_PE.fastq 24GB Illumina 1.5 101bp 150 Pair Ended
s2_R2_PE.fastq 24GB Illumina 1.5 101bp 150 Pair Ended
I'll try a couple with BWA and see what happens.
Comment
Latest Articles
Collapse
-
by seqadmin
Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...-
Channel: Articles
10-18-2024, 07:11 AM -
-
by seqadmin
Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.
Nobel Prize for MicroRNA Discovery
This week,...-
Channel: Articles
10-07-2024, 08:07 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
New Model Aims to Explain Polygenic Diseases by Connecting Genomic Mutations and Regulatory Networks
by seqadmin
Started by seqadmin, Yesterday, 05:31 AM
|
0 responses
10 views
0 likes
|
Last Post
by seqadmin
Yesterday, 05:31 AM
|
||
Started by seqadmin, 10-24-2024, 06:58 AM
|
0 responses
20 views
0 likes
|
Last Post
by seqadmin
10-24-2024, 06:58 AM
|
||
New AI Model Designs Synthetic DNA Switches for Targeted Gene Expression in Specific Cell Types
by seqadmin
Started by seqadmin, 10-23-2024, 08:43 AM
|
0 responses
48 views
0 likes
|
Last Post
by seqadmin
10-23-2024, 08:43 AM
|
||
Started by seqadmin, 10-17-2024, 07:29 AM
|
0 responses
58 views
0 likes
|
Last Post
by seqadmin
10-17-2024, 07:29 AM
|
Comment