learniup

lobstr_index.py can't build correctly index

1

ERROR: command 'intersectBed -a /tmp/tmpeUwvrC/lobSTR_ref_merged.bed -b /tmp/tmpeUwvrC/lobSTR_ref_plusflank_sorted.bed -wa -wb | awk '{print $1 "\t" $2 "\t" $3 "\t" $7"_"$8"_"$9";"}' | bedtools groupby -g 1,2,3 -c 4 -o concat > index_dir/lobSTR_mergedref.bed'...

[lsx@macaca lobSTR_data]$ python scripts/lobstr_index.py --str lobSTR_region --ref CR_genome_chr_edition.fa --out output_directory Traceback (most recent call last): File "scripts/lobstr_index.py", line 216, in if not DEBUG: GetRefFasta(genome, refkeys, merged_str_file, str_ref_fasta, str_map_file) File "scripts/lobstr_index.py",...