grep - Retrieve lines of a text that contain a certain string only in the end and not somewhere in between -

-

i have text file of taxonomic assignation of bacterias looks ( numbers indicate different bacterias): 

1130952 k__bacteria; p__acidobacteria; c__acidobacteriia; o__acidobacteriales; f__acidobacteriaceae; g__edaphobacter; s__modestum 555445  k__bacteria; p__firmicutes; c__clostridia; o__clostridiales; f__lachnospiraceae; g__clostridium; s__fimetarium 325910  k__bacteria; p__firmicutes; c__clostridia; o__clostridiales; f__ruminococcaceae; g__; s__ 744205  k__bacteria; p__proteobacteria; c__deltaproteobacteria; o__; f__; g__; s__

many of bacterias don´t have classification down specie level, lack information : "s__". see bacterias have information (as in 2 bacterias above, 1 being "s__modestum" , other "s__fimetarium"). using mac terminal (mac os x 10.9.5) , tried, 

grep -v "s__" file

but since assignation contain s__ noting (it excludes them guess..).

i have tried using * @ end in s__* doesn't work either.

what apply command , line , count of bacterias species assignation.

1 1130952 k__bacteria; p__acidobacteria; c__acidobacteriia; o__acidobacteriales; f__acidobacteriaceae; g__edaphobacter; s__modestum 1 555445    k__bacteria; p__firmicutes; c__clostridia; o__clostridiales; f__lachnospiraceae; g__clostridium; s__fimetarium

just ask grep match character after s__ not being end of line (expressed $):

$ grep 's__[^$]' file 1130952 k__bacteria; p__acidobacteria; c__acidobacteriia; o__acidobacteriales; f__acidobacteriaceae; g__edaphobacter; s__modestum 555445  k__bacteria; p__firmicutes; c__clostridia; o__clostridiales; f__lachnospiraceae; g__clostridium; s__fimetarium

to count of lines matching condition, need use awk store counter values array:

$ awk '/s__[^$]/ {a[$0]++} end {for (i in a) print a[i], i}' file 1 555445  k__bacteria; p__firmicutes; c__clostridia; o__clostridiales; f__lachnospiraceae; g__clostridium; s__fimetarium 1 1130952 k__bacteria; p__acidobacteria; c__acidobacteriia; o__acidobacteriales; f__acidobacteriaceae; g__edaphobacter; s__modestum

to make check appearing @ end of file, need checks:

grep -e 's__[^ $]+$' file awk '/s__[^ $]+$/ {a[$0]++} end {for (i in a) print a[i], i}' file

they check after s__ there set of @ least 1 character not being space or end of line. , then, end of line.

update

thank worked great! there way can sum lines know how many counts have in total of non "s__"? – isa

sure, add print length(a) see how many elements array has:

$ awk '/s__[^ $]+$/ {a[$0]++} end {for (i in a) print a[i], i; print length(a)}' 1 555445  k__bacteria; p__firmicutes; c__clostridia; o__clostridiales; f__lachnospiraceae; g__clostridium; s__fimetarium 1 1130952 k__bacteria; p__acidobacteria; c__acidobacteriia; o__acidobacteriales; f__acidobacteriaceae; g__edaphobacter; s__modestum 2

Comments

Post a Comment

Popular posts from this blog

c++ - No viable overloaded operator for references a map -

-
i trying use map can make tag names reference number. when try , use it, in code error (6 in total every time reference map): src/main.cpp:25:45: error: no viable overloaded operator[] type 'std::map<std::string, std::string>' const char* idcs = node.child_value(tagmap[3]); this code: #include "pugi/pugixml.hpp" #include <iostream> #include <string> #include <map> int main() { pugi::xml_document doca, docb; std::map<std::string, pugi::xml_node> mapa, mapb; std::map<std::string, std::string> tagmap {{"1", "data"}, {"2", "entry"}, {"3", "id"}, {"4", "content"}}; if (!doca.load_file("a.xml") || !docb.load_file("b.xml")) { std::cout << "can't find input files"; return 1; } (auto& node: doca.child(tagmap[1]).children(tagmap[2])) { co
Read more

java - Custom OutputStreamAppender not run: LOGBACK: No context given for <MYAPPENDER> -

-
i wish write custom appender on basis of outputstreamappender . wrote following class package tests; import java.io.ioexception; import java.io.outputstream; import ch.qos.logback.core.outputstreamappender; public class myappender<e> extends outputstreamappender<e> { public myappender() { system.out.println("myappender created"); setoutputstream(new outputstream() { @override public void write(int b) throws ioexception { throw new unsupportedoperationexception(); } }); } } then wrote following runner: package tests; import org.slf4j.loggerfactory; import ch.qos.logback.classic.logger; public class runner { private static final logger log = (logger) loggerfactory.getlogger(runner.class); public static void main(string[] args) { //log.addappender(new myappender<iloggingevent>()); log.info("hello world"); } } finall
Read more

java - Cannot secure connection using TLS -

-
i'm trying establish secure socket connection between java client applet (built jdk 1.7.0_75-b13) , vc++ server application. as test vehicle, used vc++ client/server sample found in msdn forums, modified use schannel , able establish socket using cipher suite tls_rsa_with_aes_128_cbc_sha. works of tls 1.0/1.1/1.2. when try opening socket java applet same server application, connection rejected server reporting following: tls 1.0 acceptsecuritycontext failed: 0x80090327 tls 1.1 acceptsecuritycontext failed: 0x80090331 tls 1.2 acceptsecuritycontext failed: 0x80090331 this java code used create socket: debugprint("setting secure connection"); sslsocketfactory sslsocketfactory = (sslsocketfactory) sslsocketfactory.getdefault(); sslsocket sslsocket = (sslsocket) sslsocketfactory.createsocket("127.0.0.1", socketnumber); debugprint("starting handshake"); sslsocket.settcpnodelay(true); sslsocket.setsolinger(false, 0); sslsocket.s
Read more