Hi All,
I have ULont reads mapped to draft assembly and would like to extract reads that span two contigs at ends/beginings. Read1 from the below example image.
Soft clipping information from the CIGAR can be used to identify such reads but is a tool or any easy way to extract such reads?
any help is much appreciated. thank you
•
written
by
popayekid55
1 answer
As crazy as it sounds, I would recommend blastn for this. Convert the read .fastq into a .fasta and index it with makeblastdb. This will be your subject. Then for each contig split into >contig_left_100bp and >contig_right_100bp prefix and suffix sequences. This will be your query. Run blastn and find subjects with hits above your desired confidence thresholds to two different edges of contigs.
This will also identify circular sequences such as plasmid.
However, if you want to do more than identify such reads (i.e., to perform long-read scaffolding) you can use a tool such as LINKS or longstitch. These may have a way to extract the reads that support scaffolding.
Log in to answer this question.
-
Issue softclipping reads when they belong and don't belong to a common amplicon
written by ManuelDB
-
samtools ampliconclip option --strand is not working?
written by ManuelDB
-
Samtools ampliconclip hard clip not giving same results as soft clip
written by Kevin •
-
Mate reads are aligned to different chromosomes
written by wormball •
-
How to extract targets for specific TF from RegEnrich results
written by Fleur •
-
extracting orphan reads from bam file
written by Assa Yeroslaviz
-
Double soft-clipped reads and how to romeve them
written by Carles Borreda •
-
understanding contig statistics
written by mary •
-
extracting matched part(M) only for the sequence from the CIGAR string
written by Varun Gupta
-
How To Explain Uneven Coverage Of A Dna Seqment Obtained Via Pcr Amplification.
written by rohan
?
Thank you. It was an example image. I have multi-contig file (~19k) which is not feasible to check this way. I want to extract reads that span any 2 contigs, something like reads supporting translocations but in my case I want the reads that spans either beginning or ends
hum... why ?
If I'm understanding correctly, the "spanning" reads need not be between adjacent contigs, so (19000 choose 2) of the above commands would need to be executed. This, of course, could be offloaded to a script.