Hi!
1) I have a bam-file (sorted, with duplicates removed using samtools). It is important that the reads do not completely cover the reference genome.
2) I also have a vcf-file obtained from this bam-file (using bcftools tools) with options (consensus calling)
3) For the reference genome there is a gtf-file with annotation (from ensembl).
I would like to find out some obvious things:
a) What genes have been sequenced in this case (number and name of genes)?
b) I want to analyze the variants for these genes. (How different are these versions from the reference ones, are there any problems - stop codons, deletions, inversions, etc.).
Tell me, please examples of specific commands, what tools is better to achieve this (bedtools? bedops?)
It is worth noting that my data is not on the human genome, so I want to know the solution in general, and not with the help of specific for humans.
That is, how to connect the variant-calling data for these reads and the gene annotation of the reference genome?
Thanks
-
The introduction of single-cell sequencing has advanced the ability to study cell-to-cell heterogeneity. Its use has improved our understanding of somatic mutations1, cell lineages2, cellular diversity and regulation3, and development in multicellular organisms4. Single-cell sequencing encompasses hundreds of techniques with different approaches to studying the genomes, transcriptomes, epigenomes, and other omics of individual cells. The analysis of single-cell sequencing data i
...01-24-2023, 01:19 PM -
Single-cell sequencing is a technique used to investigate the genome, transcriptome, epigenome, and other omics of individual cells using high-throughput sequencing. This technology has provided many scientific breakthroughs and continues to be applied across many fields, including microbiology, oncology, immunology, neurobiology, precision medicine, and stem cell research.
The advancement of single-cell sequencing began in 2009 when Tang et al. investigated the single-cell transcriptomes...01-09-2023, 03:10 PM