This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Detecting Structural Variant From Chip Seq Data

Hi,

I am newbie to NGS data analysis.

I used BWA for aligning to reference genome and then samtool. Now am having bam file, sam file and VCF file. i used IGV browser to view the bam file in than i seen IGV uses color coding to flag anomalous insert sizes.

Red for an inferred insert size that is larger than expected (deletion) Blue for an inferred insert size that is smaller than expected (insertion)

i need to extract the deletion part from the bam file or samfile. can anybody give a suggestion of tools on how to extract this particular region ?

ngs chipseq
• 3,617 views • link
updated by Sean Davis Sean Davis • written by sara • sara

1 answer

Some choices to try might include breakdancer (pe), delly (pe and sr), pindel (sr), and cnvnator (rd) where pe = paired-end, sr = split read, and rd = read depth. There are MANY other options and I would be reluctant to suggest which would be "best" for any given situation.

• 98 views • link
written by Sean Davis Sean Davis

Won't it be a problem to detect SV using chip-seq data? I though that chip-seq covers too short regions for SV calling.

• 67 views • link
written by PoGibas PoGibas

I agree that the false-negative rate may be very high, but for all but the read-depth approaches, there is at least the potential for these softwares to find something if it is present in the data. The original poster seems to believe that she can see a signal by eye.

• 69 views • link
written by Sean Davis Sean Davis

@ Sean Davis, through visualization only i seen the signal. Since am new to this field i need help from expertise.

• 0 views • link
written by sara • sara

Log in to answer this question.

More posts like this
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • Can Igv View Blast+ Output File?
    written by shl198 shl198
    <p>Hi, I aligned short reads using blast+ and generated tabular file. I also have <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> files generated by aligners, and I can visualize the …
  • Pindel And Multiple Insert Size Libraries Per Bam File
    written by john.joseph.farrell • john.joseph.farrell
    <p>When running Pindel, the expected insert size is specified for each <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file in the config file.</p> <p>What should be done when there are …
  • Removing Phix Data From Ngs Data For Analysis
    written by crivenster crivenster
    <p>I have a NGS data containing 95-98% of PHIX data,which i dont need. So after some reading, I used <a href='http://bio-bwa.sourceforge.net/'>BWA</a> to align the reads …
  • Disable Soft-Clipping Bwa Aln
    written by Irsan Irsan
    <p>Is there a way to tell <a href='http://bio-bwa.sourceforge.net/'>bwa</a> aln to do not soft-clip reads and try to align them? </p> <p>We see that many soft-clipped …
  • Problem Running Pindel With Bam Files
    written by xenophiliuslovegood • xenophiliuslovegood
    <p>Hi, I am running <strong>Pindel version 0.2.4p, March 26 2012</strong>, on a set of <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files. It chokes on them, printing out:</p> <p><strong>No currentState.Reads …
  • Counting The Whole Insert Size From Paired-End Reads As Coverage
    written by Alastair Kerr Alastair Kerr
    <p>We have updated our workflows for per base sequence coverage to use genomeCoverageBed from <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files. However for pair-end data it seems as though …
  • Preprocessing The Bam For Gatk
    written by Pierre Lindenbaum Pierre Lindenbaum
    <p>The <strong>GATK</strong> often requires that the <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> contains some extra informations in the header ( Group ID, Group Library, SO:coordinate...).</p> <p>Which commands I should …
  • Bwa What Is In .Sai File
    written by Sara • Sara
    <p>Hello, </p> <p>I would like to know what is in .sai file after alignment with <a href='http://bio-bwa.sourceforge.net/'>BWA</a>. what should be remove before convert it to …