This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extract codon between two flanking sequences from FASTA

Hi all, given two 10bp flanking sequences, I would like to extract the codon between them across all samples in a FASTA file. For example:

Flank1:
CAGGCATGCC
Flank2:
TCATCGCTGG

FASTA 

>sample1
GCGCACCATGGTCAGGCATGCCTCCTCATCGCTGGGCACAGCCCAGAGGGT
>sample2
GGCAGAACCCGCGCACCATGGTCAGGCATGCCACCTCATCGCTGGGCACAGCCCAGA
>sample3
GGCAGATTCCCCGCACCATGGTCAGGCATGCCACTTCATCGCTGGGCACA

Output

>sample1
TCC
>sample2
ACC
>sample3
ACT

I have performed the opposite of this where I extract the flanking sequences based on coordinates, as well as extracting the sequences between two coordinates using bedtools getfasta, but struggling with extracting based on flanking nucleotide sequences. Thanks for any help!

genome sequencing
• 1,356 views • link
written by renyulb • renyulb

in perl

$seq =~ /CAGGCATGCC ([ATGC]{3})+ TCATCGCTGG/x

print $1
• 0 views • link
written by Sishuo Wang Sishuo Wang

This can work if the sequence is not wrapped.

perl -n -e '{if($_ =~ /CAGGCATGCC(.+)TCATCGCTGG/){print $x, $1, "\n" };$x=$_}' Test.fa

• 0 views • link
written by microfuge microfuge

0 answers

No answers yet.

Log in to answer this question.

More posts like this