More posts like this
-
Illumina sequencing data trimming help
written by Elena 0Hi! I never posted here before but I wanted somebody's opinion. I did paired end Illumina sequencing on mosquito DNA and ran fastqc to check …
-
How to do preprocessing from fastqc
written by Lulu 0Could someone provide a detailed guide on downstream preprocessing steps based on my FASTQC report? I conducted FASTQC analysis on a paired sample, and here …
-
Distance Matrix
written by Lorenzo 0I have to plot the heatmap representing distances between beta-carbons. (considering alpha-carbons when the CB is missing). I can calculate distances between CA, but when …
-
What do the red bar and black line under the coverage track mean in the IGV screenshot?
written by jonhai 0What do the red bar and black line under the coverage track mean in the IGV screenshot? ![enter image description here][1] [1]: /media/images/a234aaea-bc17-4b13-9401-c17affbf
-
Difference in sashimi plot between single-cell versus bulk RNA-seq data
written by tien 4Hello, I'm using pysashimi to draw splicing junction plot for both single-cell and bulk RNA data but I observe that the coverage rate in bulk …
-
My PCA plot looks good but i don't have differentially expressed genes
written by druggable 6Dear All, I have the following PCA plot from microarray data: ![enter image description here][2] It looks quite good because the green samples and red …
-
About outliers and non -separated samples in PCA
written by Zahra 11Hi all, I have plotted PCA for my samples(Tumor and Normal) in some cancer types. I have used the HTSeq-counts data from TCGA. Then I've …
-
How to plot depth coverage plot?
written by dyfn1947 0Hello I want to plot whole chromosome depth coverage plot each chromosome size is 50Mb and they have 130Kb deletion. but when I made a …
-
DESeq2 R markdown plotting problem
written by Pahi 3Hello! I'm working on my DESeq2 data analysis and I want to make a report with R Markdown. But I have a problem when I …
-
Can not plot degree distribution in Cytoscape Network Analyzer on logscale
written by bhodai 0I am using Cytoscape 3.9.0. When I tried - Tools > Analyze Network, the following window poped-up ![enter image description here][1] When I clicked on …
Are you sure you mean repeat sequences? Do you not mean duplicate sequences?
Yes I'm sure that is repeat sequence
you could hard or soft mask your repeat sequence in the genome, and re-perform mapping.
Or use something like
samtools viewto remove reads overlapping these regions. That is probably better because if a read really comes from these repeats they are properly "decoyed" with full genome alignment. But if you mask the regions they might falsely align somewhere else. Make a BED file with the coordinates you want not included, make the complement against the entire genome (bedtools complement) and then use the-Loption ofsamtools viewto only keep reads that overlap the complement file (which is the genome minus the regions you do not want).i was skeptical of this but here is lh3 saying so also Which Aligners Recognize Soft-Masked Repeats In Reference Sequences? :)
He says
which is I think the best you can do for the aforementioned reasons.
Thanks so much !
Thanks so much !