multiple reference sequences
Hello , Can I use freebayes on multiple reference sequences . These reference sequences are merged in one FASTA file and each reference start with ">" if yes , there is a special parameter that I need to add to my command line ? for instance , I 'm only using one reference sequence and for that I'm running this command :
freebayes -f path/ref_sequence -L path/BAM_file --vcf path/variants --ploidy 2 --min-alternate-fraction 0.05
• 732 views
•
link
updated
by
GenoMax
•
written
by
emilydolivo97 •
•
written
by
emilydolivo97 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
FreeBayes detection capacity
written by quentinperriere •
-
freebayes calls presumably homozygous variants as heterozygotes
written by johanna.pieplow •
-
freebayes error variant calling
written by evelyn
-
freebayes problem in sample name in vcf
written by evelyn
-
how to get one vcf file for freebayes from many bam files
written by evelyn
-
Variant discrepancy between technical replicates
written by kspata
-
Find variants supported by at least 20 reads freebayes
written by kspata
-
SNPs mining using freebayes: problem in using multiple bam file of different ploidy
written by ram.mainali •
-
How to filter SNPs after Freebayes, if at all, in a fungal haploid genome alignment?
written by Lina F
-
How does FreeBayes evaluates variants that did not meet the --min-alternate-count (-C) and --min-al…
written by Joanne Lim •
It sounds like you need to build and align to a pangenome assembly first, before you can use any of the freebayes functionality. See recent human pangenome paper here that does exactly this.
Appending references together into a single file would probably be very messy, and likely not very informative.