map.err in Soap denovo

I am using Soap denovo. I keep getting this error.

********************
Map
********************
Parameters: map -s /scratch/sysusers/godwin/Soap2/soap.contig -g graph_prefix
Kmer size: 23.
Contig length cutoff: 25.
Cannot open graph_prefix.contig. Now exit to system...

How do I change Kmer size then. or contig length cutoff?

assembly

• 1,691 views • link

updated 4.0 years ago by Ram

• written 11.1 years ago by mafireyi •

0 answers

No answers yet.

More posts like this

I do not understand how to get total peaks from MACS2 shown to me. What mistake(s) am I making here?
written 2.4 years ago by Ronin •

Hi all, I am trying to compare peak files between Homer and MACS2. My goal is to compare as integers the number of peaks between …
Cannot open file.fasta. Now exit to system...
written 2.8 years ago by san96 •

Hi, I'm trying to run SOAPdenovo but I get the following: user@DESKTOP-05K13O7:~/agave/crudos/ensamble$ SOAPdenovo-Trans-1.0.5/SOAPdenovo-Trans-31mer all -s config_file -o output_prefix Version 1.04 pregraph -s config_file -o output_prefix …
EDAseq, edgeR libraries cannot be installed in M2 mac
written 3.9 years ago by Sakib Salam •

I am having such a weird problem of not being able to install EDAseq or edgeR package from TCGAbiolinks in new M2 mac. Where I …
struct.error in MACS2 callpeak
written 5.2 years ago by ty0017 •

Hi all, I am trying to callpeaks on sorted and subsampled BAM files from the results of my ChIP-seq program. When I do callpeaks, I …
Error during Spades Meta-assembly: err code: -9
written 7.6 years ago by Ming •

Hello there, I was trying to run a meta-assembly on my Hi-Seq samples and I came across this error: Command line: /scratch/shiming/tools/SPAdes-3.13.0-Linux/bin/spades.py -o /gpfs0/scratch/shiming/trimmed_sequences/fastq/Assembly-Spades/CFC280618_metaspades --meta …
Why mapping rate is low after adapter trimming paired-end data?
written 9.5 years ago by gundalav

I have the pipeline with the following trimming and mapping parameter: trim_galore --stringency 5 --dont_gzip --trim1 --length 30 -q 0 -a CTGTCTCTTATACACATCT -o $TRIMGALORE_DIR $READ1 …
Running cufflinks - gtf file not recognised
written 10.1 years ago by snp87 •

Hi, I am trying to run cufflinks on my data with following command, cufflinks -p 4 -N –g Mus_musculus.GRCm38.84.gtf -b –u Mus_musculus.GRCm38.dna.toplevel.fa accepted_hits_1.bam but keep …
Error: segment-based junction search failed with err =127
written 10.8 years ago by leenaehyeon •

tophat2 mapping error I always use this command and mapping is successfully end. tophat -o tophat_out_sample -G hg19_genes.gtf --transcriptome-index hg19_transcripts --library-type fr-firststrand hg19 Sample_1_bbduk.fastq Sample_2_bbduk.fastq …
Allpaths- Keep getting ConvertToFastbQualb.pl failed for group 'paired_ends' when I run PrepareAllP…
written 11.1 years ago by mafireyi •

I am trying to use Allpaths for denovo assembly. My data summary looks like the following. ``` Hiseq_Run12_17122014 25GB Mate-pair (Size selected to 3KB) Hiscan_Run20_12022012 …
Perl program that parses the DNASIS restriction enzyme file
written 12.0 years ago by Sidney •

Code below is what I have so far. output not working at all. # Declare and initialize variables my %rebase_hash = ( ); my @file_data …