生物信息_Call_snp_by_soapsnp_外顯子
數據：
人外顯子（5.6G左右，兩個gz文件，Agilent 50M，pe測序，90讀長）
方法：call_snp_by_soapsnp
每步估計需要投多大（僅作參考）：
SOAP（也就是Soapaligner，5.6G）：
從兩個fq.gz到分染色體之後的.soap
注：SOAP這步輸入文件有一個chrOrder，要自己生成的，一般簡單複製粘貼略作修改就好。中途會生成.soap.pe、.soap.se、.soap.unmap，把.soap.pe、.soap.se一塊cat到.soap.pe.se.gz，再對.soap.pe.se.gz分染色體就會得到.soap。最好建一個文件夾放.soap結果。
SOAPsnp（0.5G）:
從各個染色體的.soap到.snp
注：.snp結果實際上是.cns文件