• python 中实现将fastq文件转换为fasta文件


    001、

    root@PC1:/home/test# ls
    a.fastq  test.py
    root@PC1:/home/test# cat a.fastq                                     ## 测试fastq文件
    @DJB775P1:248:D0MDGACXX:7:1202:12362:49613
    TGCTTACTCTGCGTTGATACCACTGCTTAGATCGGAAGAGCACACGTCTGAA
    +
    JJJJJIIJJJJJJHIHHHGHFFFFFFCEEEEEDBD?DDDDDDBDDDABDDCA
    @DJB775P1:248:D0MDGACXX:7:1202:12782:49716
    CTCTGCGTTGATACCACTGCTTACTCTGCGTTGATACCACTGCTTAGATCGG
    +
    IIIIIIIIIIIIIIIHHHHHHFFFFFFEECCCCBCECCCCCCCCCCCCCCCC
    root@PC1:/home/test# cat test.py                                    ## 测试程序
    #!/usr/bin/python
    in_file = open("a.fastq", "r")
    out_file = open("result.txt", "w")
    dict1 = {}
    idx = 0
    
    for i in in_file:
        idx = idx + 1
        i = i.strip()
        if idx % 4 == 1:
            key = "".join(i.split("@")[1:])
            dict1[key] = ""
        elif idx % 4 == 2:
            dict1[key] = i
    for i,j in dict1.items():
        print(">{0}\n{1}".format(i,j), file = out_file)
    
    in_file.close()
    out_file.close()
    root@PC1:/home/test# python test.py                                         ## 执行程序
    root@PC1:/home/test# ls
    a.fastq  result.txt  test.py
    root@PC1:/home/test# cat result.txt                                         ## 执行程序结果
    >DJB775P1:248:D0MDGACXX:7:1202:12362:49613
    TGCTTACTCTGCGTTGATACCACTGCTTAGATCGGAAGAGCACACGTCTGAA
    >DJB775P1:248:D0MDGACXX:7:1202:12782:49716
    CTCTGCGTTGATACCACTGCTTACTCTGCGTTGATACCACTGCTTAGATCGG

    参考:https://www.jianshu.com/p/5ee54bea4cb0

  • 相关阅读:
    printf里的=、++
    线程也疯狂-----异步编程
    自己搭建node服务器环境(请求静态资源、get请求、post请求)
    React学习
    2020.10-2021-01总结
    接圈的作用和缺点
    CWnd,HWND; CDC,HDC
    Python通过requests模块处理form-data请求格式
    element-ui resetFields 无效的问题
    用python 将数字每三组分割
  • 原文地址:https://www.cnblogs.com/liujiaxin2018/p/16588749.html
Copyright © 2020-2023  润新知