有两个文件,一个文件sam是一列read值id号;另外一个是一份fastq文件。
想要写一段脚本(python或者perl都可以) 。根据文件sam中的id号,从
fastq文件中提取出相应得read值(完整的read)写在另外一个新文件上
http://www.360doc.com/content/16/0322/02/26456292_544193476.shtml
直接拿id去取不就可以了?
或者把fastq文件的数据放到数据库(sqlite或者其他数据库)里面,取起来不是更方便了?