yeastgenome / AGAPE

AGAPE (Automated Genome Analysis PipelinE) for yeast pan-genome analysis
3 stars 9 forks source link

find_inserted_intervals large negative values #4

Open kastman opened 6 years ago

kastman commented 6 years ago

Occasionally with some alignments (currently in the non_ref_contigs.sh workflow) I'm getting very large negative numbers for the beginning of my intervals; e.g. https://github.com/yeastgenome/AGAPE/blob/master/non_ref_contigs.sh#L21

$BIN/find_inserted_intervals $out_dir/temp.maf $scf_name $len >> $out_dir/$seq_name.inserted.intervals

occasionally creates intervals like this (see unplaced-504 beginning at -536870911):

$ cat $seq_name.inserted.intervals | head 

unplaced-402 2 618
unplaced-423 485 1991
unplaced-504 -536870911 206
unplaced-529 1 1724
unplaced-563 1 1521

Any idea why this might be, and how to fix it?