Biostar Beta. Not for public use.
Pileup file and order of the individuals
1
Entering edit mode
18 months ago
AP • 90

Hello,

I created a pileup file from a list containing names of several bam files:

samtools mpileup -B -f ref.fa -b bam_list.txt > ouput.mpileup

The bam_list.txt looks like that:

  • Ind1.bam
  • Ind2.bam
  • Ind3.bam

Since the pileup format does not contain the individul IDs, it is difficult to know which columns correspond to what individual.

Does samtools give information following exactly the given list? E.g. columns 4 to 6 would correspond to Ind1.bam, 7 to 9, Ind2.bam ...etc (see example of pileup format: https://en.wikipedia.org/wiki/Pileup_format)

It's hard to tell and I need to make sure of what the order is.

Thanks!

ADD COMMENTlink
0
Entering edit mode
18 months ago
AP • 90

After multiple testing, I can confirm that columns in the pileup file are placed in the same order as in the list of bam files provided for the analysis. In other words, columns 4 to 6 would correspond to Ind1.bam, 7 to 9 to Ind2.bam...etc

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1