Commit 646f8aec authored by khalid's avatar khalid
Browse files

Force the order of transcripts in /tmp/ALL_transcripts.ctgs

parent d524750f
......@@ -115,6 +115,9 @@ echo
esac
echo -e "\nIndex built : $out/$toolidx\n"
refseqALL=$out/ALL_transcripts.ctgs
echo '' > $refseqALL
#for (( j=0; j <i; j++ ))
for ref in "${orders[@]}"
do
......@@ -128,10 +131,11 @@ echo
Seqlength = length($2) ;
print a[1]"\t"Seqlength > refseqs
}' $out/$ref".fasta_mod"
# regroup all samples
cat $out/*.ctgs >> $refseqALL
done
# regroup all samples
refseqALL=$out/ALL_transcripts.ctgs
cat $out/*.ctgs > $refseqALL
# mapping successively every read sets on all transcripts.
#for (( k=0; k <i; k++ ))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment