Bioinformatics Stack Exchange is a question and answer site for researchers, developers, students, teachers, and end users interested in bioinformatics. It only takes a minute to sign up.
Sign up to join this community
I'm trying to map ONT long reads to a portion of a gene I'm looking at. The region is about 25bp long. When I search for the region in the document it pulls up the sequence in every read but when I map to reference with minimap2 it produces no contigs.
Is 25bp too small to be a reference? I can potentially extend it but it's right in an area of conservation in this region of a gene and I would rather not extend much farther.
That's not the right way to do it.
You either have reference or not. If not, you need to assemble one (that problem is called genome assembly). If yes, you map all the reads to the whole reference - the process of mapping is made for this - placing individual WGS reads to their corresponding locations and as a matter of fact, missing portions of the reference can mess it up (that was one of the reasons why there was a push for T2T for human).
Once the reads are mapped, you can use any genome browser (e.g. IGV) to zoom in to the exact location you are interested in and check the stack of reads mapping in there. Given how small the region is, manual inspection of the quality of the mapping is the best you got.
