Re: [Samtools-help] large fragmented reference genome

[Thomas K. <tk...@sa...>]

Hi Kevin - what is the command you are running?

Thomas

On 9/29/14, 7:34 PM, Koh, Chushin wrote:
> Hi all,
>
> I have a very fragmented reference genome (>million scaffolds) that I was able to call variants using Samtools-mpileup V0.19. However, V1.1 hanged on the same data set with no error produced, just never run to completion. The V1.1 still does not work if I extract alignments on a small subset of scaffolds, however, it will work if I also re-header the sub-bam files. Is there anyway to increase the maximum allow scaffolds in the header? Thanks.
>
> With Regards,
> Kevin Koh
>
> ------------------------------------------------------------------------------
> Meet PCI DSS 3.0 Compliance Requirements with EventLog Analyzer
> Achieve PCI DSS 3.0 Compliant Status with Out-of-the-box PCI DSS Reports
> Are you Audit-Ready for PCI DSS 3.0 Compliance? Download White paper
> Comply to PCI DSS 3.0 Requirement 10 and 11.5 with EventLog Analyzer
> http://pubads.g.doubleclick.net/gampad/clk?id=154622311&iu=/4140/ostg.clktrk
> _______________________________________________
> Samtools-help mailing list
> Sam...@li...
> https://lists.sourceforge.net/lists/listinfo/samtools-help



-- 
 The Wellcome Trust Sanger Institute is operated by Genome Research 
 Limited, a charity registered in England with number 1021457 and a 
 company registered in England with number 2742969, whose registered 
 office is 215 Euston Road, London, NW1 2BE.