This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Seqtk

Hello guys, Apricot has 8 chromosomes but I wanna get just 1 chromosome instead of others using Seqtk tool. Could you pls help me with this issue? Which code Can I use for it?

seqtk
• 2,016 views • link
written by bestone • bestone

What assembly are you using? You should link to it if possible or provide the names of the fasta headers with grep "^>" file.fasta

• 0 views • link
written by rpolicastro rpolicastro

1 answer

use samtools faidx instead, First index the genome

samtools faidx genome.fa

then extract the chromosome of interest

samtools faidx genome.fa chrom1 > subset.fa
• 0 views • link
written by Istvan Albert Istvan Albert

Thank you. for your reply. My file is "Prunus_armeniaca_Stella.fasta.gz" it is not "fa" . I tried it but I couldn't get any result

• 0 views • link
written by bestone • bestone

You need first to extract the FASTA file from the archive.

gzip --decompress --stdout Prunus_armeniaca_Stella.fasta.gz >Prunus_armeniaca_Stella.fasta

samtools faidx Prunus_armeniaca_Stella.fasta title_of_the_chromosome_you_need >one_chromosome.fasta
• 0 views • link
written by shelkmike shelkmike

Hello shelkmike, I tried the codes, it works but I run this code : samtools faidx Prunus_armeniaca_Stella.fasta >one_chromosome.fasta. I got two different files one them is Prunus_armeniaca_Stella.fasta.fai the other one is one_chromosome.fasta. But I have a problem with one_chromosome.fasta. It shows zero byte. Do I have any missing code?

• 0 views • link
written by bestone • bestone

I figure it out guys thank you so much for all of your help

• 0 views • link
written by bestone • bestone

Log in to answer this question.

More posts like this
  • bash script
    written by bestone • bestone
    Hello everyone, there are 32 different apricot genotypes, and I want to analyze their whole genome sequence. I tried to write a bash script code …
  • Error: the required database file annovar/AT_refGeneMrna.txt/hg18_refGene.txt does not exist
    written by bestone • bestone
    Hello guys, I am trying to run annovar using this code: ``` annovar/table_annovar.pl annovar/output.avinput \ annovar/AT_refGeneMrna.fasta \ annovar/AT_refGeneMrna.txt \ -out myanno \ -remove \ -protocol …
  • gff3ToGenePred
    written by bestone • bestone
    Hello guys, I downloaded gff3ToGenePred and tried to convert the gff3 file in Apricot, but gff3ToGenePred gave an error even though it was loaded. this …
  • snpEff.config Error
    written by bestone • bestone
    Hello guys I have a problem with the snpEff. config. I have an apricot gff3 file I converted it to gtf file and then I …
  • conda error
    written by bestone • bestone
    Hello guys, I am using Mac terminal Anaconda3. I couldn't fix the error. I am getting the error as I posted below. Could you pls …
  • CNNScoreVariants Error
    written by bestone • bestone
    Hello guys, I have a problem with CNNScoreVariants code. When I use the code : gatk CNNScoreVariants \ -V input.vcf \ -R reference.fasta \ -O …
  • gatk Hardfilter Error
    written by bestone • bestone
    Hello, I run the code but I got the Error. Could you pls help me with this issue? How can I fix it? gatk VariantFiltration …
  • SnpEff Error
    written by bestone • bestone
    Hello guys, I. run this code : snpEff Prunus_armeniaca_cv_Stella.gff3.gz output.vcf > output.txt I am getting this Error! Could you pls help me with this issue? …
  • gatk4 download issue
    written by bestone • bestone
    Hello guys, I try to download gatk4 to my Mac but I couldn't download it. it gives this error: conda env create -n gatk4 -f …
  • bowtie2
    written by bestone • bestone
    Hello, I wrote this code (test) ➜ dene bowtie2 --very-fast-local -x /Users/uguremre/Desktop/dene/bowtie2_index/bowtie2 -1 /Users/uguremre/Desktop/dene/B11_L03_57_1.fq.gz -2 /Users/uguremre/Desktop/dene/B11_L03_57_2.fq.gz -S /Users/uguremre/Desktop/dene/B11.sam But I got this ERROR. How can …