This is a test version of Biostars. For the public version, visit https://www.biostars.org.
remove the fasta sequence which doesn't contain ATGC

I have a fasta sequence in which some sequence doesn't contain "ATGC" like below so I want to remove all fasta sequences from my file which doesn't contain the proper sequence:-

>hsa-NHS_0029
partial

Thanks in advance

fasta
• 2,215 views • link
updated by cpad0112 cpad0112 • written by harry • harry

partial has two 'a's and one 't' in it. Try

$ seqkit -w 0 grep -srip '[ATGC]' test.fa
• 0 views • link
written by cpad0112 cpad0112

it doesn't remove those fasta sequences with a header in my output file.

• 0 views • link
written by harry • harry

Please post example fasta file, code you have run and output from example fasta file.

• 0 views • link
written by cpad0112 cpad0112

Below is the example fasta file in which you can see a fast sequence is "partial" which is not actually a fast sequence and I want to remove all this sequence from my fasta file with a header.

>hsa-NHS_0029
partial

>hsa-MAP3K2_0035
CAAAGCAGACACAGTCTTTTTCCTTATGGGGCTTATAGTGTGGTATTGGTGGTCTTCAAGCTTTAGGTGCATGCTCTCTAGGGGCAGCACAGGCCAGTCTTTGGGAAGTAGAAAGAAAACACTAGAATTTTAATCTTGTAAAAATATTTATTTTGGAAGTATGCCTACATGATATATGTACCCAAAAATGTTGAGTTGATAAGGATACAGAAATCAAAATTTGAATCCCTGGAGAACAACTTTATTCTTTTTAACAGTTTCATTTCATTCTGTCATACAGATGTCACATAGTTAAGTTTAATCAGTTTCCTGAATAAACATTTTGCATTGCTCTTTAACATTCATATGGTAGTTTCAGTGCAAATTAAAAACTGCAGAAAATTTTGAGAATTTCCTTTTTAAAATTTGTATTAGGATAAAACAGACTTCAGAGTTTTTGAACTTTGTGTTGAAAATGTGTGCCTGTGTTAAAAAATGTGCTGTTGTAGTGTTAAGTAAAGTGAATATTTACGTAATTTTTAGTTAATCTTGTATCACTCTAAAAGGAG

>hsa-TMOD2_0026
ATGTTGTCAAAGGTGAAAAAGTAAAGCCAGTATTTGAGGAACCACCAAATCCCACAAATGTGGAAATAAGCCTGCAGCAGATGAAAGCCAATGATCCTAGCTTGCAAGAAGTCAACCTCAACAACATTAAGTTCGTAAGAAGAGAGTTGAAGCAGACCGAAG

As you suggested a seqkit command, I ran it but it don't remove those sequences from my file.

seqkit -w 0 grep -srip '[ATGC]' test.fa >test_new.fa

So I want the output is like below :- In which the "partial" fasta sequence is removed with a header.

>hsa-MAP3K2_0035
    CAAAGCAGACACAGTCTTTTTCCTTATGGGGCTTATAGTGTGGTATTGGTGGTCTTCAAGCTTTAGGTGCATGCTCTCTAGGGGCAGCACAGGCCAGTCTTTGGGAAGTAGAAAGAAAACACTAGAATTTTAATCTTGTAAAAATATTTATTTTGGAAGTATGCCTACATGATATATGTACCCAAAAATGTTGAGTTGATAAGGATACAGAAATCAAAATTTGAATCCCTGGAGAACAACTTTATTCTTTTTAACAGTTTCATTTCATTCTGTCATACAGATGTCACATAGTTAAGTTTAATCAGTTTCCTGAATAAACATTTTGCATTGCTCTTTAACATTCATATGGTAGTTTCAGTGCAAATTAAAAACTGCAGAAAATTTTGAGAATTTCCTTTTTAAAATTTGTATTAGGATAAAACAGACTTCAGAGTTTTTGAACTTTGTGTTGAAAATGTGTGCCTGTGTTAAAAAATGTGCTGTTGTAGTGTTAAGTAAAGTGAATATTTACGTAATTTTTAGTTAATCTTGTATCACTCTAAAAGGAG

    >hsa-TMOD2_0026
    ATGTTGTCAAAGGTGAAAAAGTAAAGCCAGTATTTGAGGAACCACCAAATCCCACAAATGTGGAAATAAGCCTGCAGCAGATGAAAGCCAATGATCCTAGCTTGCAAGAAGTCAACCTCAACAACATTAAGTTCGTAAGAAGAGAGTTGAAGCAGACCGAAG
• 0 views • link
written by harry • harry

do you really have string 'partial' as your sequence? didn't understand that. Here is the solution:

$ seqkit -w 0 grep -svp "partial" test.fa

>hsa-MAP3K2_0035
CAAAGCAGACACAGTCTTTTTCCTTATGGGGCTTATAGTGTGGTATTGGTGGTCTTCAAGCTTTAGGTGCATGCTCTCTAGGGGCAGCACAGGCCAGTCTTTGGGAAGTAGAAAGAAAACACTAGAATTTTAATCTTGTAAAAATATTTATTTTGGAAGTATGCCTACATGATATATGTACCCAAAAATGTTGAGTTGATAAGGATACAGAAATCAAAATTTGAATCCCTGGAGAACAACTTTATTCTTTTTAACAGTTTCATTTCATTCTGTCATACAGATGTCACATAGTTAAGTTTAATCAGTTTCCTGAATAAACATTTTGCATTGCTCTTTAACATTCATATGGTAGTTTCAGTGCAAATTAAAAACTGCAGAAAATTTTGAGAATTTCCTTTTTAAAATTTGTATTAGGATAAAACAGACTTCAGAGTTTTTGAACTTTGTGTTGAAAATGTGTGCCTGTGTTAAAAAATGTGCTGTTGTAGTGTTAAGTAAAGTGAATATTTACGTAATTTTTAGTTAATCTTGTATCACTCTAAAAGGAG
>hsa-TMOD2_0026
ATGTTGTCAAAGGTGAAAAAGTAAAGCCAGTATTTGAGGAACCACCAAATCCCACAAATGTGGAAATAAGCCTGCAGCAGATGAAAGCCAATGATCCTAGCTTGCAAGAAGTCAACCTCAACAACATTAAGTTCGTAAGAAGAGAGTTGAAGCAGACCGAAG
• 0 views • link
written by cpad0112 cpad0112

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • extract the all columns from 1 file giving a list of another file
    written by harry • harry
    I have 1 text file like this below : X 119614423 119614883 intergenic 52 NA hsa-intergenic_023278 X 123867152 123873856 XIAP 25 NA hsa-XIAP_0027 X 129496750 …
  • Fasta header contain multiple string
    written by harry • harry
    I have large fasta file. As you see below there are > sign present in some fasta header like >exon2_ENST00000218032|>exon2_ENST00000218032 >exon17_ENST00000253024|>exon17_ENST00000253024 I want to remove …
  • Remove fasta sequence on the basis of header name
    written by harry • harry
    >exon18_ENST00000194900|exon19_ENST00000194900 GCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGCCATTTTCATCAAGCCCAAGTCCATTGAAGCCCTTATTACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATG >exon18_ENST00000194900|exon21_ENST00000194900 TCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAATACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATGG >exon18_ENST00000194900|exon20_ENST00000194900 CAGACATATGAACAAGCAAATAAGATCTATGACAAAGCCATGAAACTGGAGCAGGAATTTGGAGAGTACTTTACAGTACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATGG >exon18_ENST00000194900|exon21_ENST00000194900 CTCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAATACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATG >exon17_ENST00000194900|>exon17_ENST00000194900 ATGATGACCTGATCTCCGAATTTCCACATAAATTTGGATCCTGTGTGCCACTCACTATGCAAGGCCTGTGATCATCCTGGGCCCAATGAAGGACCGAGTCA >exon19_ENST00000194900|exon21_ENST00000194900 TCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAAGCAAGCACTGCATCTTAGATGTTTCCGGCAATGCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGCC >exon19_ENST00000194900|exon20_ENST00000194900 CAGACATATGAACAAGCAAATAAGATCTATGACAAAGCCATGAAACTGGAGCAGGAATTTGGAGAGTACTTTACAGGCAAGCACTGCATCTTAGATGTTTCCGGCAATGCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGCC >exon19_ENST00000194900|exon21_ENST00000194900 CTCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAAGCAAGCACTGCATCTTAGATGTTTCCGGCAATGCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGC In this example of fasta sequence, you see there …
  • Want to reverse the fasta sequence
    written by harry • harry
    >ENST00000014935 TAGGTTTTTAATTACCTTTATTTTATTTTGCCAAACATACCTGGGAATACCTTTTATTTTTTTTTTACCTTGGGGTGATGGTTCCAAACCATAAATGTGATTATAGTTAACACATGACCCTTCTAGCGTCCCAGCCAGTGTTTTTCCTGACCTCTCTTCTTTGGAGAGGAGGATGGAAGGGAGGGGTCCGGCATGCTGCTGGCATTTTGCTGTGTCCTGCAGCCCCTTTCCGGGACACCTGGGTTCACACAGCTTTTTAGCTTACATAACTGGTGCAGATTTTCTGTGTGGAGATGTTGCCTTGACCAGCCTTGGCTGGACTTTACCAGGCATGCAGAAGCCTGTACCAACACAGACTACAGCACCCAGGAGGTGCGAGTGTGGCTGCTCAGCGGTTATAACAGGCCTGACTGCATTGTTCACCGGATTATAATGAGCCAAAATGTTTCCCGGTGTTTGCTGGTTTCAGGGAAGGAGTTTGATATAGCAGATTAACCACCCTCCTTGTAGCTATTGGGGCTTAATGGTTTCCTGGTGATTCTTACCAATCCACAATAAACATGGCCCATTGGCATATCTGCTGCACAAGTGTCCTATCTCACCAATCTGGGTTTTTGTTCTCAGTAACTTTCCTTCTTGTCATACAACATCTTCATTCCTCTTTCTGAACCCTCCCTTCCCCTACCCCAACCCAGAGCCCACTTTGTCTCCACTCCTGATACTACACTACCTGGCAGGTGGCATGAGTGCAGGGCCCCTGGCTTCCTCTCCTAATCTAGGCACAAGCCCAACCAAAGAACAAGAGCCAAATCAAACAAGGCAGGCAGGGGTGGACTACAGTCACAGGGCAACTATAGTTGAAGCCCCCCAGCCCCAGGGCTGGATGGACGGGGGAGGCTGGGGTTTAAGTCCCAAAAGGCAGCAGGCCCTGGGGGGGTAGGGGGACGCTCAGGCAGCAGGGCACAGCTGAGGGGACAGGAGTGATAGCAGCAACAGAACAGTGAGGCTGAGAGGCTGGACGCTGTGCGCCTGGCTCAGCTTCAGCTCCACCTCCACGGGGTAGTGGTCACTGATGTTGAGGGC >ENST00000014944 TCCTCCTCGGTGAGCTGGAAGCTCGTGGGGAAGTCAAAGGCAGCCGCAGTGTGCAGCAGACTCCGGCAGCGCTCCCCGTGCAGCACGACGCGGTCATAGGTGCAGTGGGTGCTGGCCCGCACTGTGGTGTCCTCCCCATCGGCAATCACCCAGTGGAAGCCTGGCTCAGTCCGCAGCTCCAGCTTGTCCAGGCGCTTTTTGGTCAGTGAAGCGCAGTCAGCATTGAAGTCCCCAAGCAGGATCACGTC Likewise, I have a huge fasta sequence and I want to reverse the sequence of those fasta sequences like below:- >ENST00000014935 …
  • creating virtual library preparation of all possible exons in a gene in human
    written by harry • harry
    Hi everyone, I want to make a virtual library fasta file for all the possible exons in a gene and also all possible exon-intron fasta …
  • trimming of fasta file
    written by harry • harry
    As I have multiple fasta sequences like see below: >read1 CAGGTCTGGCTGGATGAAGGGCACGGCATAGGTCTGACCTGCCAGGGAGTGCTGCATCCTCACAGGAGTCATGGTGCTGCTGAAGATGTCTCCAGAGACCTTCTGCAGGTACTGCAGGGCATCCGCCATCTGCTGGACGGCCTCCTCTCGCCG >read2 CATCTGCGGAGGCTGCCGTGACGTAGGGTATGGGCCTAAATAGGCCATTGTGAGTCATGAGCTTGGTCTGTAGAGGCTGACTGGAGAAAGTTCTGGGCCTGGAGAGGCTGCCGGGAGGTAGGAGTGGTGAGGTCGACTTGAGAAAGTTCAGGGCCTGGAGAGAAGGCTGGGAGGCAGGAGCTGGGTCTAAAGAGGCCATTGTAACGATGGAGCTGTGCCTGTGGAGGCTGTTGTGAGGCAGTAGCCT >read3 TTGAGGTGGGAGGATCGCTTCAGCCTGGAAGGTTGAGGCTGCAGTCAGCTGCGATAGCACTACTACACTCCAGCCTTGGACAACAGAGGGAGACCTTTCGCTGTCACCCCTCTAGAATCCACGTATACGAAAATTCCAAATGTTAGTTGGGCATAGTGGCAAGCACCTGTAGTCTCAGCCACGTGGGAGG These are in different lengths so I want to isolate …
  • Processing corresponding fasta records from separate multi-fasta files
    written by Digsby • Digsby
    I am trying to align corresponding fasta records from separate multi-fasta files which all contain the same number of fasta records. Each multi-fasta file contains …
  • Extract Sequences from Multiple sequence alignment containing special character like gap?
    written by MSRS MSRS
    I have multiple sequence alignment files of DNA in FASTA which contain deletions and insertions. I want to extract sequences that only contain a gap …
  • Parsing dataframe and sort a fasta file
    written by Chvatil Chvatil
    Hi all, I need some help. In fact in the contexte of my work I have a dataframe of seq names paired together, in each …
  • How to remove some fasta sequences by header information from a large fasta file, any command and s…
    written by seta seta
    Dear all, I would like to remove some fasta sequences within a large fasta file in term of header information (sequence name), could anybody please …