I have a query sequence with 1935 bp and a reference sequence with 2817 bp. My query aligns fully with the reference sequence. I need to extract the regions from the reference sequence(2817bp) which are not a part of my query i.e.approximately 882 bp. I need to do this for many such files (nearly 200). Please help me with a script if possible. I couldn't find any tools online for the same. Kindly help!
•
written
by
pradyumna Jayaram
2 answers
what if i do not know the coordinate regions? I mean as I have mentioned I have to perform this nearly 200 times so it will be difficult to check the coordinates of the reference file which do not match with my query each time. Can you help me with a solution which will directly extract the unaligned region from the reference from a SAM/BAM file?
Log in to answer this question.
-
Mapping to a template sequece
written by QX •
-
Extracting regions around sites in VCF file and writing them as fasta sequences/files
written by Aiswarya •
-
Post mapping extract query sequence mapped to each gene feature
written by L_bioinfo •
-
which tool to use for Transposon structure finding CENSOR or repeat masker ??
written by manaswiniparija3 •
-
How to obtain 200 bp flanking sequence from both sides of BLAST search?
written by anikcropscience •
-
Aligning to ultra-short reference sequence
written by spaladug •
-
Extract alignment blocks of a given region from MAF file having whole genome alignments
written by urjaswita
-
Different results between local and online BLAST (RefSeq database)
written by philipp •
-
Insert a pattern in a sequence
written by jaqy •
-
How to extract reads that map within a defined interval from each other
written by dgimode •
Look at the samtools faidx solution: Extract User Defined Region From An Fasta File
what if i do not know the coordinate regions? I mean as I have mentioned I have to perform this nearly 200 times so it will be difficult to check the coordinates of the reference file which do not match with my query each time. Can you help me with a solution which will directly extract the unaligned region from the reference from a SAM/BAM file?
extract the region with 0 coverage bedtools: extracting no coverage regions
Can you please explain in detail? I am a beginner.
The link that Pierre provided gives the exact command. What additional explanation do you require?