前情提要：

是这样的，服务器里的raid10阵列有一块一直识别不了，所以通过热拔插确定了是那一块之后（不要学习这个操作），重启时进入了initramfs。

以BRCA1为例，分别从三大数据库出发写几个低通量的方法。

实际上是抄来的，然后加了点如果不满足3的倍数时候报Warning。

得到一个需求，截断fastq里的引物序列，为什么是截断而不是去除，目的是为了交付数据时不被同行得到引物同时读长看起来也没那么短。一开始想象中是很简单的，只要替换fastq中的引物序列就好了，但是后来发现有些目标区域其实也含有引物的序列，如果直接替换了会导致这些目标区域变成缺失之类的判读。