Methods and systems for detecting sequence variants
The invention provides methods for identifying rare variants near a structural variation in a genetic sequence, for example, in a nucleic acid sample taken from a subject. The invention additionally includes methods for aligning reads (e.g., nucleic acid reads) to a reference sequence construct accounting for the structural variation, methods for building a reference sequence construct accounting for the structural variation or the structural variation and the rare variant, and systems that use the alignment methods to identify rare variants. The method is scalable, and can be used to align millions of reads to a construct thousands of bases long, or longer.
1. A method for identifying a mutation in proximity to a structural variation in a sequence, the method comprising the steps of:
obtaining a plurality of nucleic acid sequence reads, wherein at least one nucleic acid read comprises a mutation;
comparing said reads to a reference sequence construct, wherein said reference sequence construct is stored in computer memory as a directed acyclic graph comprising at least two alternative sequences at a position in the reference sequence construct, one of which is a structural variation,
scoring sequence overlaps for each nucleic acid read against the reference sequence construct;
aligning each read to a location on the construct such that the score for each read is maximized; and
identifying the mutation as being aligned within 100 bp or fewer of the structural variation.
2. The method of claim 1 , further comprising assembling the nucleic acid reads to each other based upon the alignment of the nucleic acid reads with respect to the reference sequence construct.
3. The method of claim 1 , wherein the structural variation is at least 100 bp long.
4. The method of claim 1 , wherein the reference sequence construct further comprises at least two additional alternative sequences at a second position in the reference construct, and one of the additional alternative sequences comprises a sequence matching the mutation.
5. The method of claim 4 , wherein the first and second positions are separated by 100 bp or fewer.
6. The method of claim 1 , wherein the reference sequence construct further comprises at least two additional alternative sequences at a second position in the reference construct, and neither of the additional alternative sequences comprises a sequence matching the mutation.
7. The method of claim 6 , wherein the first and second positions are separated by 100 bp or fewer.
8. The method of claim 1 , wherein the structural variation is 1 kilobase to 3 megabases in length.
9. The method of claim 1 , wherein the mutation was not previously identified in a variant call format (VCF) file.
10. The method of claim 1 , wherein the mutation was previously identified in a variant call format (VCF) file.
11. The method of claim 1 , wherein the structural variation is selected from the group consisting of deletions, duplications, copy-number variations, insertions, inversions, and translocations.
12. The method of claim 1 , wherein the mutation is selected from the group consisting of a deletion, a duplication, an inversion, an insertion, and a single nucleotide polymorphism.
13. The method of claim 1 , wherein the mutation does not comprise a sequence matching the reference construct.
14. The method of claim 1 , wherein the reference sequence construct comprises a genome of an organism.
15. The method of claim 1 , wherein the reference sequence comprises a chromosome of an organism.