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

Unknown error in raw_reads_filter_v0.5.pl #43

Closed
srimku opened this issue Jan 21, 2020 · 8 comments
Closed

Unknown error in raw_reads_filter_v0.5.pl #43

srimku opened this issue Jan 21, 2020 · 8 comments

Comments

@srimku
Copy link

srimku commented Jan 21, 2020

I am getting this error for my initial run using MitoZ. I request you to kindly help me to run MitoZ successfully.


python3 MitoZ.py all --genetic_code 1 --clade Chordata --thread_number 13 --outprefix duck-01 --fastq1 TS-COR_R1.fq.gz --fastq2 TS-COR_R2.fq.gz --fastq_read_length 300 --insert_size 300 --run_mode 2 --filter_taxa_method 1 --requiring_taxa 'Chordata'
MitoZ.py:1942: DeprecationWarning: time.clock has been deprecated in Python 3.3 and will be removed from Python 3.8: use time.perf_counter or time.process_time instead
program_begin_clock = time.clock()
2020-01-21 16:59:57
perl /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl -1 /home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R1.fq.gz -2 /home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R2.fq.gz -3 /home/sridhar/mito/release_MitoZ_v2.3/tmp/duck-01.cleandata/clean.1.fq.gz -4 /home/sridhar/mito/release_MitoZ_v2.3/tmp/duck-01.cleandata/clean.2.fq.gz -m 3 -l 15 -n 10 -q 55,20 -z

input files:
/home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R1.fq.gz
/home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R2.fq.gz
output files:
/home/sridhar/mito/release_MitoZ_v2.3/tmp/duck-01.cleandata/clean.1.fq.gz
/home/sridhar/mito/release_MitoZ_v2.3/tmp/duck-01.cleandata/clean.2.fq.gz
filter parameters:
k, keep only bases between BEG and END for each read: full length read
n, the maximun N allowed in single end reads: 10
q, the maixmum percentage of low quality (ASCII code's integer) bases allowed in single end reads: 55,20
d, Filtering duplications: False

Use of uninitialized value $r01 in chomp at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 181, line 25348.
Use of uninitialized value $r02 in chomp at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 182, line 25348.
Use of uninitialized value $r03 in chomp at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 183, line 25348.
Use of uninitialized value $r04 in chomp at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 184, line 25348.
Use of uninitialized value $r01 in exists at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 195, line 25348.
Use of uninitialized value $r02 in transliteration (tr///) at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 199, line 25348.
Use of uninitialized value $seq_len in subtraction (-) at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 261, line 25348.
Use of uninitialized value $seq_len in division (/) at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 265, line 25348.
Illegal division by zero at /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl line 265, line 25348.
Error occured when running command:
perl /home/sridhar/mito/release_MitoZ_v2.3/bin/filter/raw_reads_filter_v0.5.pl -1 /home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R1.fq.gz -2 /home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R2.fq.gz -3 /home/sridhar/mito/release_MitoZ_v2.3/tmp/duck-01.cleandata/clean.1.fq.gz -4 /home/sridhar/mito/release_MitoZ_v2.3/tmp/duck-01.cleandata/clean.2.fq.gz -m 3 -l 15 -n 10 -q 55,20 -z


@srimku
Copy link
Author

srimku commented Jan 21, 2020

I am getting error with raw_reads_filter_v05.pl

Please give your suggestions to rectify this error

@linzhi2013
Copy link
Owner

linzhi2013 commented Jan 21, 2020

Please give a few examples of your input fastq files. I can not determine what was going on there based on the message you posted.


gzip -dc /home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R1.fq.gz | head -n 20 > test.1.fq
gzip -dc /home/sridhar/mito/release_MitoZ_v2.3/TS-COR_R2.fq.gz | head -n 20 > test.2.fq

Then send me the test.1.fq and test.2.fq files.

PS.

  1. The mitochondrial genetic code for vertebrate is 2, not 1.
  2. If you prefer to install MitoZ from source code, please run the test data first (https://github.com/linzhi2013/MitoZ/tree/master/test), to make sure everything is fine. Otherwise, please use Docker version.

@srimku
Copy link
Author

srimku commented Jan 22, 2020 via email

@linzhi2013
Copy link
Owner

Hi Dr. Srimku,

There was a typo in your command.

The correct one is guanliangmeng/mitoz:2.3

Not guangliangmeng/mitoz:2.3

@srimku
Copy link
Author

srimku commented Jan 22, 2020 via email

@srimku
Copy link
Author

srimku commented Jan 22, 2020 via email

@srimku
Copy link
Author

srimku commented Jan 22, 2020 via email

@linzhi2013
Copy link
Owner

linzhi2013 commented Jan 22, 2020 via email

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