PacificBiosciences / FALCON_unzip

Making diploid assembly becomes common practice for genomic study
BSD 3-Clause Clear License
30 stars 18 forks source link

task_hasm() dies on missing ctgs #28

Closed flowers9 closed 7 years ago

flowers9 commented 8 years ago

Perhaps this is proper behavior, but in the deduping loop, if a ctg directory doesn't exist, the script exits:

for f incat ../reads/ctg_list ;do cd $WD/$f; fc_dedup_h_tigs.py $f; done

To get past this, I simply rewrote it as

for f in `cat ../reads/ctg_list`;do
if [ -d "$WD/$f" ]; then
cd $WD/$f
fc_dedup_h_tigs.py $f
fi
done

I'm not sure what a missing ctg directory implies - for example, if that's a critical error of some sort, and needs to be dealt with. If not, this patch a least allows the run to complete.

jfass commented 8 years ago

Check the "Possible MUMmer issue ..." issue ... doesn't answer your question, but there's other discussion of workarounds.

yingzhang121 commented 7 years ago

I had the same issue. So I just manually created the directory for all the names in ctg_list. It seems this way can push fcunzip.py forward. However, I do suspect whether the missing contig will be correctly incorporated into the multiple all* output files.