摘要: 方法1 awk'{if(NR%4 == 1){print ">" substr($0, 2)}}{if(NR%4 == 2){print}}'fastq>fasta方法2: awk'BEGIN{P=1}{if(P==1||P==2){gsub(/^[@]/,">");print}; if(P==4)... 阅读全文
posted @ 2015-11-16 15:40 大雪无痕2 阅读(1348) 评论(0) 推荐(1) 编辑
摘要: bam(sam)格式文件转化为fasta格式bam2fasta的转变方式:samtools view SL3003_SL3004_100122-hg18.bam | \awk '{OFS="\t"; print ">"$1"\n"$10}' - > SL3003.fastasam2fasta的转变方... 阅读全文
posted @ 2015-11-16 15:05 大雪无痕2 阅读(1207) 评论(0) 推荐(0) 编辑
摘要: vi 文件名 #编辑方式查看,可修改cat 文件名 #显示全部文件内容more 文件名 #分页显示文件内容tail 文件名 #仅查看尾部,还可以指定行数head 文件名 #仅查看头部,还可以指定行数sed '1,$'p 文件名 #这个可以指定起始行和结束行的,把1和$换了就可以awk .......... 阅读全文
posted @ 2015-11-16 14:51 大雪无痕2 阅读(178) 评论(0) 推荐(0) 编辑