Commit a88025d2 authored by Mathieu Giraud's avatar Mathieu Giraud

vidjil.cpp: display number of reads and relative concentration in .fa headers of 'out/clones.fa'

parent 53defc7f
......@@ -1052,7 +1052,10 @@ int main (int argc, char **argv)
cout << representative.sequence << endl;
// ... and in out_clones
out_clones << ">clone-" << setfill('0') << setw(WIDTH_NB_CLONES) << num_clone << "-representative" << " " << seg.info << setfill(' ') << endl ;
out_clones << ">clone-" << setfill('0') << setw(WIDTH_NB_CLONES) << num_clone << "-representative"
<< "-" << setfill('0') << setw(WIDTH_NB_READS) << clone_nb_reads
<< "-" << setprecision(3) << 100 * (float) clone_nb_reads / nb_segmented << "%"
<< " " << seg.info << setfill(' ') << endl ;
out_clones << representative.sequence << endl;
break ;
......
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