Getting the read names/ IDs?

[ghost]

Hello,

I was wondering if there was a way to get the read names for each base. For example, let's say at the position 237 there are 20 reads supporting A and 10 supporting C, is it possible to have the names of the reads supporting the A and the names of the reads supporting the C?

[dridk]

I would like this feature too... @aderzelle Did you find another way to do this ?

[ghost]

@dridk I found bamql https://github.com/BoutrosLaboratory/bamql it does have an option, nt_exact, to specify what nucleotide value you want at position x, relative to the chromosome coordinates. It is slow, though

It could be that https://github.com/vibansal/HapCUT2 also gives that, one would need to parse the compact frag file I guess (I haven't looked into it yet).

[dridk]

Whaoo ! Bamql is awesome ! I added your suggestion on biostar ! https://www.biostars.org/p/322664/#450406

[ghost]

yes, I don't really understand why it does not have more use ^^'

could be very few people need to make such specific queries, I guess.