This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problems of Pilon Running with Multiple Samples

Hi, all. Recently I tried Pilon wishing to get the fixed vcf and the sequences. But there are many samples, each of which was sequenced in two lanes. I've run Pilon with commands like below

For one sample:
java -Xmx10G -jar pilon-xx.jar --genome ref.fa --bam A_L1.bam --bam A_L2.bam --output pilon.1 --outdir ./ --vcf --fix all,breaks --threads 20 --minmq 30 --minqual 30

For two samples: java -Xmx10G -jar pilon-xx.jar --genome ref.fa --bam A_L1.bam --bam A_L2.bam --bam B_L1.bam --bam B_L2.bam --output pilon.1 --outdir ./ --vcf --fix all,breaks --threads 20 --minmq 30 --minqual 30

But the output was similar that there is only one column for SAMPLE. I wonder is Pilon couldn't verify those group tags in the bam files like samtools? Besides, I got one prefix.fasta file apart from the vcf file. Is that the corresponding sequence of the sample?

I've searched the community and found one closed question. So I hope someone familiar with Pilon could help me with that. Thanks in advance.

snp snp

0 answers

No answers yet.

Log in to answer this question.