[Bioperl-l] standalone blast outformat

dimitark at bii.a-star.edu.sg dimitark at bii.a-star.edu.sg
Wed Oct 30 03:45:29 UTC 2013


Hi guys,
i was wondering how can i get the query coverage when blasting with  
Bioperl. So i figure out one way is to specify the outformat to  
tabular(7) and the specify that i want the qcovs(query coverage per  
subject). From BLAST manual i see one can do: -outfmt "7 qcovs ..."  
but in Bioperl i can only say: -outformat => 7. If i do: outformat =>  
"7 qcovs" i get an error.

Well another option i suppose is to get the lengths of all hsps from a  
hit and estimate the query coverage that way but im not exactly sure  
if this is correct way.

Any idea how can i get the query coverage?

Thank you very much for your time and any hints

Regards
Dimitar




More information about the Bioperl-l mailing list