private static final float DELTA = 0.00001f;
// columns 9717- of gene P30419 from Ensembl (much modified)
- private static final String FASTA = ""
- +
- /*
- * forward strand 'gene' and 'transcript' with two exons
- */
+ private static final String FASTA = "" +
+ /*
+ * forward strand 'gene' and 'transcript' with two exons
+ */
">gene1/1-25 chromosome:GRCh38:17:45051610:45051634:1\n"
+ "CAAGCTGGCGGACGAGAGTGTGACA\n"
+ ">transcript1/1-18\n--AGCTGGCG----AGAGTGTGAC-\n"
* reverse strand gene and transcript (reverse complement alleles!)
*/
+ ">gene2/1-25 chromosome:GRCh38:17:45051610:45051634:-1\n"
- + "TGTCACACTCTCGTCCGCCAGCTTG\n"
- + ">transcript2/1-18\n" + "-GTCACACTCT----CGCCAGCT--\n"
+ + "TGTCACACTCTCGTCCGCCAGCTTG\n" + ">transcript2/1-18\n"
+ + "-GTCACACTCT----CGCCAGCT--\n"
/*
* 'gene' on chromosome 5 with two transcripts
+ ">transcript4/1-18\n-----TGG-GGACGAGAGTGTGA-A\n";
private static final String[] VCF = { "##fileformat=VCFv4.2",
- // fields other than AF are ignored when parsing as they have no INFO definition
+ // fields other than AF are ignored when parsing as they have no INFO
+ // definition
"##INFO=<ID=AF,Number=A,Type=Float,Description=\"Allele Frequency, for each ALT allele, in the same order as listed\">",
"##INFO=<ID=AC_Female,Number=A,Type=Integer,Description=\"Allele count in Female genotypes\"",
"##INFO=<ID=AF_AFR,Number=A,Type=Float,Description=\"Allele Frequency among African/African American genotypes\"",
SequenceI gene1 = alignment.findName("gene1");
int[] to = new int[] { 45051610, 45051634 };
int[] from = new int[] { gene1.getStart(), gene1.getEnd() };
- gene1.setGeneLoci("homo_sapiens", "GRCh38", "17", new MapList(from, to,
- 1, 1));
+ gene1.setGeneLoci("homo_sapiens", "GRCh38", "17",
+ new MapList(from, to, 1, 1));
/*
* map 'transcript1' to chromosome via 'gene1'
to = new int[] { 45051612, 45051619, 45051624, 45051633 };
SequenceI transcript1 = alignment.findName("transcript1");
from = new int[] { transcript1.getStart(), transcript1.getEnd() };
- transcript1.setGeneLoci("homo_sapiens", "GRCh38", "17", new MapList(
- from, to,
- 1, 1));
+ transcript1.setGeneLoci("homo_sapiens", "GRCh38", "17",
+ new MapList(from, to, 1, 1));
/*
* map gene2 to chromosome reverse strand
SequenceI gene2 = alignment.findName("gene2");
to = new int[] { 45051634, 45051610 };
from = new int[] { gene2.getStart(), gene2.getEnd() };
- gene2.setGeneLoci("homo_sapiens", "GRCh38", "17", new MapList(from, to,
- 1, 1));
+ gene2.setGeneLoci("homo_sapiens", "GRCh38", "17",
+ new MapList(from, to, 1, 1));
/*
* map 'transcript2' to chromosome via 'gene2'
to = new int[] { 45051633, 45051624, 45051619, 45051612 };
SequenceI transcript2 = alignment.findName("transcript2");
from = new int[] { transcript2.getStart(), transcript2.getEnd() };
- transcript2.setGeneLoci("homo_sapiens", "GRCh38", "17", new MapList(
- from, to,
- 1, 1));
+ transcript2.setGeneLoci("homo_sapiens", "GRCh38", "17",
+ new MapList(from, to, 1, 1));
/*
* add a protein product as a DBRef on transcript1
SequenceI gene3 = alignment.findName("gene3");
to = new int[] { 45051610, 45051634 };
from = new int[] { gene3.getStart(), gene3.getEnd() };
- gene3.setGeneLoci("homo_sapiens", "GRCh38", "5", new MapList(from, to,
- 1, 1));
+ gene3.setGeneLoci("homo_sapiens", "GRCh38", "5",
+ new MapList(from, to, 1, 1));
/*
* map 'transcript3' to chromosome
SequenceI transcript3 = alignment.findName("transcript3");
to = new int[] { 45051612, 45051619, 45051624, 45051633 };
from = new int[] { transcript3.getStart(), transcript3.getEnd() };
- transcript3.setGeneLoci("homo_sapiens", "GRCh38", "5", new MapList(
- from, to,
- 1, 1));
+ transcript3.setGeneLoci("homo_sapiens", "GRCh38", "5",
+ new MapList(from, to, 1, 1));
/*
* map 'transcript4' to chromosome
to = new int[] { 45051615, 45051617, 45051619, 45051632, 45051634,
45051634 };
from = new int[] { transcript4.getStart(), transcript4.getEnd() };
- transcript4.setGeneLoci("homo_sapiens", "GRCh38", "5", new MapList(
- from, to,
- 1, 1));
+ transcript4.setGeneLoci("homo_sapiens", "GRCh38", "5",
+ new MapList(from, to, 1, 1));
/*
* add a protein product as a DBRef on transcript3
@Test(groups = "Functional")
public void testLoadVCFContig() throws IOException
{
- VCFLoader loader = new VCFLoader(
- "test/jalview/io/vcf/testVcf2.vcf");
+ VCFLoader loader = new VCFLoader("test/jalview/io/vcf/testVcf2.vcf");
SequenceI seq = loader.loadVCFContig("contig123");
assertEquals(seq.getLength(), 15);