How to jvarkit/backlocate multiple mutations using a list file?
1
0
Entering edit mode
19 months ago
TraPS-VarI • 0

The following is a command to query a single mutation in FGFR4 -->"FGFR4\tG388R" using jvarkit

echo -e "FGFR4\tG388R" | java -jar dist/backlocate.jar -R hg19.fa --gtf Homo_sapiens.GRCh37.87.gtf

The output is as follows:

FGFR4   Gly 388 Arg FGFR4   ENST00000292408 +   G   1161    GGG AGG|CGG G   chr5    176520242   ENST00000292408.Exon9   .   .
FGFR4   Gly 388 Arg FGFR4   ENST00000292408 +   G   1162    GGG AGG|CGG G   chr5    176520243   ENST00000292408.Exon9   .   .
FGFR4   Gly 388 Arg FGFR4   ENST00000292408 +   G   1163    GGG AGG|CGG G   chr5    176520244   ENST00000292408.Exon9 

. .

I am unable to figure a way to query a long list containing thousands of mutations. Can anyone please suggest a way to query the list (containing ~1000s mutations) (as -i infile) instead of individual mutations? Thanks

java shell echo ubuntu jvarkit • 680 views
ADD COMMENT
0
Entering edit mode

Here is the solution as suggested by the @lindenb

$ cat input_list.txt | java -jar dist/backlocate.jar -R hg19.fa --gtf Homo_sapiens.GRCh37.87.gtf

where input_list contains

ICOS\tF141X CDH10\tL614X KCNE2\tV49X MMP24\tV603X

ADD REPLY
0
Entering edit mode
19 months ago

\t should be tabulations, not a litteral \ followed by t... of course...

that's why it's echo -e "FGFR4\tG388R" and not just echo "FGFR4\tG388R"

ADD COMMENT
0
Entering edit mode

Thanks for the solution. It works now.

ADD REPLY

Login before adding your answer.

Traffic: 1077 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