There are two ways that users can download DNA sequence from JBrowse: in relation to a feature, like a gene, or by specific sequence coordinates.
Getting Sequence for a Gene or Transcript
To get the sequence for a gene, open any of the âCurated Genesâ tracks and right click or control click on the feature and select âView Sequenceâ from the resulting popup menu. If you are using the main âCurated Genesâ track, youâll get a dialog box asking which transcript you want to view. It does this because it needs to know what subfeatures like exons and introns can be shown. If youâre using one of the âtype-specificâ tracks (like âCurated Genes (noncoding)â), you wonât see this dialog because these tracks show individual transcripts rather than grouping the transcripts together like the âmainâ gene track does. Itâs possible (in fact, likely) that youâll be displayed a warning message about overlapping subfeatures (since exons and CDS overlap) and then be presented with a dialog box containing a list of subfeatures with buttons and the sequence below them. The buttons allow modifying the sequence display to highlight, hide or change the case of the subfeatures selected. For example, you could do an âin silico mRNA processingâ by hiding the introns. In addition, there are drop down menus for showing a set range of up or downstream sequence, up to 4kb. This sequence can then be selected, copied and pasted into word processing program with highlighting preserved. Of course, hiding or changing the case will be preserved regardless of where it is pasted into.
Two Important Notes
- The highlighting tool doesnât always get the subfeatures right! It sometimes gets confused when calculating boundaries and their overlaps. If itâs important, please take the time to check the results before using them.
- Highlighting doesnât always get preserved when cutting and pasting, and it seems to depend on what browser is being used. In my hands, it works in Firefox but not in Safari or Chrome.
Getting Sequence for Other (non-gene) Features
If there are features other than those that are gene related, the process is similar. Control or right click on a feature and select âView Detailsâ from the popup menu. The resulting dialog contains all of the information that JBrowse âknowsâ about the feature and, as long as it isnât too long (ie, less than 250kb), will contain the sequence of the feature. In this instance you donât even need to select and copy to get the sequence, since there is a âsave as FASTAâ button.
Getting Sequence of a Specific Range
If you need some other range of sequence not directly associated with a feature, we will make use of the âSave track dataâ functionality that is built into every track. In this instance, we want to use the reference sequence track to get its data (ie, the DNA sequence). To do this, turn on the âReference Sequence (DNA)â track. Then zoom, pan, search for a landmark or enter exact coordinates in the search field to the exact region that you want the sequence for. The coordinates are in the format âChromosome:start..endâ (like GBrowse, âII:20000..30000â for example). When youâre there, mouse over the track label (ie, the words âReference Sequence (DNA)â in the JBrowse view) and the track label will become more opaque and a down triangle will appear, indicating a menu is available.
Click on that triangle and select âSave track dataâ from the menu. You will be presented with a dialog menu that will let you either view or save the visible region as FASTA. There isnât an option to pick where it gets saved, so look for it in your default download location, typically a Downloads folder. A useful side note is that this option is available for every track, so you can always download the data for the track you are looking at regardless of the data type.
Dear Wormbase,
I was trying to get the sequence for a specific region by using Jbrowse, I cannot see any FASTA format for download , it was gff, sqn and bed.
Please help me if , I am missing something, there are no scroll downbars, so I think , it’s just these 3 formats.
Please let me know the solution.
Regards
Bikash