从bam文件中提取fasta文件
2021-03-25 本文已影响0人
白米饭睡不醒
samtools view GPR210305-17-VI.Hisat_aln.sorted.bam | awk '{OFS="\t"; print ">"$1"\n"$10}' - > output.fasta
samtools view GPR210305-17-VI.Hisat_aln.sorted.bam | awk '{OFS="\t"; print ">"$1"\n"$10}' - > output.fasta