How to identify structural variants manually?
1
2
Entering edit mode
9.7 years ago
mangfu100 ▴ 800

Hi.

I have a question from using pipeline that find the structural rearrangements like SVDetect or BreakDancer.

I ran above software and now I am having a results from SVDetect and BreakDancer.

(Input files is a tumor bam file based on paired end and output files is just a txt file.)

Results have many columns that describe the sequence information like start position, end position, class of

type(indel, translocation,duplication and etc...).

I almost understand why these types are classified by software.

However, I want to see manually that really my tumor bam file has the structural variation which is previously

classified by software.

For example, I used SVDetect which find structural variants and got a output result and its content is below.

chr10     255485     255896     chr10     255738     256047           (F,F,F,F,F,F)     (R,R,R,R,R,R)     TRANSLOC

As I mentioned, I already understood the information above conceptually. but my ultimate goal is just want to identify whether TRANSLOC really happened at above position in tumor bam files or not.

Because many software has mis classification. Therefore, It is required for me to see variants manually not only seeing the output file.

Anyway my question is how you validate the output files ?

If you have your own validate process.. could you explain me in detail?

I heard the IGV program and I already used this tool but I can't find anything information that i want.

I am looking forward to you reply.

Thank you

alignment gene next-gen • 2.3k views
ADD COMMENT
5
Entering edit mode
9.7 years ago

IGV is definitely a good tool for this task. There are some details that might be applicable here

Just FYI, running the structural variant detection on the "normal" sample from the same patient or a pool of "normals" can very greatly reduce your false positive results if you exclude any SV in the tumor that is overlapping one in the normal (or normals).

ADD COMMENT

Login before adding your answer.

Traffic: 2918 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6