bcftools for plant SNPs
1
0
Entering edit mode
5 months ago
Ric ▴ 350

Hi, I ran bcftools in the following way:

bcftools mpileup -O b -o ${3}-${basename}.bcf -f $asm $bam 
bcftools call --ploidy 2 -m -o ${3}-${basename}.vcf ${3}-${basename}.bcf

for a plant genome with --ploidy 2. However, I got this message:

Note: The maximum per-sample depth with -d 250 is 125.0x

PRE-DEFINED PLOIDY FILES

 * Columns are: CHROM,FROM,TO,SEX,PLOIDY
 * Coordinates are 1-based inclusive.
 * A '*' means any value not otherwise defined.

GRCh37
   .. Human Genome reference assembly GRCh37 / hg19

GRCh38
   .. Human Genome reference assembly GRCh38 / hg38

X  
   .. Treat male samples as haploid and female as diploid regardless of the chromosome name

Y  
   .. Treat male samples as haploid and female as no-copy, regardless of the chromosome name

1  
   .. Treat all samples as haploid

Run as --ploidy <alias> (e.g. --ploidy GRCh37).
To see the detailed ploidy definition, append a question mark (e.g. --ploidy GRCh37?).

How can I increase the ploidy level?

Thank you in advance,

SNP sequence • 224 views
ADD COMMENT
1
Entering edit mode
5 months ago

The option --ploidy is not defined to change the ploidy of the reference sequence as you can see in the manual. It's used to define the haploid (sexual) regions of the genome.

--ploidy 2

just like human. So you don't have to change the ploidy

ADD COMMENT
0
Entering edit mode

What should I set for --ploidy value for plants?

ADD REPLY

Login before adding your answer.

Traffic: 2013 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6