Menu

#47 vcf-consensus handling of case

v1.0_(example)
open
nobody
None
1
2014-03-05
2014-03-05
travcollier
No

vcf-consensus throws the "The fasta sequence does not match the REF at..." error if the VCF file and the REF don't agree case wise. This is a bug IMO.

Trivial to fix:
perl/vcf-consensus: line 194

  • if ( $$vline[3] ne substr($$opts{fa_buf},$pos,$ref_len) )
  • if ( uc($$vline[3]) ne uc(substr($$opts{fa_buf},$pos,$ref_len)) )

Discussion


Log in to post a comment.

MongoDB Logo MongoDB