有两个fastaq格式的DNA序列文件,想写一个perl程序完成!
问题描述：
有两个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
答
请给 精度 下个定义@ERR013180.1 HWI-EAS-249_38:2:1:2:857/1
TTTTCTTGTTCTTGACTCTTCTGCATAAGTANTTAAATCC
+
BBBBBCBB=BCBB6BBBB6BB>.6,9@1'6,!3B6'3,6?
第二行是断片第四行是对应第二行的精度转成ASCII码再减33 就是精度了！ 谢谢