Question: Tools for integrating the long reads transcriptomic data into reference genome with gaps in specific region
gravatar for divinomas
5 months ago by
divinomas0 wrote:

We want to investigate the epigenomic state of a specific region (several million base pairs) on human chromosome, which is very polymorphic and notoriously difficult to assemble using short reads. We notice that some long reads transcriptomic data but not long reads genomic data are available for this region. So is there any tool to extract the DNA sequence at all the related polymorphic sites and "patch" the human reference genome with the cell-specific transcriptomic data?

tool assembly genome • 175 views
ADD COMMENTlink modified 5 months ago by ATpoint15k • written 5 months ago by divinomas0

Perhaps you could do SNP/indel calling by aligning the transcriptomic reads with an aligner such as BBMap (maybe using then calling variants with BBTools', then using BCFtools to get the consensus sequence?

ADD REPLYlink modified 5 months ago • written 5 months ago by jean.elbers970
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 683 users visited in the last hour