有两个fastaq格式的DNA序列文件,想写一个perl程序完成!
有两个fastaq格式的DNA序列文件,想同时把每个文件中每组断片的两端精度小于10的基因删除,之后再将每组平均精度低于30的断片删除,想写一个perl程序完成,无奈写不出,大体是这样的断片
@ERR013180.1 HWI-EAS-249_38:2:1:2:857/1
TTTTCTTGTTCTTGACTCTTCTGCATAAGTANTTAAATCC
+
BBBBBCBB=BCBB6BBBB6BB>.6,9@1'6,3B6'3,
@ERR013180.2 HWI-EAS-249_38:2:1:2:474/1
AACTGGCATTCACCCAAACCCCACCGACACCNGGTTTTAT
+
B7-=;;BB27?BBC@.33BBBBBB