eliminate hits in a blastn file
0
1
Entering edit mode
8.7 years ago
Daniel James ▴ 10

How to remove a row with an e value less than 95% in a blastn file out put file with a python script

python blast • 1.2k views
ADD COMMENT
1
Entering edit mode

Daniel James, Having seen your previous post, I assume it is hits with an e value of greater than 0.1 that you wish to remove.

ADD REPLY
0
Entering edit mode

Yes e value must be greater than 0.1

ADD REPLY
0
Entering edit mode

what have you tried ?

ADD REPLY
0
Entering edit mode

if(l.split("\t")[2]>=95 & l.split("\t")[10]<=0.1): lin=' '.join((l.split("\t")[0],l.split("\t")[1],l.split("\t")[2],l.split("\t")[3],l.split("\t")[10]))

ADD REPLY
0
Entering edit mode

Use BioPython. Bio.Seq is much better than manual parsing.

ADD REPLY
0
Entering edit mode

if(l.split("\t")[2]>=95 & l.split("\t")[10]<=0.1): lin=' '.join((l.split("\t")[0],l.split("\t")[1],l.split("\t")[2],l.split("\t")[3],l.split("\t")[10]))

ADD REPLY

Login before adding your answer.

Traffic: 2296 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6