有两个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
人气:464 ℃ 时间:2020-03-21 14:30:54
解答
请给 精度 下个定义@ERR013180.1 HWI-EAS-249_38:2:1:2:857/1
TTTTCTTGTTCTTGACTCTTCTGCATAAGTANTTAAATCC
+
BBBBBCBB=BCBB6BBBB6BB>.6,9@1'6,!3B6'3,6?
第二行是断片第四行是对应第二行的精度转成ASCII码再减33 就是精度了! 谢谢
推荐
猜你喜欢
- 将4个数abcd排成2行2列将4个数a,b,c,d排成2行,2列,两边各加一条竖直线,记成 abcd ,定义 adcb =ad-bc ,
- 求名人孝敬长辈的故事一则
- 怎么理解冰醋酸、液氨、乙醇作为溶剂具有质子自递反应?
- 初一新学期开学感言
- 什么气候用蒸发浓缩冷却结晶,什么时候用趁热...
- 超市运来一批水果糖和巧克力糖,水果糖的颗数比巧克力糖的2倍还多40颗.售货员将5颗巧克力糖和12粒水果糖包成一个小包.最后巧克力糖全部装完,水果糖还剩8颗.问:水果糖共有几颗,巧克力糖共几颗?
- 无脊椎动物属于哪一级分类单位?
- 在口里填上合适的数.