Table of Contents:

When and Why should I submit an alignment to GenBank?

If you are submitting multiple sequences from the same locus or region, you may submit the sequences to GenBank as an alignment.

Reasons to submit sequences as an alignment:

  1. You will be given the option to use Feature Propagate to annotate features in your submission. Feature Propagate allows you to annotate just one sequence and then features are applied to the other sequences in the alignment automatically. This can be a useful tool to assist with annotation.

  2. Alignments that pass quality checks are made available in the PopSet database after processing by NCBI staff. Sequences are also retrievable in the Nucleotide database by individual Accession numbers.

Note that individual sequences in the alignment receive Accession numbers after review by NCBI staff. The alignment itself does not receive an Accession number.

Which alignment formats are accepted?

BankIt currently accepts the following alignment formats:

Currently, BankIt does not yet support submission of alignments that include accessioned sequences already present in the GenBank/ENA/DDBJ database. Contact [email protected] if you have questions about this.

How do I submit an alignment?

Start a submission in BankIt. Proceed through the forms, providing the requested information until you arrive at the Nucleotide page. On the Nucleotide page, specify that you are importing one of the accepted alignment types by selecting the Alignment radio button. Upload your alignment file in one of the acceptable formats.

What is a sequence_ID and how do I format my sequence_IDs?

The sequence_ID identifies the same sample throughout all steps of the submission. The sequence_IDs in your alignment must conform to the following rules:

What are the guidelines for each alignment format?

FASTA+GAP Format for Aligned Nucleotide Sequences

The sequence alignment software that you are using may have an option to output your alignment in the FASTA format. To align the sequences, the software may insert gaps, thereby creating the FASTA+GAP format. The gaps will only show up in the alignment, not in the individual sequence in the database. The gaps in this example are represented by the character.

Sequences in FASTA+GAP format resemble FASTA sequences. See the page on FASTA format help for instructions on formatting FASTA sequences.

The following is an example of FASTA+GAP format without source information:

>A-0V-1-A
TCACTCTTTGGCAACGACCCGTCGTCATAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-2-A
TCACTCTTTGGCAAC---GCGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-3-A
TCACTCTTTGGCAAC---GCGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-4-A
TCACTCTTTGGCAACGACCCGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-7-A
TCACTCTTTGGCAACGACCAGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

You may add source information to the definition lines so that BankIt can determine the correct organism and any other modifiers for each sequence, however it is not required. If you do not provide source information in the alignment file you will be prompted for the information with instructions on the Organism and Source Modifiers pages in BankIt.

The following is an example of FASTA+GAP format with source information:

>A-0V-1-A [organism=Gallus gallus] [clone=C]
TCACTCTTTGGCAACGACCCGTCGTCATAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-2-A [organism=Drosophila melanogaster] [strain=D]
TCACTCTTTGGCAAC---GCGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-3-A [organism=Caenorhabditis elegans] [strain=E]
TCACTCTTTGGCAAC---GCGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-4-A [organism=Rattus norvegicus] [strain=F]
TCACTCTTTGGCAACGACCCGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

>A-0V-7-A [organism=Aspergillus nidulans] [strain=G]
TCACTCTTTGGCAACGACCAGTCGTCACAATAAAGATAGAGGGGCAACTAAAGGAAGCTCTA
TTAGATACAGGAGCAGATGATACAGTATTAGAAGAAAT

If you opt to include source information with your alignment, you must include it for each sequence.

After the “>” character, is the sequence_ID. Organism name follows in brackets. Optional modifiers also follow in brackets. BankIt will not be able to correctly interpret the organism name and the source modifiers unless you correctly format them within the square brackets. For each modifier, use the value appropriate for your samples, do not copy the values present in the above example. See the list of valid source modifiers.

NEXUS Format for Aligned Nucleotide Sequences

The sequence alignment software that you are using may have an option to output your alignment in the NEXUS interleaved format.

NEXUS files can contain ? for “missing” at the 5’ and 3’ ends of sequences, as long as this parameter is properly defined within the header of the NEXUS file. BankIt will replace the “?” characters in the sequences with “N”s since they are defined as “missing” data in the header. Gaps in the alignment are represented by the - character, as specified in the header of the NEXUS file. The gaps will only show up in the alignment, not in the individual sequence in the database.

The following is an example of NEXUS Interleaved format.

#NEXUS

begin data;
   dimensions ntax=5 nchar=100;
   format datatype=dna missing=? gap=- interleave;
   matrix

A-0V-1-A   ????TCTTTG GCAACGACCC GTCGTCATAA TAAAGATAGA GGGGCAACTA
A-0V-2-A   TCACTCTTTG GCAAC---GC GTCGTCACAA TAAAGATAGA GGGGCAACTA
A-0V-3-A   TCACTCTTTG GCAAC---GC GTCGTCACAA TAAAGATAGA GGGGCAACTA
A-0V-4-A   TCACTCTTTG GCAACGACCC GTCGTCACAA T????ATAGA GGGGCAACTA
A-0V-7-A   TCACTCTTTG GCAACGACCA GTCGTCACAA TAAAGATAGA GGGGCAACTA


A-0V-1-A   AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
A-0V-2-A   AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
A-0V-3-A   AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
A-0V-4-A   AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
A-0V-7-A   AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAG????
;
End;

In this example, the first few lines provide information about the data in the sequence alignment. The following five lines contain the Sequence IDs, followed by the sequences. In this example, the sequence_ID for the first sequence is A-0V-1-A. Note that subsequent blocks of sequence also contain the sequence_ID.

You may add the organism names and source modifiers to the alignment as shown in the example, however it is not required. If you do not provide source information in the alignment file you will be prompted for the information with instructions on the Organism and Source Modifiers pages in BankIt. The following is an example of NEXUS with source information optionally added to the end of the file.

NEXUS Interleaved EXAMPLE with SOURCE information

#NEXUS

begin data;
        dimensions  ntax=3 nchar=100;
        format datatype=dna  missing=? gap=-  interleave ;
        matrix

[     1                                                   50]
ABC_1 ???ATTGCGT TATGGAAATT CGAAACTGCC AAATACTATG TCACCATCAT
ABC_2 GATATTGCTT TATGGAAATT CGAAACTGCC AAATACTATG TCACCATCAT
ABC_3 ???ATTGCTT TATGGAAATT CGAAACTGCC AAATACTATG TTA-------

[     51                                                 100]
ABC_1 TGATGCACCT GGACACAGAG ATTTCATCAA GAACATGATC ACTGGTACTT
ABC_2 TGATGCACCT GGACACAGAA ATTTCATCAA GAACATGATC ACTGGTACTT
ABC_3 TGATGCACCT GGACACAGAG ATTTCATCAA AAACATGATC ACTGGTACTT
;
END;

begin ncbi;
sequin
>[organism=Saccharomyces cerevisiae][strain=ABC][clone=1]
>[organism=Saccharomyces cerevisiae][strain=ABC][clone=2]
>[organism=Saccharomyces cerevisiae][strain=ABC][clone=3]
;
end;

If you opt to include source information in a nexus file, source information must be included for all sequences in the alignment and it must be formatted as shown in the example with “begin ncbi;”, “sequin”, source information, “;” and “end;”.

The source information begins with a “>” character. The first line of source information applies to the first sequence (ABC_1), the second line to the second sequence (ABC_2), and so on. If you opt to include source information with your alignment, you must have one line of source information for each sequence. These inserted lines contain modifiers formatted like in the FASTA definition line, but do not begin with the sequence_ID. Instead, the sequence_ID is present at the beginning of the sequence lines as shown in the example.

After the “>” character, the organism name follows in brackets. Optional modifiers also follow in brackets. BankIt will not be able to correctly interpret the organism name and the source modifiers unless you correctly format them within the square brackets. For each modifier, use the value appropriate for your samples, do not copy the values present in the above example. See the list of valid source modifiers.

PHYLIP Format for Aligned Nucleotide Sequences

The sequence alignment software that you are using may have an option to output your alignment in the PHYLIP format.

The following is an example of PHYLIP format:

     5    100
A-0V-1-A   TCACTCTTTG GCAACGACCC GTCGTCATAA TAAAGATAGA GGGGCAACTA
A-0V-2-A   TCACTCTTTG GCAAC---GC GTCGTCACAA TAAAGATAGA GGGGCAACTA
A-0V-3-A   TCACTCTTTG GCAAC---GC GTCGTCACAA TAAAGATAGA GGGGCAACTA
A-0V-4-A   TCACTCTTTG GCAACGACCC GTCGTCACAA TAAAGATAGA GGGGCAACTA
A-0V-7-A   ----TCTTTG GCAACGACCA GTCGTCACAA TAAAGATAGA GGGGCAACTA


           AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
           AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
           AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
           AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT
           AAGGAAGCTC TATTAGATAC AGGAGCAGAT GATACAGTAT TAGAAGAAAT

In this example, the first line indicates that there are 5 sequences, each with 100 nt of sequence. The following five lines contain the sequence_IDs, followed by the sequences. In this example, the sequence_ID for the first sequence is A-0V-1-A. Note that subsequent blocks of sequence do not contain the Sequence ID. Gaps in the alignment are represented by the - character. The gaps will only show up in the alignment, not in the individual sequence in the database.

You may add the organism names and source modifiers to the alignment as shown in the example, however it is not required. If you do not provide source information in the alignment file you will be prompted for the information with instructions on the Organism and Source Modifiers pages in BankIt. The following is an example of PHYLIP with source information optionally added to the end of the file.

The following is an example of PHYLIP format with SOURCE information:

      3  100
ABC-1      ---ATTGCGT TATGGAAATT CGAAACTGCC AAATACTATG TCACCATCAT
ABC-2      GATATTGCTT TATGGAAATT CGAAACTGCC AAATACTATG TCACCATCAT
ABC-3      ---ATTGCTT TATGGAAATT CGAAACTGCC AAATACTATG TTA-------

           TGATGCACCT GGACACAGAG ATTTCATCAA GAACATGATC ACTGGTACTT
           TGATGCACCT GGACACAGAA ATTTCATCAA GAACATGATC ACTGGTACTT
           TGATGCACCT GGACACAGAG ATTTCATCAA AAACATGATC ACTGGTACTT

>[organism=Saccharomyces cerevisiae][strain=ABC][clone=1]
>[organism=Saccharomyces cerevisiae][strain=ABC][clone=2]
>[organism=Saccharomyces cerevisiae][strain=ABC][clone=3]

If you opt to include source information with your alignment, you must include it for each sequence.

The first line of the source information begins with a “>” character. The first line of source information applies to the first sequence (ABC-1), the second line to the second sequence (ABC-2), and so on. If you opt to include source information with your alignment, you must have one line of source information for each sequence. These inserted lines contain modifiers formatted like in the FASTA definition line, but do not begin with a sequence_ID. Instead, the sequence_ID is present at the beginning of the sequence lines as shown above.

After the “>” character, the organism name follows in brackets. Optional modifiers also follow in brackets. BankIt will not be able to correctly interpret the organism name and the source modifiers unless you correctly format them within the square brackets. For each modifier, use the value appropriate for your samples, do not copy the values present in the above example. See the list of valid source modifiers.

Clustal(w) Format for Aligned Nucleotide Sequences

The sequence alignment software that you are using may have an option to output your alignment in the Clustal(w) format.

The following is an example of Clustal(w) format:

CLUSTAL W

seq1	------ACTAGACTGGGTGCGGTACCTAAGTG-TTACTGGCGGTGTGTTGCTCTATGATT
seq2	CCCCTTTCTGGAAAAGTTCACGGTACTATTCG-TTTCTGCCTGTGAGCTGCTATACGATT
seq3	------ATTAGGCTAG--ACGGTACCAATAGTCGCGCTGACTGTGGGTTGCTCTACGACT
seq4	-----AATTAGGCTGGTGTCAGTATCAAG--GTTCCCTGGATGTTAGTTGCTCTACGACT
seq5	-----AAATAGGCTGGTGTCAGTATCAAGAGGATCCTGGAATGTTAGATGCTCTAAGACT
	        *.*.:: .   ::*:: *        . ::* ::**. * ****:*  **.*

seq1	CGGAATGTTACCAGGATGATATACCTAGTTGCCTAGATGCACACCTTGAATTTGTCGAAA
seq2	CGGCATGTTCCCCGGTTGATGTTGCTAGTTGCAAAGTTGCAAAGACTCAATTTGCTGAAA
seq3	CGGTAAGATCCCTGGTTGGTTGTACTCGTTGCAAATAGGCACCCATTGTAATAGTTGCAA
seq4	AGGTAAGTTCCCTGGGTGGTGTAACTCGTTGCCTATGCGCAACTACTTTAGTGGAAGAAA
seq5	AGGTAAGTTCCCTGGGTGGTATAACTCGTTGCCTAGGTGCACCGAATCCAGTGGATGAAA
	:** *:* *:**  *.**.* .: **.*****: *. .*** :.: *  * * *  *:**

seq1	CACCCTATTTTCGGGGTATGGGTGCAGGCCAGGAAGTAG--
seq2	CACCCTATTTTCGCGGTATGGTTGCAGACCAGGAAGTAGGC
seq3	CTTCCTATTTTCGGGGCATGGTTGCAGACTAGGCAGATGC-
seq4	CTACCTAGTTTCGCGGCAAGTGTGCAGGCTAGGCAGAAG--
seq5	CTACCTAGATTCGGGGCAAGTGTGCAGGCTAGGCAGAAA--
	*::*** . ****:** * *..*****.*:***:**:

In this example, there are 5 sequences in the alignment. Each line contains the sequence_ID followed by the sequence for that sequence_ID. For example, the first sequence is seq1. Each subsequent block of sequence contains the sequence_IDs. Gaps in the alignment are represented by the - character. The gaps will only show up in the alignment, not in the individual sequence in the database. After each alignment block are the sequence conservation characters included in the Clustal(w) output.

At this time, you may not modify the Clustal(w) output to include source information. You will be prompted for source information in the BankIt forms as you continue with your submission.