This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I create R8 (homopolymer repeat) filter without using illumina pipeline?

Illumina instruments have built-in -or online- analysis software for variant analysis (CASAVA). This software can filter out the false positive variants near the homopolymer repeats (AAAAAAAA) and filter them with "R8" tag.

Is it possible to make homopolymer repeat filter without using Illumina's own pipeline? (with another software?)

P.S I noticed that, GATK have HomopolymerRun script that makes similar job. But it is no longer supported and recommended.

Edit: I posted same question on StackExchange. However, I couldn't find answer yet.

homopolymer-filters repeat
• 4,719 views • link
updated by Ram Ram • written by Can Holyavkin Can Holyavkin

1 answer

I wrote a tool to find the number of homopolymers around a variation: https://github.com/lindenb/jvarkit/wiki/VCFPolyX

$ java  -jar dist/vcfpolyx.jar -R reference.fa input.vcf
(...)
2   1133956 .   A   G   2468.84 .   POLYX=23
2   1133956 .   A   AG  3604.25 .   POLYX=23
2   2981671 .   T   G   47.18   .   POLYX=24
(...)
• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Thank you. It will help much. :)

• 0 views • link
written by Can Holyavkin Can Holyavkin

no sequence dictionary in the reference.

I've created one with picard CreateSequenceDictionary and is in sam format. But still the error persists. Your issue tracker link was broken, that is why I came here. What might be the problem?

• 0 views • link
written by venu venu

Initially I've indexed with samtools but this error appeared saying create with picard, so I've done with picard. But still it is not resolved. Screenshot.

• 0 views • link
updated by Ram Ram • written by venu venu

what the name of the index created by picard. It should be hg19.dict

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thank you. Working well. I've used different name for reference file other hg19 but the VCF file is containing hg19, I thought this was the problem. Anyways thank you once again.

• 0 views • link
written by venu venu

Log in to answer this question.

More posts like this

No similar posts found.