摘要:
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 阅读全文
摘要:
# 用于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/ 阅读全文