This is a test version of Biostars. For the public version, visit https://www.biostars.org.
obtain reference and base from read in a region from a bam file

Hi

I want to obtain the reference base and the base of the reads that align in a certain region from a bam file. And get something like this:

      Reference           Read1             Read2
          A                 A           indel or mismatch
          C                 G                  C
          T                 T                  T

Do you know if there is a software to do it? samtools?

alignment

1 answer

sam2tsv : https://github.com/lindenb/jvarkit/wiki/SAM2Tsv

Log in to answer this question.