This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Removing a chromosome from a fasta file

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 remove it. I know its repeated because when I use grep to look for '>' it shows the following which is the list of all the chromosomes (20 appears 2 times):

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

Any help on removing the repeated chromosome 20 would be appreciated.

Thanks in advance!

fasta chromosome sequence
• 3,700 views • link
written by t86dan • t86dan 
awk 'BEGIN {write=1;} /^>/ { if ($1 in seen) {write=0; next;} else {seen[$1]=1; write=1;}} {if (write==1) print $0;}' your_fasta_file.fa
• 0 views • link
written by cschu181 cschu181

Is there a reason why this happens? dups sequences are understandable, but having one chromosome dups? are you working with cancer data?

• 0 views • link
written by H.Hasani H.Hasani

I would also double check that the sequence is the same.

• 0 views • link
written by st.ph.n st.ph.n

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • How can I get this sorted by ascending order? (when the number and character is together)
    written by woojoy14 • woojoy14
    samtools view NA12878_dedup.bam | cut -f 3 | sort -u > chr.list8 > chr1 chr10 chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr2 …
  • 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 …
  • 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 <- …
  • Splitting chromosomes in bins of 100kb
    written by ChIP ChIP
    Hi! I have the genomic table file of hg19, describing the lengths of each chromosomes. How can I split it up into bins of 100kb …
  • 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 …