摘要: then switch to the seqtk directory and make it: Now you can invoke it as such: Convert FASTQ to FASTA: seqtk seq -a in.fq.gz > out.fa Convert ILLUMINA 阅读全文
posted @ 2017-01-19 13:17 Datapotumas 阅读(4719) 评论(0) 推荐(0) 编辑
摘要: # 用于fasta格式文件的碱基数目和GC含量的统计 grep -v '>' input.fa| perl -ne '{$count_A=$count_A+($_=~tr/A//);$count_T=$count_T+($_=~tr/T//);$count_G=$count_G+($_=~tr/G/ 阅读全文
posted @ 2017-01-19 13:14 Datapotumas 阅读(7329) 评论(1) 推荐(1) 编辑