pb.panigrahi86 +1 for that. I love one liners!
• 1 views
•
link
written
by
lakhujanivijay
Script for replacing selenocysteine in a phylip file
Hi, I need help. I do not have too much experience in coding. I am trying to replace selenocysteine (U) in a phylip file with X. How can I do that? I tried to use Bioperl but still not successful. My trial is written below.
#!/usr/bin/perl
use Getopt::Std;
use Bio::Seq;
use Bio::SeqIO;
use Bio::AlignIO;
use Bio::SimpleAlign;
my %opts = ();
getopts ('f:', \%opts);
my $file = $opts{'f'};
my $alnin = Bio::AlignIO->new (-format=>'phylip', -file=>"$file");
my $alnout = Bio::AlignIO->new (-format=>'phylip', -file=>"$file");
while (my $aln = $alnin->next_aln()){
my $id = $aln->id_linebreak();
my $seq = $aln->seq();
$seq =~s/U/X/g;
$seq =~s/O/X/g;
print "$id\n$seq\n";
}
• 2,753 views
•
link
updated
by
finswimmer
•
written
by
ahmedmagds •
•
written
by
ahmedmagds •
2 answers
Since you have not given any input file, but mentioned the type to be phylip. I took an example phylip format file from here and try to substitute U-> X and O->X randomly and tried the code below.
cat phylip.txt | perl -ne 'my @x=split /( +)/,$_; print $x[0]; map{ my $y=$_; $y=~s/U/X/g;$y=~s/O/X/g; print $y;}@x[1..$#x];'
What the script does
It loops over each line, skips first column, and then look for U/O in rest of the data and substitute with X.
It is a quick solution, not well tested. You can check with your file if this works.
Thanks
Priyabrata
• 1 views
•
link
written
by
pbpanigrahi
Now with the example, here a solution with sed and awk as cpad0112 suggested.
$ sed 's/ /\t./g' phylip.txt|awk -v OFS="\t" '{id=$1; gsub("[OU]", "X", $0); $1=id; print}'|sed 's/\t./ /g'
fin swimmer
EDIT:
Pure awk solution:
$ awk '{match($0, "(^[^ ]*)(.*)", s); gsub("[OU]", "X", s[2]); print s[1]s[2]}' phylip.txt
• 1 views
•
link
written
by
finswimmer
Hello cpad,
if I understood this crazy phylip file format correct, it is important how many whitespaces are between the ID and the sequences. There seems to be two possibilities. The ID column is exactly 10 characters width or all sequences in the file start at the same position.
This is why I choose the long way round with sed.
fin swimmer
• 1 views
•
link
written
by
finswimmer
then let us change after 10 characters and skip first line (or whatever number of lines).
May be this should work with one id and one sequence per line:
output: skip first line, then insert a tab after 10 characters, replace characters in all lines except first line and remove the tab. Even if headers contain tabs, it works.
sed '1! s/^.\{10\}/&\t/g' sed.text | awk -v OFS="" 'NR!=1 {gsub("[OU]", "X",$2)}1'
abcdef OUjiOU
abcdefOUjiXX
input:
$ cat sed.text
abcdef OUjiOU
abcdefOUjiOU
Log in to answer this question.
More posts like this
-
Tajima-Nei Distance estimate with BioPerl
written by mauricio.1313 •
-
Trying to show the gaps of each seq on bio::Graphics after converting clustalw
written by jeffrosario92 •
-
How to convert .XFMA file to .nexus file
written by jeccy.J
-
How to convert mugsy (maf) alignment to phylip for RaxML input
written by jeccy.J
-
how XFMA file can be converted to nexus or phylip for further phylogenetic analysis
written by jeccy.J
-
How to write relaxed AND sequential Phylip format using AlignIO?
written by jon.brate
-
How to use and generate consensus sequence from input multiple sequence alignemnt file
written by rkumar •
-
Counting Gaps In A Column Of Alignment
written by sksweety24 •
-
Extracting Matching Subsequence From Pairwise Alignment
written by Woa
-
Taking Out Specific Columns From An Alignment
written by figo
Hello ahmedmagds,
Please use the formatting bar (especially the
codeoption) to present your post better. I've done it for you this time.Thank you!
Hello ahmedmagds,
could you please post an example of your input and desired output? What do you mean by "I tried to use Bioperl but still not successful"? And is perl mandatory or are other solutions fine as well?
fin swimmer
If it's simple replacement, you can try
sed 's/[U,O]/X/g' test.txtBut this will replace any possibility of U or O in sequence header. e.g. a sequence with name AOUxyz. So if you can modify, then it world be great.
I agree with that @pb. However, I am looking for example input from OP. I think most of the new bie's posts do not post data, expected output and error. If i can catch hold of example file with selenocysteine in phylip format, then I can work it out.
Thanks everyone for the help! Sorry for the bad format and not posting a file example. The sed command is easy to use and quick but the problem that sometimes you have U/O in the names which happens one time in the current file but in other examples would be more. so I will put small part of the file here so we can discuss more:
Those can be done away with awk and sed. @OP