$ sh classification.sh SAMPLE_NAME reads.fa
$ sh concoct.sh SAMPLE_NAME contigs.fa PEreads1 PEreads2 species_number
$ sh reads_mapping.sh SAMPLE_NAME contigs_cutting.fa PEreads1 PEreads2
$ sh target_find.sh contigs_cutting.fa marker.fa clustering_gt1000.csv SAMPLE_NAME
$ sh blast_megahit.sh contigs_cutting.fa blast_result_file clustering_gt1000.csv out.F DIR_NAME
$ sh target_evaluation.sh clusterN(target cluster) contig_species cluster_sN(target cluster) out.F DIR_NAME
$ sh dataloss.sh contigs_list clustering_gt1000.csv out.F DIR_NAME
$ sh findclusterreads.sh cluster_sN(target cluster) out.F SAMPLE_NAME
$ python coverage_percent.py
Huazhong University of Science & Technology
If you have any problems while using this pipeline, please contact us at once and provide some details about problems, we will fix them as soon as possible.
You can follow us on Github , most of our softwares are opensource!
© 2017 QC-bind. All rights reserved | HUST Shawn