導(dǎo)讀
fasta恬总、sam、bam相互轉(zhuǎn)換肚邢∫佳撸總結(jié)一下拭卿，可能不是最好用的方法。
一贱纠、fasta <=> sam
1. fasta => sam
# bowtie2比對獲得sam文件
bowtie2 -1 *_1.fastq -2 *_2.fastq -p 16 -x prefix -S *.sam
2. fasta <= sam
cat *.fasta | awk -F"\t" '{print ">"$1"\n"$10}' > *.sam
二峻厚、sam <=> bam
1. sam => bam
# -@：線程
# -b：輸出格式為BAM
# -S：自動檢測輸入格式
# -o：輸出文件
samtools view -@ 16 -b -S final.sam -o final.bam
2. sam <= bam
samtools view *.bam -O SAM > *.sam
三、fasta <=> bam
1. fasta => bam
# bam -> fq -> fa
bam2fastq --aligned *.bam -o *.fq
2. fasta <= bam
samtools view *.bam -O SAM | awk -F"\t" '{print ">"$1"\n"$10}' > *.fa
