How would I find out if a viral gene has been integrated into a sequence?
I am trying to determine if a patient population has integrated HPV into their genome by looking at presence of certain gene sequences. The data is from the TCGA website, which come in BAM files, so I was thinking of doing HTSeq-Count to get the counts, and then mapping the sequences that are counted back to the HPV genome. I'm not sure if this would work, however, and the BAM files from TCGA are both mapped and unmapped reads.
Any guidance? Thanks.
• 905 views
•
link
updated
by
Biostar
•
written
by
idkwhy90 •
•
written
by
idkwhy90 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Changing the reference allele
written by evoecogen •
-
Extracting Unmapped Reads from BAM files using Rsamtools
written by dk0319 •
-
Samtools view Exec format error Help
written by aman-coasab •
-
HTSeq-counts can not map any reads to gtf file
written by nazaninhoseinkhan
-
htseq ignores contiges with zero reads
written by statfa
-
How to extract all reads but concordantly paired from bam file?
written by valerie
-
please help with quantifier.pl error
written by zizigolu
-
HTSeq-count no feature is too high
written by sandeepmarla •
-
How to determine reference genome that was used for alignment given BAM with obscure or no header
written by Irene@Sequencing.com
-
genes counts using htseqcounts
written by bolbolman2000 •
You would have to convert BAM files back to fastq and then remap the reads to HPV alone or HPV+Human genome to be safe. You can then use counting program to get the counts.