multiple reference sequences
0
0
Entering edit mode
18 days ago

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
freebayes • 164 views
ADD COMMENT
2
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

Traffic: 1149 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