Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • surachit
    Junior Member
    • Oct 2015
    • 1

    #1

    Combining fasta file and bed file into a different file using biopython,python

    hi, I want to create a file using biopython or python or java .for each bed file range, the particular ATGC sequence need to be retrieved and all the result is to be put in this file.
    for example, bed file contains like this:

    chr1 14362 14829 chr1:14363-14829:WASH5P 467 +
    chr1 14969 15038 chr1:14970-15038:WASH5P 69 +
    chr1 15795 15947 chr1:15796-15947:WASH5P 152 +
    chr2 14362 14829 chr2:14363-14829:WASHP 467 +
    chr3 14969 15038 chr3:14970-15038:WASH 69 +
    chr10 15795 15947 chr10:15796-15947:WASHOP 152 +
    ..........................................................................................
    ........................................................................................

    and fasta file contains like this:

    >chr1 dna:chromosome chromosome:GRCh37:1:1:249250621:1
    NNNNNGCCAAGTnggggctaaNNNNGGGGCCCCCCCCCCCCCcCCC

    >chr2 dna:chromosome chromosome:GRCh37:1:1:249250621:1
    NNNNNGCCAAGNNNNGCCAAGT
    nggggctaaNNNNGCCAAGT
    nggggctaaNNNNGCCAAGT
    nggggctaa

    >chr3 dna:chromosome chromosome:GRCh37:1:1:249250621:1
    AGTACNNNNGCCAAGT
    nggggctaaNNNNGCCAAGT
    nggggctaa
    ................................................
    ...................................

    Now I want the result file like the following:

    chr1 14362 14829 chr1:14363-14829:WASH5P 467 +
    ATTTGCCCCCNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
    ATTTTTTTTTTTTTTTTTTTTTTTTTTTTTTGGGGGGGGGGGGGGGGGGGGG
    chr1 14969 15038 chr1:14970-15038:WASH5P 69 +
    NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNATTTTTTTTTT
    AAAAAAAAAAAATTTTTTTTTTTTTTTTTTTTTt

    itmeans for a particular range I have to retrieve the particular fasta sequence and put them into the file
    I did it in java but it is taking very long time...may be I am doing some mistake please help me such that I can use some library function and run the program easily

Latest Articles

Collapse

  • SEQadmin2
    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
    by SEQadmin2



    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
    ...
    10-14-2026, 09:55 AM
  • SEQadmin2
    New Genomics Technologies Take Aim at Long-Standing Limits
    by SEQadmin2


    Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

    We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
    ...
    09-28-2026, 10:25 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, Yesterday, 11:13 AM
0 responses
19 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, Yesterday, 10:48 AM
0 responses
18 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-15-2026, 10:52 AM
0 responses
36 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-14-2026, 05:39 AM
0 responses
42 views
0 reactions
Last Post SEQadmin2  
Working...