Hi all,

I have a dataset (see example below) that I would like to go through and copy all lines containing a certain string ("LGIG") plus the line immediately following that line to a new file. I have no problem grepping lines containing the string LGIG but I'm lost how to translate that to line number and shift up one line number for each instance of that string.

Thanks!
Kevin


Example input file:
>14219|CAP|227704
MEFHLPLNRDDLRQRSALNHYVVEEIYPIRIIPSKLHEFNAALRAQGAKCILSHFDVLFSVLLHGANLQPELRSQAWDYLMKVVLKLNTPLVEALEGDTI SSDNCDLLNILKMSVYLLCQFIETFEGEVSKPTVVGATGKGRQKKAQKQAFAELAHLDWDDERERAVRALLQLLQLHLQQLWSPPVVEEDFVNLVTCCCY KMLENQDVVKNKTTRDTIFQVLGVLVKKYNHALGCSLKFIQLLQHFEHLASPLALAVSVFATDYGIKSVAAEIMREIGNMDSKDLARDTSATRAYATFLV ELAEKIPCVMLPSISVLLCLLDGESYSMRNSVLGVLGEMVIRVLSKEELDAKQKCTRDQFLDKLEDHLHDVHAFTRSKTLQIWLAIVNEKALPLPRQHQL LDLVIGRLQDKSSSVRKCALQLLTAVLRMNPFAAKLPLEELKEGYEKESAKLREMQPEEPQLTAADIAKQELEKDWMRMHRGIKKAIALREEKEEDEEEE KGSDEEKVISDEDTVESVKSQIIAALKDGDYLRCVLLFTAAREEWPTDPAFVCSIPMQLDDEEDEGERKNQEVIKCIRNLFFDEVLHSNAFKEALQPQDL DSSRRSSHDTSLVNELTKQQVLVQYLKDSMNFVGQVQQSVPIVCQLLGSKNISDVLESVNFFVTGFKFGVSNSMTGIRRMLVLVWSKEQGVRDAVVEAYR NLYLNLEENNPRLRALAVVNNLTALSLGASLGDLTSLEELVCEFVRSDDLDTHVVQLLWERFAMKIPNTSAAESRAALVLLGMAAGAKVDIVRSNVDVLV KEGLGPRGELDLMLVRDTCSALCKLVPKNKDKTGVAKEPFRFPQDHEAFKRLEFILQNSVSHLESRFWVPMAEQAVNVIYSLAEHPDAITASIIKNVAKQ VVLLERFPLRRETEGEGPLCSPESPTKGVKCPTAYVTRLLSLVGHTALRQLIHLDVSTFGEMKRRHQMQEGKGAKQGTSASARNRSKNNGSSSDDEEDDM GIGGAIAEDAEAEFILKVTEKEVVGGEGLLAALQPLLVGICSNQSKYPDPELQAAASLALAKYMLVSSEFCESHLQLLFTILERSPHAVIRANTIIAMGD LTFRFPNLIEPWTPNLYARLRDPSPQVRKNTLMVLTHLILNDMVKVKGQISELASCIVDDDTRITGLAKLFFHELSRKGNAIYNIMPDMVSRLSDTEVGV DEGNFRVIMKQVECNSEIQKFILCFFCYRYLFSFIQKDRQCESLVEKLCHRFRVTRVERQWRDLAFCLSMLSYSDKSIRKLQENVGCFSDKMADEDVYNS FVTIMNNAKKFAKPETKSLVEEFEQKLEQYHSKGVEDADAEDKAAKSKRSPGRRKAKGRTPGTSRARRQRRSGGDDERDFLGNDSHPTEAGDKPRPKPAI TFDSDDSDIELFKVQEDKAPTQTSQLPSDLENSDPNVSFESPGLRRLKRRHKSVNNKNIVASSSQSSSRTTRARHRNQ
>14219|LGIG|61640
MSFEFTIPINLDCLLSKTNVSQYVVEEVLPLRIIPGAVQDFKFAVRNDNFAVLNHFDTAYSLLSLQKDFEDAIKEEMWDILLQVGQCVTNEIAHGLEDPE LTPDFKLQLLNTLKMTCYLICQFIDMFEVEDTKPGIQINGRGKSKKTTVNKSGRDWEKEKKRGVQTLLNIIQPNLNRLWDPPIAEEEFVNLVSNCCYRLL ENPAIVRTKDIRDVISQLLGVLIKKYNHSLGASLKIMQLLQHFEHLVVPIVQILEIFVNQYQDKSIISELMREIGRLDSRDIAKDTSGTRILAQFLVELS HQLPAAMLPTISVLIGHLDGESYTLRNGVLSVIGEMLVKVLSQENLEDKLKSTRECFLDKLEDHIHDVNAFVRSRVLQIWLHIVNEKCLPLPRQENLVNL ILGRLQDKSSQVRKYAIQLMIALMKNNPFASKLPVEELQASYEKEKEKLKEMTPEETDISDLEECWEAVEKKLRDHVFNHDEDEAETPEEEPSTCTENTQ EELNAEKLCLYKILIMVFREKDLNSEEELEEESERELEIENEGIISCLKDIYLVQKKSEVLSSDPASQPDTSVLNDITKQQVLVQYLKDSTTFADHIQQA VPIICQLLGSKTTSDVFEAIEFFVTGFEFGVTATMLGIRRMLVLIWSSEETIKESVVNAYKRLYLNPNAGGNQRTVALAIVKNLSALLQGASLGDITSFD ALIQQFVKSDDIGQTVIQVLWEKFIPKIPNTTIEESRAALLLLSMIAGAKPEIVKSNIDVLVNEGLGMRAETDMLLAQITCRTLLKLAAGKKTKGEVAAE PFKFSESHEMFQRLSFLLVHCKNTETKVWVPFAEQAINVIYKLSEQPDIIAGEIIKNIAKEVVKSYQSIDIDQPISTVSTLVLTRLLSVSGHVALRHLVH LDSNVFGEMKRRRAIQEEKKEKEQANKKSTRISENIEDELGLAGAAAEDAEAEYIRRICETDIVTGENLLSTLHPLIVAVCTDSTKYPDTQLRTAATLAL AKFTMVSSEFCDAHLQLLFTILEKSPNPAIRANTIIALGDLSFRFPNLIEPWTPHLYGRLRDESAQVRKNTLQVLTHLILNDMVKVKGQISEL
>14237|CAP|174526
MEGPLAVNHHSSSSSSSQEPAPPPAPSLAPAPLAIPQQQVNAASLCRIGQELVHDIVLKASEIFVLLKNMQLPNGATTSQPQNYPETKAKLDDNMKQMLM NFKKLKILYIKVHEHTANLDSRPIEEMLPIEVNGEVKGEGPKLTDEVKYASEEHKEVISQLRNRNRELKQIIDKLRLTVWEINTMLATRKS
>14237|LGIG|86853
PPAGPQQPMVSPNKIVNAATFCRFGQEYIHEIITKATEIFGSRGSFWKCSKFLFQLTNYQDRKTKLEEELKTLCVTFKKLRAIYDKVRETFGGEVETLPV SEVYPYVLEEHRELVETFAHRPRLEPTTTVNMVFFFIFQQVRVKNQQLKEIIDQIRSINWEINTMITMR
>14237|HROB|164180
MSLIVKFENQFHETGVTDVFTEPVDIIPWKFFLLHLFSVQLDSIFDEYLCARINVRDLNLMVIYRSPNSSNDNNDNLIGVLNEFGDLSGRHLVLGYFNFP DIDWKLRVCAGSDYKIENRFLQLIDDRFWLQHVNKPTRYGTNSCPHVLDLIITSEDCVSDLQYSSPLGRKPKKMMSNSQSITNMSQAASSAPPGQTPNQS SSIAINFESPSKKLTPVAMCKYGQELVQEIMHKSYDIFGQLKSIQLVPNERKSRLEDSLHKIEIAFQRLHYVYDSVCNMTKYLEHKPIESYLPMKDKCQD FTTQNIRSDEESKLTEQIQAKNIELKELIDLMRTILWDIDTMMALRKT
>14237|NVEC|200470
MSPLTYTFIAGSSMDSWPSLQSSVMSALKEEEKHLPQPKPPAKPVENPITLTIQGQRTVQEIADKAIELFRKLQNARLSTDTTSRQSQAQTKDIRDQLNT LKGLLAKLRNVYNDTKRAVTIPAGENVENLIPLECSPRNEMDTSDGGGPVEQERAELQKKLKEKNDQLKEVIDKLRIAVWDINTMIMLKPS
>14286|LGIG|234779
MYIASFVLKMVSNRFLVKVAIGGAIFTLTSISGMKIYIENKFQRQDFYLKSMDLLRNYEPAEERLGKPIIDRTINIGELHHNYTDGINARLRIPLKGSIT NGALYVLASRETPKESWHIDKLDLEIPSHFQRWTFYYDPRDSSKVKITGGSDEADETSIIDNNAVQDTQAS

Desired output file for this example:
>14219|LGIG|61640
MSFEFTIPINLDCLLSKTNVSQYVVEEVLPLRIIPGAVQDFKFAVRNDNFAVLNHFDTAYSLLSLQKDFEDAIKEEMWDILLQVGQCVTNEIAHGLEDPE LTPDFKLQLLNTLKMTCYLICQFIDMFEVEDTKPGIQINGRGKSKKTTVNKSGRDWEKEKKRGVQTLLNIIQPNLNRLWDPPIAEEEFVNLVSNCCYRLL ENPAIVRTKDIRDVISQLLGVLIKKYNHSLGASLKIMQLLQHFEHLVVPIVQILEIFVNQYQDKSIISELMREIGRLDSRDIAKDTSGTRILAQFLVELS HQLPAAMLPTISVLIGHLDGESYTLRNGVLSVIGEMLVKVLSQENLEDKLKSTRECFLDKLEDHIHDVNAFVRSRVLQIWLHIVNEKCLPLPRQENLVNL ILGRLQDKSSQVRKYAIQLMIALMKNNPFASKLPVEELQASYEKEKEKLKEMTPEETDISDLEECWEAVEKKLRDHVFNHDEDEAETPEEEPSTCTENTQ EELNAEKLCLYKILIMVFREKDLNSEEELEEESERELEIENEGIISCLKDIYLVQKKSEVLSSDPASQPDTSVLNDITKQQVLVQYLKDSTTFADHIQQA VPIICQLLGSKTTSDVFEAIEFFVTGFEFGVTATMLGIRRMLVLIWSSEETIKESVVNAYKRLYLNPNAGGNQRTVALAIVKNLSALLQGASLGDITSFD ALIQQFVKSDDIGQTVIQVLWEKFIPKIPNTTIEESRAALLLLSMIAGAKPEIVKSNIDVLVNEGLGMRAETDMLLAQITCRTLLKLAAGKKTKGEVAAE PFKFSESHEMFQRLSFLLVHCKNTETKVWVPFAEQAINVIYKLSEQPDIIAGEIIKNIAKEVVKSYQSIDIDQPISTVSTLVLTRLLSVSGHVALRHLVH LDSNVFGEMKRRRAIQEEKKEKEQANKKSTRISENIEDELGLAGAAAEDAEAEYIRRICETDIVTGENLLSTLHPLIVAVCTDSTKYPDTQLRTAATLAL AKFTMVSSEFCDAHLQLLFTILEKSPNPAIRANTIIALGDLSFRFPNLIEPWTPHLYGRLRDESAQVRKNTLQVLTHLILNDMVKVKGQISEL
>14237|LGIG|86853
PPAGPQQPMVSPNKIVNAATFCRFGQEYIHEIITKATEIFGSRGSFWKCSKFLFQLTNYQDRKTKLEEELKTLCVTFKKLRAIYDKVRETFGGEVETLPV SEVYPYVLEEHRELVETFAHRPRLEPTTTVNMVFFFIFQQVRVKNQQLKEIIDQIRSINWEINTMITMR
>14286|LGIG|234779
MYIASFVLKMVSNRFLVKVAIGGAIFTLTSISGMKIYIENKFQRQDFYLKSMDLLRNYEPAEERLGKPIIDRTINIGELHHNYTDGINARLRIPLKGSIT NGALYVLASRETPKESWHIDKLDLEIPSHFQRWTFYYDPRDSSKVKITGGSDEADETSIIDNNAVQDTQAS

Try this:
The -n option quiets all unwanted output. The expression uses /text/ to match any line with "LGIG" in it and +1 to add the following line. Finally 'p' is added to print the matched lines.

Edit: P.S. You should use [code][/code] tags around your code and data to protect the formatting, both of it and the page it's posted on.

grep -A1 pattern file > newfile

I thought about that, but grep puts a "--" group separator between blocks of matches, and I'm not aware of any way to change or turn it off either. I doubt he'd want any of those in his file.

Of course, you could always run it through grep again to filter them out.

Well spotted; I hadn't thought of that, and my sample data only had one
match for the search =}

Your sed solution is cleaner & faster

You guys are great!

Just a note, I had to change the search term to LGIG| so that it didn't also count lines that had amino acid sequences LGIG.

Thanks,
Kevin