This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Detection of chimeric HIV-Host transcripts

I am trying to detect the presences of chimeric HIV-host transcripts in my bulk RNA-seq datasets (single-end reads). The current protocols I have found online seem to align reads to both a human reference and HIV reference at the same time using Bowtie. Then they take reads that aligned to both HIV and human then BLAST those reads against a concatenated Human+HIV genomes.

Does anyone know how I can align to two references at once and pick out the sequences that align to both HIV and Human references? I have used HISAT2 to align my reads to hg38 and HIV reference separately but I do not know how to do this at the same time and then pick out these sequences.

Thank you in advance.

single-end hiv chimeric-transcripts
• 1,033 views • link
updated by Ram Ram • written by Mrist01 • Mrist01

1 answer

Does anyone know how I can align to two references at once and pick out the sequences that align to both HIV and Human references? I have used HISAT2 to align my reads to hg38 and HIV reference separately but I do not know how to do this at the same time and then pick out these sequences.

cat human.fa hiv.fa > merge.fa

index merge.fa

re-map the reads on merge.fa

and then something like: Extracting chimeric reads from mapping

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I did this but then I am un able to index the concatenated genomes. I get hisat2-build command not found. Haven't found a good answer on how to fix this problem.

Is there another way for me to index this file?

• 0 views • link
written by Mrist01 • Mrist01

Log in to answer this question.

More posts like this