xiaochuanle / MECAT

MECAT: an ultra-fast mapping, error correction and de novo assembly tool for single-molecule sequencing reads
107 stars 26 forks source link

problem with step 3 extract_sequence,no results file... #32

Open fhu2008-tio opened 6 years ago

fhu2008-tio commented 6 years ago

Only success once after install the soft, (step1 to step4 are all ok!)

but when repeat again ,fail at step 3 with no results file generated,

command: extract_sequences ecoli_filtered.fastq.corrected.fasta corrected_ecoli_25x.fasta 4800000 25)

no corrected_ecoli_25x.fasta.fasta.qual , no corrected_ecoli_25x.fasta.fasta.qv no corrected_ecoli_25x.fasta.fasta.frg no corrected_ecoli_25x.fasta.fasta

also, there were no errors, no warnings,as show below?

does anyone else can fix it ?

[root@localhost Genome_Assembly]# extract_sequences ecoli_filtered.fastq.corrected.fasta corrected_ecoli_25x.fasta 4800000 25

step 1: convert fasta to fastq step 2: convert fastq to CA step 3

Starting file 'ecoli_filtered.fastq.corrected.fasta.fastq.libname.frg'.

Processing SINGLE-ENDED SANGER QV encoding reads from: '/home/workplace/Genome_Assembly/ecoli_filtered.fastq.corrected.fasta.fastq'

GKP finished with no alerts or errors.

step 4 Longest picked cutoff: 12691 Scanning store to find libraries used. Added 0 reads to maintain mate relationships. Dumping 0 fragments from unknown library (version 1 has these) Dumping 7639 fragments from library IID 1 step 5 Longest picked cutoff: 12691 Scanning store to find libraries used and reads to dump. Added 0 reads to maintain mate relationships. Dumping 0 fragments from unknown library (version 1 has these) Dumping 7639 fragments from library IID 1

xiaochuanle commented 6 years ago

@fhu2008-tio This is the problem with the name of the file suffix, we have fixed this problem in the updated mecat. Thanks for using mecat!