Hi all,
I am quite new to RNA-Seq, hoping someone could guide me here.
I am mapping my sequencing reads against the HIV virus + human genome with a reference genome I made by adding the hiv-1 genome to the hg19 reference genome. Mapping was done using STAR and I wanted to view the chimeric junctions using IGV.
However, when I load the chimeric SAM file into IGV, I could only see the mates mapped to the virus genome and I could not visualize any reads mapped to the human genome. Even selecting view position of mate in the hiv view, it will go to the position on the said chromosome but there is no resulting read seen.
Any suggestion on how can I visualize the chimeric reads?
Thank you very much.
•
written
by
heir_of_isildur88 •
1 answer
Using the Sashimi plot option in IGV is a good solution. You need to import a bam/sam file that contains chimeric reads. Using HISAT2 to align the reads did enable me to see chimeric reads.
Not sure about the SA tag but definitely works with the N cigar value.
Log in to answer this question.
-
Detection of chimeric HIV-Host transcripts
written by Mrist01 •
-
Mapping of CDS sequence of genes using minimap2
written by harisanjanas •
-
Different mapping results between HISAT and TopHat
written by concetta •
-
visualize my .BAM file in IGV
written by chrisgr •
-
HIV NL4-3 transcriptome fasta
written by xiaoleiusc
-
Detect Transgene Insertion Site using SAM format
written by kspata
-
Extract chimeric reads
written by Ram
-
Extracting chimeric reads from mapping
written by Anushka •
-
Why reads not match to the genome? too many wrong bases to the genome
written by maple964 •
-
Confirming Viruses detected from NGS data
written by Ron
I think you could use savant It may help in better visualization
It will give you some arch picture that you can easily follow