This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I get this sorted by ascending order? (when the number and character is together)
samtools view NA12878_dedup.bam | cut -f 3 | sort -u > chr.list8

>

chr1
chr10
chr11
chr12
chr13
chr14
chr15
chr16
chr17
chr18
chr19
chr2
chr20
chr21
chr22
chr3
chr4
chr5
chr6
chr7
chr8
chr9
chrM
chrX
chrY

I run the code above and got this list as chr1 chr10 chr11 .... I would like to get this sorted as chr1 chr2 chr3 ....

Can anyone help me out with this?

sort
• 1,590 views • link
updated by Chris Miller Chris Miller • written by woojoy14 • woojoy14

1 answer

If you're using GNU sort, then "sort -V" will do what you want.

In general, though, it's not necesary. The community has (somewhat) standardized on character sort order, which is what you have there.

• 0 views • link
written by Chris Miller Chris Miller

Log in to answer this question.

More posts like this
  • Removing a chromosome from a fasta file
    written by t86dan • t86dan
    Hello, I have a fasta file for human Genome GRCh37 Reference Assembly, but for some reason chromosome 20 is repeated and I would like to …
  • MAFTOOLS HELP - Mutational Signature
    written by David_emir David_emir
    Hi All, I am trying to run mutational signatures on my MAF file, using maftools (https://github.com/PoisonAlien/maftools). while trying to run this i am constantly My …
  • Error with correctGCBias in deepTools
    written by clmin • clmin
    Hello, I'm trying to correct the GC bias with deepTools. I successfully used computeGCBias > computeGCBias -b file.bam --effectiveGenomeSize 2150570000 -g mm9.2bit -l 200 --GCbiasFrequenciesFile …
  • Why the nucleotide sequences at each end of human chromosomes are "NNNNN..."
    written by jinkuozhang • jinkuozhang
    When I was fetching Human Genomics DNA sequence using Coordinates, I found that the nucleotide sequences at each end of human chromosomes(except chr17) are "NNNNN...". …
  • diferent contigs chrM
    written by cristina_sabiers cristina_sabiers
    Hi! I try to call snps to learn how to do it. I got myt bam file, I hadn't made it. I have the snp, …
  • coverageBed with -g option ?
    written by bioguy24 bioguy24
    I am running a coverageBed and getting the below error. ``` cmccabe@DTV-A5211QLM:~/Desktop/NGS$ coverageBed -d -sorted -a /home/cmccabe/Desktop/NGS/bed/bedtools/xgen_targets_sorted.bed -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ERROR: Database file /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam …
  • problem with methylKit of the annotation
    written by hua.peng1314 hua.peng1314
    I am dealing with RRBS data follow [the manual][1]. I met some problem when doing annotation.The problem showed below. ``` > read.transcript.features("mm9.refseq.genes.bed",remove.unsual=TRUE,up.flank=1500,down.flank=500,unique.prom=TRUE) > diffAnn <- …
  • CONTRA problem. It does not stop.
    written by mafonso • mafonso
    Hi, I am trying to run an algorithm for CNV detection and was testing CONTRA. The test file works but using my own files the …
  • How To Retrieve Chromosome Lengths Using Bsgenome Package?
    written by James Ashmore • James Ashmore
    <p>Hello, I am following along with the <a href='http://www.ebi.ac.uk/training/online/sites/ebi.ac.uk.training.online/files/user/42/documents/embo_chipseq_course_1101026.pdf'>EBI's online course for ChIP-Seq analysis</a> and have encountered an error when trying to retrieve chromosome lengths:</p> …
  • Error Import File With "Import" Of Rtracklayer
    written by ska007 • ska007
    <p>Hi, This is related to a the question found here (http://biostar.stackexchange.com/questions/3414?sort=newest#sort-top). I am trying to import a gff file using import("file.gff") which contains a set …