Yiwei Niu

Hi Lilin, Thank you for your quick reply. I tried this method, but it did not work. In `highlight.text.yadj=rep(1, length(N_ht)), highlight.text.xadj=rep(0, length(N_ht))`, `N_ht` should be the highlighted SNPs? Otherwise it...

Hi Lilin, Thanks a lot for your prompt reply. There were about six SNPs I wanted to highlight. Using `highlight.text.yadj = rep(1, N_ht), highlight.text.xadj = rep(0, N_ht)` could not generate...

Thanks for your reply! I'm sorry for the delayed response. As you recommended, I included several genomes of the same order of my target insect. Then I ran MashMap like...

Thank you for your reply! Sorry, which suggestion do you mean? include an insect genomes in the reference list? I already done that. I include 12 other insects of the...

I counted the number of reads in different categories again. Since I 've used `-s 500 --pi 80`, I used `identity < 80 and length < 1000` as the thereshold....

Hi @biobug16 I used the following code to get log2(CPM/10+1) count from Seurat object ```R expr = CONICSmat::normMat(as.matrix(Seurat::GetAssayData(seurat_obj, assay = "RNA", slot = "counts"))) ``` It was done by extracting...

I am also interested in this issue. Did you find a solution? @dr-ashu-geno

Thank you very much for your reply @dr-ashu-geno. Since xTea cannot merge calls from multiple samples. Recently I tried a new tool to call MEI: [MEGAnE](https://github.com/shohei-kojima/MEGAnE), which showed better performance...

Hi, Serghei. Thank you for your quick reply! Sorry I didn't make myself clear. I've already done the alignment using `STAR` and saved unmapped reads to `FASTQ` files. I don't...

OK! Thank you for your help! Also looking forward to the new version of `ImRep`. BTW, the link in your paper is broken.