bam(sam)格式文件转化为fasta格式
bam2fasta的转变方式:
samtools view SL3003_SL3004_100122-hg18.bam | \
awk '{OFS="\t"; print ">"$1"\n"$10}' - > SL3003.fasta
sam2fasta的转变方式:
cat *.sam | awk '{print ">"$1"\n"$10}' > *.fasta
samtools view SL3003_SL3004_100122-hg18.bam | \
awk '{OFS="\t"; print ">"$1"\n"$10}' - > SL3003.fasta
sam2fasta的转变方式:
cat *.sam | awk '{print ">"$1"\n"$10}' > *.fasta