This is a test version of Biostars. For the public version, visit https://www.biostars.org.
poly ploid consensus seq generation

Hi everyone, I want to generate polyploid consensus sequence for wheat from BAM file for purposes of retrieving the specific sequence of one transcript/ gene.i want to fetch the sequence of transcript/gene abc from rna seq of wheat.what should i do for that?

my command line is: samtools mpileup -uf reference_cdna.fa Batis_Aligned.bam | bcftools call -c --ploidy | vcfutils.pl vcf2fq > Batis_cnpolyploids.fastq but i found this error

call: option '--ploidy' requires an argument

rna-seq snp assembly genome rna-seq

0 answers

No answers yet.

Log in to answer this question.