eliminate hits in a blastn file
How to remove a row with an e value less than 95% in a blastn file out put file with a python script
• 1,689 views
•
link
updated
by
Ram
•
written
by
Daniel James •
•
written
by
Daniel James •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
I need help with a blast command, please help
written by neto.pvp •
-
How read a CSV file skipping the first column using python
written by mandimunari •
-
How to loop through a file to get a value for the coverage in each line of the samtools pileup file…
written by M.O.L.S •
-
Error passing generated fasta to blastn in biopython
written by Daniel
-
eliminating fields in a blastn file
written by Daniel James •
-
eliminating fields in a blastn file
written by Daniel James •
-
batch blast programme
written by Eva_Maria
-
E-values calculation (MEME)
written by anahochmanova •
-
Extract several sequences from genome in FASTA format with genomic coordinates.
written by Ander •
-
Finding Pdb Codes Of Sequences Homologous To An Input Sequences
written by amenabanu •
what have you tried ?
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]))
Use BioPython. Bio.Seq is much better than manual parsing.
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.
Yes e value must be greater than 0.1
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]))