Created
February 6, 2024 21:56
-
-
Save cmdcolin/08e47c53641b5d2152a08933bbd1e546 to your computer and use it in GitHub Desktop.
Convert JBrowse 1 NCList on disk to GFF
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import { LocalFile } from 'generic-filehandle' | |
import NCList from '@gmod/nclist' | |
const filter = new Set([ | |
'Start', | |
'End', | |
'Strand', | |
'Source', | |
'Seq_id', | |
'Name', | |
'Score', | |
'Id', | |
'Type', | |
'Subfeatures', | |
'Parent', | |
]) | |
;(async () => { | |
const store = new NCList({ | |
baseUrl: '', | |
urlTemplate: process.argv[2] + '/{refseq}/trackData.json', | |
readFile: (url: string) => new LocalFile(url).readFile(), | |
}) | |
for await (const feature of store.getFeatures({ | |
refName: process.argv[3], | |
start: 0, | |
end: 50_000_000, | |
})) { | |
const parent = feature.get('parent') | |
const name = feature.get('name') | |
const rest = feature | |
.tags() | |
.filter((f: string) => !filter.has(f)) as string[] | |
if (rest.length !== 0) { | |
console.error(rest.length, rest) | |
} | |
console.log( | |
[ | |
feature.get('seq_id'), | |
feature.get('source'), | |
feature.get('type'), | |
feature.get('start'), | |
feature.get('end'), | |
feature.get('score') || '.', | |
feature.get('strand') === -1 ? '-' : feature === 1 ? '+' : '.', | |
feature.get('phase') || '.', | |
[ | |
`ID=${feature.get('ID')}`, | |
parent ? `Parent=${parent}` : '', | |
name ? `Name=${name}` : '', | |
...rest.map(r => `${r}=${feature.get(r)}`), | |
] | |
.filter(f => !!f) | |
.join(';'), | |
].join('\t'), | |
) | |
} | |
})() |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
usage
run it in the data directories "tracks" folder