X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=test%2Fjalview%2Fext%2Fhtsjdk%2FTabixFeatureReaderTest.java;fp=test%2Fjalview%2Fext%2Fhtsjdk%2FTabixFeatureReaderTest.java;h=0000000000000000000000000000000000000000;hb=f943e5b7a7a5ce2b819495f83dbad28028a9a956;hp=48a932e3a995adf771b84b3091e303a361fe350c;hpb=9cfc8fa0bb0118066ce0406f701b533ca8e84ff6;p=jalview.git diff --git a/test/jalview/ext/htsjdk/TabixFeatureReaderTest.java b/test/jalview/ext/htsjdk/TabixFeatureReaderTest.java deleted file mode 100644 index 48a932e..0000000 --- a/test/jalview/ext/htsjdk/TabixFeatureReaderTest.java +++ /dev/null @@ -1,51 +0,0 @@ -package jalview.ext.htsjdk; - -import htsjdk.tribble.CloseableTribbleIterator; -import htsjdk.tribble.TabixFeatureReader; -import htsjdk.variant.variantcontext.VariantContext; -import htsjdk.variant.vcf.VCFCodec; - -import java.io.IOException; - -import org.testng.annotations.Test; - -public class TabixFeatureReaderTest -{ - // gnomAD exome variant dataset - private static final String VCF_PATH = "/Volumes/gjb/smacgowan/NOBACK/resources/gnomad/gnomad.exomes.r2.0.1.sites.vcf.gz"; - - // "https://storage.cloud.google.com/gnomad-public/release/2.0.1/vcf/exomes/gnomad.exomes.r2.0.1.sites.vcf.gz"; - - /** - * A 'test' that demonstrates querying an indexed VCF file for features in a - * specified interval - * - * @throws IOException - */ - @Test - public void testQuery() throws IOException - { - /* - * if not specified, assumes index file is filename.tbi - */ - TabixFeatureReader reader = new TabixFeatureReader<>( - VCF_PATH, VCF_PATH, new VCFCodec()); - - /* - * gene NMT1 (human) is on chromosome 17 - * GCHR38 (Ensembl): 45051610-45109016 - * GCHR37 (gnoMAD): 43128978-43186384 - * CDS begins at offset 9720, first CDS variant at offset 9724 - */ - CloseableTribbleIterator features = reader.query("17", - 43128978 + 9724, 43128978 + 9734); // first 11 CDS positions - while (features.hasNext()) - { - System.out.println(features.next().toString()); - } - - features.close(); - reader.close(); - } - -}