Collecting specific genome data from a file and collect it in the same title

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
6
down vote

favorite












I have genomes data in a file, genomes-seq.txt. The titles of the sequences begin with >, and then the genome name:



>genome.1
atcg
atcg
atcggtc

>genome.2
atct
tgcgtgctt
attttt

>genome.
sdkf
sdf;ksdf
sdlfkjdslc
edsfsfv

>genome.3
as;ldkhaskjd
asdkljdsl
asdkljasdk;l

>genome.4
ekjfhdhsa
dsfkjskajd
asdknasd


>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.234
efijhusidh
siduhygfhuji

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


I want to collect the similar data for genome.1 in one file so it looks like this:



>genome.1
atcg
atcggtc

iruuwi
sdkljbh
sdfljnsdl
ljhdcj
sdljhsdil
fweusfhygc


But every time I do it using sed I get:



>genome.1
atcg
atcg
atcggtc

>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


That is, multiple genome.1s. How can I do it correctly so on large data set I don't need to remove all the repetitions?






bash







share|improve this question









New contributor




paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.















  • 2




    Hi @paul, what is your sed command that you used?
    – Goro
    10 hours ago











  • I tried but it didn't work
    – paul
    9 hours ago






  • 2




    Show what you tried and we can help fix your errors.
    – glenn jackman
    9 hours ago






  • 1




    Re "...every time I do it using sed...": You ought to include the sed line in the question. Otherwise, this amounts to a work order (using this site as a script-writing service).
    – Peter Mortensen
    7 hours ago











  • The file format is FASTA.
    – Peter Mortensen
    7 hours ago














up vote
6
down vote

favorite












I have genomes data in a file, genomes-seq.txt. The titles of the sequences begin with >, and then the genome name:



>genome.1
atcg
atcg
atcggtc

>genome.2
atct
tgcgtgctt
attttt

>genome.
sdkf
sdf;ksdf
sdlfkjdslc
edsfsfv

>genome.3
as;ldkhaskjd
asdkljdsl
asdkljasdk;l

>genome.4
ekjfhdhsa
dsfkjskajd
asdknasd


>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.234
efijhusidh
siduhygfhuji

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


I want to collect the similar data for genome.1 in one file so it looks like this:



>genome.1
atcg
atcggtc

iruuwi
sdkljbh
sdfljnsdl
ljhdcj
sdljhsdil
fweusfhygc


But every time I do it using sed I get:



>genome.1
atcg
atcg
atcggtc

>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


That is, multiple genome.1s. How can I do it correctly so on large data set I don't need to remove all the repetitions?






bash







share|improve this question









New contributor




paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.















  • 2




    Hi @paul, what is your sed command that you used?
    – Goro
    10 hours ago











  • I tried but it didn't work
    – paul
    9 hours ago






  • 2




    Show what you tried and we can help fix your errors.
    – glenn jackman
    9 hours ago






  • 1




    Re "...every time I do it using sed...": You ought to include the sed line in the question. Otherwise, this amounts to a work order (using this site as a script-writing service).
    – Peter Mortensen
    7 hours ago











  • The file format is FASTA.
    – Peter Mortensen
    7 hours ago












up vote
6
down vote

favorite









up vote
6
down vote

favorite











I have genomes data in a file, genomes-seq.txt. The titles of the sequences begin with >, and then the genome name:



>genome.1
atcg
atcg
atcggtc

>genome.2
atct
tgcgtgctt
attttt

>genome.
sdkf
sdf;ksdf
sdlfkjdslc
edsfsfv

>genome.3
as;ldkhaskjd
asdkljdsl
asdkljasdk;l

>genome.4
ekjfhdhsa
dsfkjskajd
asdknasd


>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.234
efijhusidh
siduhygfhuji

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


I want to collect the similar data for genome.1 in one file so it looks like this:



>genome.1
atcg
atcggtc

iruuwi
sdkljbh
sdfljnsdl
ljhdcj
sdljhsdil
fweusfhygc


But every time I do it using sed I get:



>genome.1
atcg
atcg
atcggtc

>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


That is, multiple genome.1s. How can I do it correctly so on large data set I don't need to remove all the repetitions?






bash







share|improve this question









New contributor




paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.











I have genomes data in a file, genomes-seq.txt. The titles of the sequences begin with >, and then the genome name:



>genome.1
atcg
atcg
atcggtc

>genome.2
atct
tgcgtgctt
attttt

>genome.
sdkf
sdf;ksdf
sdlfkjdslc
edsfsfv

>genome.3
as;ldkhaskjd
asdkljdsl
asdkljasdk;l

>genome.4
ekjfhdhsa
dsfkjskajd
asdknasd


>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.234
efijhusidh
siduhygfhuji

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


I want to collect the similar data for genome.1 in one file so it looks like this:



>genome.1
atcg
atcggtc

iruuwi
sdkljbh
sdfljnsdl
ljhdcj
sdljhsdil
fweusfhygc


But every time I do it using sed I get:



>genome.1
atcg
atcg
atcggtc

>genome.1
iruuwi
sdkljbh
sdfljnsdl

>genome.1
ljhdcj
sdljhsdil
fweusfhygc


That is, multiple genome.1s. How can I do it correctly so on large data set I don't need to remove all the repetitions?






bash




bash






share|improve this question









New contributor




paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.











share|improve this question









New contributor




paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.









share|improve this question




share|improve this question








edited 5 hours ago









Peter Mortensen

82758




82758






New contributor




paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.









asked 10 hours ago









paul

333




333




New contributor




paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.





New contributor





paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.






paul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.







  • 2




    Hi @paul, what is your sed command that you used?
    – Goro
    10 hours ago











  • I tried but it didn't work
    – paul
    9 hours ago






  • 2




    Show what you tried and we can help fix your errors.
    – glenn jackman
    9 hours ago






  • 1




    Re "...every time I do it using sed...": You ought to include the sed line in the question. Otherwise, this amounts to a work order (using this site as a script-writing service).
    – Peter Mortensen
    7 hours ago











  • The file format is FASTA.
    – Peter Mortensen
    7 hours ago












  • 2




    Hi @paul, what is your sed command that you used?
    – Goro
    10 hours ago











  • I tried but it didn't work
    – paul
    9 hours ago






  • 2




    Show what you tried and we can help fix your errors.
    – glenn jackman
    9 hours ago






  • 1




    Re "...every time I do it using sed...": You ought to include the sed line in the question. Otherwise, this amounts to a work order (using this site as a script-writing service).
    – Peter Mortensen
    7 hours ago











  • The file format is FASTA.
    – Peter Mortensen
    7 hours ago







2




2




Hi @paul, what is your sed command that you used?
– Goro
10 hours ago





Hi @paul, what is your sed command that you used?
– Goro
10 hours ago













I tried but it didn't work
– paul
9 hours ago




I tried but it didn't work
– paul
9 hours ago




2




2




Show what you tried and we can help fix your errors.
– glenn jackman
9 hours ago




Show what you tried and we can help fix your errors.
– glenn jackman
9 hours ago




1




1




Re "...every time I do it using sed...": You ought to include the sed line in the question. Otherwise, this amounts to a work order (using this site as a script-writing service).
– Peter Mortensen
7 hours ago





Re "...every time I do it using sed...": You ought to include the sed line in the question. Otherwise, this amounts to a work order (using this site as a script-writing service).
– Peter Mortensen
7 hours ago













The file format is FASTA.
– Peter Mortensen
7 hours ago




The file format is FASTA.
– Peter Mortensen
7 hours ago










3 Answers
3






active

oldest

votes

















up vote
10
down vote



accepted










$sed -nr />genome.1/,/^$/p file | sed '2,$/^>genome.1$/d'

>genome.1
atcg
atcggtc

iruuwi
sdkljbh
sdfljnsdl
ljhdcj
sdljhsdil
fweusfhygc


genome.1 is the key word, change depending on the list you would like to generate.






share|improve this answer






















  • Hi Goro, can I be your sed friend? Where can I find such sed knowledge?
    – schweik
    8 hours ago










  • 1) > matches end-of-word, not the literal character >. 2) Please use . to escape the literal dot in the regular expression that is supposed to match >genome.1. 3) You should also anchor it at the line start and end to avoid false matches: /^>genome.1$/. 4) The -r flag on the first sed command is not required but harmless since you don't use any characters affected by it. 5) As a rule of thumb you should escape sed commands provided through a shell to avoid easily overlooked issues.
    – David Foerster
    6 hours ago


















up vote
6
down vote













With perl



perl -00 -ne 'if (/^>genome.1n/) s/// if $. > 1; print' file





share|improve this answer



























    up vote
    0
    down vote













    With Awk:




     if (/^>/)
     in_section = 0;
     if ($0 == ">genome.1") 
     in_section = 1;
     if (!section_count++)
     print;
     else if (in_section)
     print;



    Usage:



    awk ' if (/^>/) in_section = 0; if ($0 == ">genome.1") in_section = 1; if (!section_count++) print; else if (in_section) print; ' genome.txt





    share|improve this answer




















      Your Answer







      StackExchange.ready(function()
      var channelOptions =
      tags: "".split(" "),
      id: "106"
      ;
      initTagRenderer("".split(" "), "".split(" "), channelOptions);

      StackExchange.using("externalEditor", function()
      // Have to fire editor after snippets, if snippets enabled
      if (StackExchange.settings.snippets.snippetsEnabled)
      StackExchange.using("snippets", function()
      createEditor();
      );

      else
      createEditor();

      );

      function createEditor()
      StackExchange.prepareEditor(
      heartbeatType: 'answer',
      convertImagesToLinks: false,
      noModals: false,
      showLowRepImageUploadWarning: true,
      reputationToPostImages: null,
      bindNavPrevention: true,
      postfix: "",
      onDemand: true,
      discardSelector: ".discard-answer"
      ,immediatelyShowMarkdownHelp:true
      );



      );






      paul is a new contributor. Be nice, and check out our Code of Conduct.









       

      draft saved


      draft discarded


















      StackExchange.ready(
      function ()
      StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2funix.stackexchange.com%2fquestions%2f474268%2fcollecting-specific-genome-data-from-a-file-and-collect-it-in-the-same-title%23new-answer', 'question_page');

      );

      Post as a guest

















      3 Answers
      3






      active

      oldest

      votes








      3 Answers
      3






      active

      oldest

      votes









      active

      oldest

      votes






      active

      oldest

      votes








      up vote
      10
      down vote



      accepted










      $sed -nr />genome.1/,/^$/p file | sed '2,$/^>genome.1$/d'

      >genome.1
      atcg
      atcggtc

      iruuwi
      sdkljbh
      sdfljnsdl
      ljhdcj
      sdljhsdil
      fweusfhygc


      genome.1 is the key word, change depending on the list you would like to generate.






      share|improve this answer






















      • Hi Goro, can I be your sed friend? Where can I find such sed knowledge?
        – schweik
        8 hours ago










      • 1) > matches end-of-word, not the literal character >. 2) Please use . to escape the literal dot in the regular expression that is supposed to match >genome.1. 3) You should also anchor it at the line start and end to avoid false matches: /^>genome.1$/. 4) The -r flag on the first sed command is not required but harmless since you don't use any characters affected by it. 5) As a rule of thumb you should escape sed commands provided through a shell to avoid easily overlooked issues.
        – David Foerster
        6 hours ago















      up vote
      10
      down vote



      accepted










      $sed -nr />genome.1/,/^$/p file | sed '2,$/^>genome.1$/d'

      >genome.1
      atcg
      atcggtc

      iruuwi
      sdkljbh
      sdfljnsdl
      ljhdcj
      sdljhsdil
      fweusfhygc


      genome.1 is the key word, change depending on the list you would like to generate.






      share|improve this answer






















      • Hi Goro, can I be your sed friend? Where can I find such sed knowledge?
        – schweik
        8 hours ago










      • 1) > matches end-of-word, not the literal character >. 2) Please use . to escape the literal dot in the regular expression that is supposed to match >genome.1. 3) You should also anchor it at the line start and end to avoid false matches: /^>genome.1$/. 4) The -r flag on the first sed command is not required but harmless since you don't use any characters affected by it. 5) As a rule of thumb you should escape sed commands provided through a shell to avoid easily overlooked issues.
        – David Foerster
        6 hours ago













      up vote
      10
      down vote



      accepted







      up vote
      10
      down vote



      accepted






      $sed -nr />genome.1/,/^$/p file | sed '2,$/^>genome.1$/d'

      >genome.1
      atcg
      atcggtc

      iruuwi
      sdkljbh
      sdfljnsdl
      ljhdcj
      sdljhsdil
      fweusfhygc


      genome.1 is the key word, change depending on the list you would like to generate.






      share|improve this answer














      $sed -nr />genome.1/,/^$/p file | sed '2,$/^>genome.1$/d'

      >genome.1
      atcg
      atcggtc

      iruuwi
      sdkljbh
      sdfljnsdl
      ljhdcj
      sdljhsdil
      fweusfhygc


      genome.1 is the key word, change depending on the list you would like to generate.







      share|improve this answer














      share|improve this answer



      share|improve this answer








      edited 9 hours ago

























      answered 9 hours ago









      Goro

      8,14153878




      8,14153878











      • Hi Goro, can I be your sed friend? Where can I find such sed knowledge?
        – schweik
        8 hours ago










      • 1) > matches end-of-word, not the literal character >. 2) Please use . to escape the literal dot in the regular expression that is supposed to match >genome.1. 3) You should also anchor it at the line start and end to avoid false matches: /^>genome.1$/. 4) The -r flag on the first sed command is not required but harmless since you don't use any characters affected by it. 5) As a rule of thumb you should escape sed commands provided through a shell to avoid easily overlooked issues.
        – David Foerster
        6 hours ago

















      • Hi Goro, can I be your sed friend? Where can I find such sed knowledge?
        – schweik
        8 hours ago










      • 1) > matches end-of-word, not the literal character >. 2) Please use . to escape the literal dot in the regular expression that is supposed to match >genome.1. 3) You should also anchor it at the line start and end to avoid false matches: /^>genome.1$/. 4) The -r flag on the first sed command is not required but harmless since you don't use any characters affected by it. 5) As a rule of thumb you should escape sed commands provided through a shell to avoid easily overlooked issues.
        – David Foerster
        6 hours ago
















      Hi Goro, can I be your sed friend? Where can I find such sed knowledge?
      – schweik
      8 hours ago




      Hi Goro, can I be your sed friend? Where can I find such sed knowledge?
      – schweik
      8 hours ago












      1) > matches end-of-word, not the literal character >. 2) Please use . to escape the literal dot in the regular expression that is supposed to match >genome.1. 3) You should also anchor it at the line start and end to avoid false matches: /^>genome.1$/. 4) The -r flag on the first sed command is not required but harmless since you don't use any characters affected by it. 5) As a rule of thumb you should escape sed commands provided through a shell to avoid easily overlooked issues.
      – David Foerster
      6 hours ago





      1) > matches end-of-word, not the literal character >. 2) Please use . to escape the literal dot in the regular expression that is supposed to match >genome.1. 3) You should also anchor it at the line start and end to avoid false matches: /^>genome.1$/. 4) The -r flag on the first sed command is not required but harmless since you don't use any characters affected by it. 5) As a rule of thumb you should escape sed commands provided through a shell to avoid easily overlooked issues.
      – David Foerster
      6 hours ago













      up vote
      6
      down vote













      With perl



      perl -00 -ne 'if (/^>genome.1n/) s/// if $. > 1; print' file





      share|improve this answer
























        up vote
        6
        down vote













        With perl



        perl -00 -ne 'if (/^>genome.1n/) s/// if $. > 1; print' file





        share|improve this answer






















          up vote
          6
          down vote










          up vote
          6
          down vote









          With perl



          perl -00 -ne 'if (/^>genome.1n/) s/// if $. > 1; print' file





          share|improve this answer












          With perl



          perl -00 -ne 'if (/^>genome.1n/) s/// if $. > 1; print' file






          share|improve this answer












          share|improve this answer



          share|improve this answer










          answered 9 hours ago









          glenn jackman

          48.6k366105




          48.6k366105




















              up vote
              0
              down vote













              With Awk:




               if (/^>/)
               in_section = 0;
               if ($0 == ">genome.1") 
               in_section = 1;
               if (!section_count++)
               print;
               else if (in_section)
               print;



              Usage:



              awk ' if (/^>/) in_section = 0; if ($0 == ">genome.1") in_section = 1; if (!section_count++) print; else if (in_section) print; ' genome.txt





              share|improve this answer
























                up vote
                0
                down vote













                With Awk:




                 if (/^>/)
                 in_section = 0;
                 if ($0 == ">genome.1") 
                 in_section = 1;
                 if (!section_count++)
                 print;
                 else if (in_section)
                 print;



                Usage:



                awk ' if (/^>/) in_section = 0; if ($0 == ">genome.1") in_section = 1; if (!section_count++) print; else if (in_section) print; ' genome.txt





                share|improve this answer






















                  up vote
                  0
                  down vote










                  up vote
                  0
                  down vote









                  With Awk:




                   if (/^>/)
                   in_section = 0;
                   if ($0 == ">genome.1") 
                   in_section = 1;
                   if (!section_count++)
                   print;
                   else if (in_section)
                   print;



                  Usage:



                  awk ' if (/^>/) in_section = 0; if ($0 == ">genome.1") in_section = 1; if (!section_count++) print; else if (in_section) print; ' genome.txt





                  share|improve this answer












                  With Awk:




                   if (/^>/)
                   in_section = 0;
                   if ($0 == ">genome.1") 
                   in_section = 1;
                   if (!section_count++)
                   print;
                   else if (in_section)
                   print;



                  Usage:



                  awk ' if (/^>/) in_section = 0; if ($0 == ">genome.1") in_section = 1; if (!section_count++) print; else if (in_section) print; ' genome.txt






                  share|improve this answer












                  share|improve this answer



                  share|improve this answer










                  answered 6 hours ago









                  David Foerster

                  948616




                  948616




















                      paul is a new contributor. Be nice, and check out our Code of Conduct.









                       

                      draft saved


                      draft discarded


















                      paul is a new contributor. Be nice, and check out our Code of Conduct.












                      paul is a new contributor. Be nice, and check out our Code of Conduct.











                      paul is a new contributor. Be nice, and check out our Code of Conduct.













                       


                      draft saved


                      draft discarded














                      StackExchange.ready(
                      function ()
                      StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2funix.stackexchange.com%2fquestions%2f474268%2fcollecting-specific-genome-data-from-a-file-and-collect-it-in-the-same-title%23new-answer', 'question_page');

                      );

                      Post as a guest
































































                      -

                      Popular posts from this blog

                      How to check contact read email or not when send email to Individual?

                      Image
                      Clash Royale CLAN TAG #URR8PPP up vote 1 down vote favorite I'm using WordPress 4.9.8, CiviCRM to 5.5.1, I usually send email to contact by Search> Find contacts View contact details Action> Send email Send email ok, Contact received mail ok like picture But status only Email sent though contact read email or not. So, can CiviCRM can change status to Email read when contact read email? wordpress email share | improve this question asked Sep 26 at 0:12 ToanLuong 49 9 add a comment  |  up vote 1 down vote favorite I'm using WordPress 4.9.8, CiviCRM to 5.5.1, I usually send email to contact by Search> Find contacts View contact details Action> Send email Send email ok, Contact received mail ok like picture But status only Email sent though contact read email or not. So, can CiviCRM can change status to Email read when contact read email? wordpress email share | improve this questi...
                      Read more

                      How many registers does an x86_64 CPU actually have?

                      Image
                      Clash Royale CLAN TAG #URR8PPP up vote 2 down vote favorite I am currently learning reverse engineering and am studying the flags register. I had in my mind that rflags was just another name for one of the 16 general purpose registers, for example rax or rbx . But it looks like rflags is actually an additional register. So that makes 17 registers in total... how many more could there be? I have spent at least an hour on this and found numerous different answers. The best answer so far is this, which says that there are 40 registers in total. 16 General Purpose Registers 2 Status Registers 6 Code Segment Registers 16 SSE Registers 8 FPU/MMX Registers But if I add that up, I get 48. Could anybody provide an official answer on how many registers an x86_64 CPU has (e.g. an Intel i7). Additionally, I have seen references to 'hardware' and 'architectural' registers. What are those registers and how many are there? register x86-64 share | improve this...
                      Read more

                      Displaying single band from multi-band raster using QGIS

                      Image
                      Clash Royale CLAN TAG #URR8PPP 1 How can I extract a single band from multi-band raster in QGIS? I have an remote sensed image which has 6 bands (including NDVI band), I want to display each band separately, but have no idea how to do. I have seen some questions similar here but none worked for me. The original image (has 6 bands) is: I want to display the band 6 which should be like this: But I tried gdal_translate, and couldn't get the correct result. What I have got is: qgis raster multi-band share | improve this question edited Mar 5 at 0:53 Summer asked Mar 4 at 6:42 Summer Summer 23 6 Is this any help gis.stackexchange.com/questions/220658/… ? if not gis.stackexchange.com/questions/62133/… might help. – Michael Stimson Mar 4 at 6:46 Thanks for answering but when I used gdal_translate, qgis showed that 'Error 4: Kayena.tif: No such file or directory". Would you know how to fi...
                      Read more