Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 6f2a73cb authored by BOULLE Olivier's avatar BOULLE Olivier
Browse files

fixed expected length, removed start sequence

parent 67b55964
No related branches found
No related tags found
No related merge requests found
...@@ -179,7 +179,8 @@ then ...@@ -179,7 +179,8 @@ then
call_function python3 "$consensus_script" "$clusters_frag_dir" "$consensus_path" $coded_frag_length call_function python3 "$consensus_script" "$clusters_frag_dir" "$consensus_path" $coded_frag_length
else #ordered assembly else #ordered assembly
start_sequence=$(get_doc_param $container_path $document_index "start_sequence") #TODO remove start sequence -> use primers instead start_primer=$(get_doc_param $container_path $document_index "start_primer")
stop_primer=$(get_doc_param $container_path $document_index "stop_primer")
consensus_script="$project_dir"/sequencing_modules/reads_consensus.py consensus_script="$project_dir"/sequencing_modules/reads_consensus.py
...@@ -187,10 +188,10 @@ else #ordered assembly ...@@ -187,10 +188,10 @@ else #ordered assembly
then then
expected_length=$(($n_frag * $frag_length * 2)) # for a redundancy of 100% expected_length=$(($n_frag * $frag_length * 2)) # for a redundancy of 100%
else else
expected_length=$(($n_frag * $frag_length -66 )) # primers are not expected expected_length=$(($n_frag * $frag_length)) # primers are not expected but overhangs are (except on last fragments ->(-4))
fi fi
call_function python3 "$consensus_script" -i "$sequenced_reads_path" -o "$consensus_path" -s "$start_sequence" -e "$expected_length" call_function python3 "$consensus_script" -i "$sequenced_reads_path" -o "$consensus_path" --start $start_primer --stop $stop_primer -e $expected_length
fi fi
consensus_time=$(date +"%s") consensus_time=$(date +"%s")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment