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

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 - UML - How would you draw a try catch in a sequence diagram? -

-
Image
is there standard way of drawing try catch sequence diagrams? i've made attempt based on this feel end result doesn't feel right. code diagram based around: public static void save () { try { filehandle filehandle = gdx.files.external(file); filehandle.writestring(boolean.tostring(constantshandler.soundenabled)+"\n", false); (int = 0; < 5; i++) { filehandle.writestring(integer.tostring(constantshandler.highscores[i])+"\n", true); } } catch (throwable e) { } } my attempt note: know still need add loop. that's ok. fragments meant show conditional control flow in sequence diagram. superstructures puts limitations on use of fragments (see pp. 467 of ss2.1.1). should use critical region try part above. , option catch . however, keep telling can take quite freedom in using uml. it's language , such changes in life time. long reader ge...
Read more

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