From: Tom B. <tb...@um...> - 2011-01-29 22:15:58
|
Before indexing, the records in a .bam file must be sorted by position ("coordinate"). Use 'samtools sort'. That's probably the issue here. This email just now came through to the list. - tom blackwell - On Fri, 28 Jan 2011, Christos Kozanitis wrote: > Hi > > I am trying to index a 6GB bam file of Illumina data, which was originally in export format and was converted to sam by export2sam.pl. > > However the bai file that I obtain has a size of 40 bytes and when I run: "view sam_sorted.bam chr3" I get only 27 reads in the output (out of the 21M reads that I expect). > > Is there any step that I have skipped in the procedure? > > My version of samtools is 0.1.6. > > > Thanks > Christos > > ------------------------------------------------------------------------------ > Special Offer-- Download ArcSight Logger for FREE (a $49 USD value)! > Finally, a world-class log management solution at an even better price-free! > Download using promo code Free_Logger_4_Dev2Dev. Offer expires > February 28th, so secure your free ArcSight Logger TODAY! > http://p.sf.net/sfu/arcsight-sfd2d > _______________________________________________ > Samtools-help mailing list > Sam...@li... > https://lists.sourceforge.net/lists/listinfo/samtools-help > > > |