I have a BWA alignment of reads to a small (50 kb) reference sequence. It is a transgenic sequence inserted into a host cell genome. I want to be able to locate the insert's position in the host cell genome. There are reads at the ends (pointing outwards) which have their pairs unmapped. These mates would presumably be in the flanking genomic sequence that I want to identify. Is there an easy way to get the unmapped mates? I suppose I could make a list of the reads and write a script to parse the original fastQ files, but I am hoping there is a tool already available for this (seemingly common) purpose. Any help would be greatly appreciated.
-
-
See this thread for the info you need: http://seqanswers.com/forums/showthread.php?t=12283 -
You could parse the .bam for unmapped reads whose mates mapped close to your boundaries in the correct orientation.
You could also align the fastqs to the 50kb genome and the host genome, then filter for reads that aligned to the host whose mates aligned to the insert. That's probably the best solution. You'd want the mapping position of the reads that aligned to host anyway, so this way you'd have them.Comment
-
Thanks for the suggestions. I'll look at the other thread and try your suggestions.Comment
-
Covid-19
I am confused on the issue:
the service provider company provide AmpliSeq for Illumina On-Demand, Custom, and Community Panels. for COVID diagnostic, the library was prepared, the issue started @ sample sheet, manifest file - covid- successfully added, genome (we try our level best to integrate the genome file but no use, after creating multipath the genome was integrated in sample sheet and run started, output was 93.8=Q score, Cluster passing 96.7%, Cluster density 774K) but analysis failed (Sunday) till now we try all possible methods with illumina support but no use. initially RNA amplicon was downloaded and added in sample sheet, the sample sheet was headed by DNA amplicon (no use) & now PCR amplicon was added in sample sheet but same error. plz guide.
Comment
-
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...10-14-2026, 09:55 AM -
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment