Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error correcting raw reads #1

Closed
jagos01 opened this issue Apr 1, 2019 · 2 comments
Closed

Error correcting raw reads #1

jagos01 opened this issue Apr 1, 2019 · 2 comments

Comments

@jagos01
Copy link

jagos01 commented Apr 1, 2019

I am trying to correct some nanopore reads. My config file is called BaS_LR.txt and the following lings were modified:

PROJECT=BaS_LR
ONT_READ_LIST=BaS_LR_Seq.txt
GENOME_SIZE=5409120
THREADS=12
MIN_READ_LENGTH=500

BaS_LR.txt contains one line:
/home/data/MinION/HybridAssembly/BaS_S_L/Analysis/BaS_Guppy_merged_Porechop.fastq
which is a single fastq file containing all my reads (207,000).

The following is output when I run necat.pl correct BaS_LR.txt :

Smartmatch is experimental at /home/NECAT/Linux-amd64/bin/Plgd/Project.pm line 232.
2019-03-31 20:38:29 [Info] Start correcting rawreads.
2019-03-31 20:38:29 [Info] Start filtering rawreads for consensus.
2019-03-31 20:38:29 [Info] Run script: /home/data/MinION/HybridAssembly/BaS_S_L/Analysis/NECAT/BaS_LR/scripts/cns_pprr.sh 2>&1 |tee /home/data/MinION/HybridAssembly/BaS_S_L/Analysis/NECAT/BaS_LR/scripts/cns_pprr.sh.log
[Sun Mar 31 20:43:36 2019] ERROR: File '' does not exist (common/packed_db.c, 68)
[Sun Mar 31 20:38:29 2019] INFO: process fastq file /home/data/MinION/HybridAssembly/BaS_S_L/Analysis/BaS_Guppy_merged_Porechop.fastq (preprocess_raw_reads/main.c, 57)
2019-03-31 20:43:36 [Error] Failed to run script, 1, /home/data/MinION/HybridAssembly/BaS_S_L/Analysis/NECAT/BaS_LR/scripts/cns_pprr.sh

I am not sure how to correct these errors. Any help would be appreciated.
Thanks,
Scott

@xiaochuanle
Copy link
Owner

I guess your BaS_LR.txt contains two lines. But the second is an empty line. To validate this, please run the following command

$ wc -l BaS_LR.txt

If the output is something like

2 Bas_LR.txt

Then you should delete the second empty line.

@jagos01
Copy link
Author

jagos01 commented Apr 1, 2019

There was an empty line. Thank you for your help.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants