ADR, ADF and dp4 counts do not count strands correct
Entering edit mode
2.2 years ago
g.t.los • 0

I work on mutation detection of paired end illumina NGS data, for this I use Varscan and Lofreq, but I have some isseus with strand bias, I hoped I could use the dp4 counts and the ADR and ADF counts but those count the forward and reverse reads instead of the forward and reverse original strands - based on the read pair orientation.

To illustrate an IGV screenshot below, reads are colored by first of pair strand - the forward strand is red, and reverse is blue There are 12 reads with a mutation G > T - 6 forward and 6 reverse - but those are all on the forward strand - So if I want to filter on strandbias the counts of Lofreq and Varscan are useless.

Does somebody know a solution?

IGV screenshot strand bias

sequencing mutation Varscan Lofreq • 524 views
Entering edit mode

You may have to explain the problem a bit more. Be aware that IGV filters your reads, and so do [I presume] VarScan and LoFreq. It is next to impossible to get all of these programs to agree on things.


Login before adding your answer.

Traffic: 1817 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6