Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • dineshkumarsrk
    Junior Member
    • Mar 2019
    • 5

    How to reverse complement the nuclotides, if the coordinates are in reverse order?

    I have the series of coordinates in id.txt file, whose coordinates sequences are in genome.fasta file. The coordinates of id.txt file is shown below,
    Contig3:15-7
    Contig2:5-10
    Contig1:12-3

    The genome.fasta file is shown below,
    >Contig1
    AAGGCCATCAAGGACGTGGATGAGGTCGTCAAG
    >Contig2
    ACGAAGCGCTCGCCAAGGCCGAAGAAGAAGGCC
    >Contig3
    GCTGCGGCGCTGATCCTGGCGGCCCGCGCCGAG

    I have used the following codes to extract the sequences from genome.fasta based on the coordinates in <id.txt> by using the following command <xargs samtools faidx genome.fasta <id.txt > result.fasta> It has extracted only <Contig2:5-10sequence>. Because, the coordinates of <Contig2:5-10> is in proper order. But <Contig3:15-7> and <Contig1:12-3> coordinates are in reverse order, So samtools could not fetch those sequences.

    I need to extract those sequences from <Contig2:5-10> and <Contig3:15-7> coordinates and also I want to reverse complement them. I have larger fasta file and plenty of coordinates sequences to be extracted, So I need to automate this process. So, please help me to automate this process. Thank you in advance.

Latest Articles

Collapse

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 06-09-2026, 11:58 AM
0 responses
22 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-05-2026, 10:09 AM
0 responses
29 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-04-2026, 08:59 AM
0 responses
39 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 12:03 PM
0 responses
61 views
0 reactions
Last Post SEQadmin2  
Working...