trim.seqs(fasta=CN.fa,qfile=CN.quala, maxambig=0, maxhomop=8, flip=T, bdiffs=1, pdiffs=2, qwindowaverage=35, qwindowsize=50, processors=7) system(./grupp_CN.pl CN.trim.fasta > CN.groups) unique.seqs(fasta=CN.trim.fasta) align.seqs(fasta=CN.trim.unique.fasta, reference=gg.ref.fasta, processors=7,flip=T) summary.seqs(fasta=CN.trim.unique.align) screen.seqs(fasta=CN.trim.unique.align, name=CN.trim.names, group=CN.groups, end=4951,start=4655,minlength=65,processors=7) summary.seqs(fasta=current) filter.seqs(fasta=CN.trim.unique.good.align, vertical=T, trump=., processors=7) unique.seqs(fasta=CN.trim.unique.good.filter.fasta, name=CN.trim.good.names) pre.cluster(fasta=CN.trim.unique.good.filter.unique.fasta, name=CN.trim.unique.good.filter.names, group=CN.good.groups, diffs=1) chimera.uchime(fasta=CN.trim.unique.good.filter.unique.precluster.fasta, name=CN.trim.unique.good.filter.unique.precluster.names, group=CN.good.groups, processors=7) remove.seqs(accnos=CN.trim.unique.good.filter.unique.precluster.uchime.accnos, fasta=CN.trim.unique.good.filter.unique.precluster.fasta, name=CN.trim.unique.good.filter.unique.precluster.names, group=CN.good.groups) system(mv CN.trim.unique.good.filter.unique.precluster.pick.names CN.final.names) system(mv CN.trim.unique.good.filter.unique.precluster.pick.fasta CN.final.fasta) system(mv CN.good.pick.groups CN.final.groups) system(./gr_puhas CN.final.groups > CN.proovikaupa.groups) remove.groups(fasta=CN.final.fasta, […]
Recent Comments