Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • fahmida
    replied
    Originally posted by Votinhkiem90 View Post
    Thank all, How to sort to max to min lengthcontigs?
    Assume you ran GenoMax's awk code to create the output. In that case the following would sort the contigs based on their length.
    Code:
    sort -t $'\t' -k2 -n -r INPUT_FILE > OUTPUT_FILE

    Leave a comment:


  • Votinhkiem90
    replied
    Originally posted by GenoMax View Post
    I assume you have access to a computer running linux/unix? If you do not then the following will not work.

    In a terminal window run the following command while you are in the directory where you contig file is (in multi-FASTA format). Command reproduced here from the original link in post #3.

    Code:
    $ awk 'BEGIN{flag=0;print "Contig ID\tContig Length\tA\tT\tG\tC\tN\tOtherCharacters"}{if($0~/^>/){if(flag==1){tot=aCount+tCount+gCount+cCount+nCount+xCount;print id"\t"tot"\t"aCount"\t"tCount"\t"gCount"\t"cCount"\t"nCount"\t"xCount;}id=$0;flag=1;aCount=gCount=cCount=tCount=nCount=xCount=0;}else{aCount+=gsub(/[aA]/,"A",$0);tCount+=gsub(/[tT]/,"T",$0);gCount+=gsub(/[gG]/,"G",$0);cCount+=gsub(/[cC]/,"C",$0);nCount+=gsub(/[nN]/,"N",$0);xCount+=gsub(/[^ATGCNatgcn]/,"X",$0);}}END{tot=aCount+tCount+gCount+cCount+nCount+xCount;print id"\t"tot"\t"aCount"\t"tCount"\t"gCount"\t"cCount"\t"nCount"\t"xCount;}' INPUT_FILE >OUTPUT_FILE
    Replace the INPUT_FILE with your contig file name. Replace OUTPUT_FILE with a file name you want to use for the result.
    Thank all, How to sort to max to min lengthcontigs?

    Leave a comment:


  • GenoMax
    replied
    Originally posted by Votinhkiem90 View Post
    Can you more explained for me?
    I assume you have access to a computer running linux/unix? If you do not then the following will not work.

    In a terminal window run the following command while you are in the directory where you contig file is (in multi-FASTA format). Command reproduced here from the original link in post #3.

    Code:
    $ awk 'BEGIN{flag=0;print "Contig ID\tContig Length\tA\tT\tG\tC\tN\tOtherCharacters"}{if($0~/^>/){if(flag==1){tot=aCount+tCount+gCount+cCount+nCount+xCount;print id"\t"tot"\t"aCount"\t"tCount"\t"gCount"\t"cCount"\t"nCount"\t"xCount;}id=$0;flag=1;aCount=gCount=cCount=tCount=nCount=xCount=0;}else{aCount+=gsub(/[aA]/,"A",$0);tCount+=gsub(/[tT]/,"T",$0);gCount+=gsub(/[gG]/,"G",$0);cCount+=gsub(/[cC]/,"C",$0);nCount+=gsub(/[nN]/,"N",$0);xCount+=gsub(/[^ATGCNatgcn]/,"X",$0);}}END{tot=aCount+tCount+gCount+cCount+nCount+xCount;print id"\t"tot"\t"aCount"\t"tCount"\t"gCount"\t"cCount"\t"nCount"\t"xCount;}' INPUT_FILE >OUTPUT_FILE
    Replace the INPUT_FILE with your contig file name. Replace OUTPUT_FILE with a file name you want to use for the result.

    Leave a comment:


  • Votinhkiem90
    replied
    Originally posted by GenoMax View Post
    There is an awk script example on http://www.researchgate.net/post/Ful...ng_contig_file that should produce the info you need. You will have to run it on a unix system and on command line.
    Can you more explained for me? or Is there a software for it
    Last edited by Votinhkiem90; 09-28-2013, 06:46 AM.

    Leave a comment:


  • GenoMax
    replied
    There is an awk script example on http://www.researchgate.net/post/Ful...ng_contig_file that should produce the info you need. You will have to run it on a unix system and on command line.

    Leave a comment:


  • rhinoceros
    replied
    HMMER ships with some excellent little sequence analysis/conversion/etc. tools. I think for this task esl-seqstat would be perfect..

    Leave a comment:


  • Votinhkiem90
    started a topic calcultate Max, min of contigs length

    calcultate Max, min of contigs length

    Hi Everyone
    I want to calculate max, min length of contig file. but i dont know about perl, and i want to mems share script in perl to calculate max, min.
    Thanks.

Latest Articles

Collapse

  • seqadmin
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    06-06-2024, 07:15 AM
  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 02:20 PM
0 responses
9 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-07-2024, 06:58 AM
0 responses
181 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-06-2024, 08:18 AM
0 responses
228 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-06-2024, 08:04 AM
0 responses
184 views
0 likes
Last Post seqadmin  
Working...
X