+/*
+ * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
+ * Copyright (C) $$Year-Rel$$ The Jalview Authors
+ *
+ * This file is part of Jalview.
+ *
+ * Jalview is free software: you can redistribute it and/or
+ * modify it under the terms of the GNU General Public License
+ * as published by the Free Software Foundation, either version 3
+ * of the License, or (at your option) any later version.
+ *
+ * Jalview is distributed in the hope that it will be useful, but
+ * WITHOUT ANY WARRANTY; without even the implied warranty
+ * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+ * PURPOSE. See the GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with Jalview. If not, see <http://www.gnu.org/licenses/>.
+ * The Jalview Authors are detailed in the 'AUTHORS' file.
+ */
package jalview.ext.ensembl;
import jalview.analysis.AlignmentUtils;
EnsemblFeatureType[] features = getFeaturesToFetch();
AlignmentI geneFeatures = gffFetcher.getSequenceRecords(accId,
features);
- if (geneFeatures.getHeight() > 0)
+ if (geneFeatures != null && geneFeatures.getHeight() > 0)
{
genomicSequence = geneFeatures.getSequenceAt(0);
}
proteinSeq.createDatasetSequence();
querySeq.createDatasetSequence();
- MapList mapList = AlignmentUtils.mapCdsToProtein(querySeq, proteinSeq);
+ MapList mapList = AlignmentUtils
+ .mapCdsToProtein(querySeq, proteinSeq);
if (mapList != null)
{
// clunky: ensure Uniprot xref if we have one is on mapped sequence
for (DBRefEntry up : uprots)
{
// locate local uniprot ref and map
- List<DBRefEntry> upx = DBRefUtils.searchRefs(upxrefs, up.getAccessionId());
+ List<DBRefEntry> upx = DBRefUtils.searchRefs(upxrefs,
+ up.getAccessionId());
DBRefEntry upxref;
if (upx.size() != 0)
{
else
{
upxref = new DBRefEntry(DBRefSource.UNIPROT,
- getEnsemblDataVersion(), up.getAccessionId());
+ getEnsemblDataVersion(), up.getAccessionId());
}
Mapping newMap = new Mapping(ds, mapList);
// add the new uniprot ref
querySeq.getDatasetSequence().addDBRef(upxref);
}
-
+
}
}
-
+
/*
* copy exon features to protein, compute peptide variants from dna
* variants and add as features on the protein sequence ta-da
*/
- AlignmentUtils.computeProteinFeatures(querySeq, proteinSeq, mapList);
+ AlignmentUtils
+ .computeProteinFeatures(querySeq, proteinSeq, mapList);
}
} catch (Exception e)
{
if (fr.getSeqs().size() > 0)
{
- AlignmentI seqal = new Alignment(
- fr.getSeqsAsArray());
- for (SequenceI sq:seqal.getSequences())
+ AlignmentI seqal = new Alignment(fr.getSeqsAsArray());
+ for (SequenceI sq : seqal.getSequences())
{
if (sq.getDescription() == null)
{
int mappedLength = 0;
int direction = 1; // forward
boolean directionSet = false;
-
+
for (SequenceFeature sf : sfs)
{
/*
// abort - mix of forward and backward
System.err.println("Error: forward and backward strand for "
+ accId);
- return null;
- }
- direction = strand;
- directionSet = true;
-
- /*
- * add to CDS ranges, semi-sorted forwards/backwards
- */
- if (strand < 0)
- {
- regions.add(0, new int[] { sf.getEnd(), sf.getBegin() });
- }
- else
- {
+ return null;
+ }
+ direction = strand;
+ directionSet = true;
+
+ /*
+ * add to CDS ranges, semi-sorted forwards/backwards
+ */
+ if (strand < 0)
+ {
+ regions.add(0, new int[] { sf.getEnd(), sf.getBegin() });
+ }
+ else
+ {
regions.add(new int[] { sf.getBegin(), sf.getEnd() });
}
mappedLength += Math.abs(sf.getEnd() - sf.getBegin() + 1);
}
}
}
-
+
if (regions.isEmpty())
{
System.out.println("Failed to identify target sequence for " + accId
* (havana / ensembl_havana)
*/
Collections.sort(regions, new RangeSorter(direction == 1));
-
+
List<int[]> to = Arrays.asList(new int[] { start,
start + mappedLength - 1 });
-
+
return new MapList(regions, to, 1, 1);
}
int start = sf.getBegin();
int end = sf.getEnd();
int[] mappedRange = mapping.locateInTo(start, end);
-
+
if (mappedRange != null)
{
SequenceFeature copy = new SequenceFeature(sf);
// long start = System.currentTimeMillis();
SequenceFeature[] sfs = sourceSequence.getSequenceFeatures();
- MapList mapping = getGenomicRangesFromFeatures(sourceSequence, accessionId,
- targetSequence.getStart());
+ MapList mapping = getGenomicRangesFromFeatures(sourceSequence,
+ accessionId, targetSequence.getStart());
if (mapping == null)
{
return false;
String type, String parentId)
{
List<SequenceFeature> result = new ArrayList<SequenceFeature>();
-
+
SequenceFeature[] sfs = sequence.getSequenceFeatures();
- if (sfs != null) {
+ if (sfs != null)
+ {
SequenceOntologyI so = SequenceOntologyFactory.getInstance();
- for (SequenceFeature sf :sfs) {
+ for (SequenceFeature sf : sfs)
+ {
if (so.isA(sf.getType(), type))
{
String parent = (String) sf.getValue(PARENT);