Incorrect annotation of exon/intron boundary deletion

Tái hiện báo cáo bằng VariantAnnotation_1.49.1, cơ sở dữ liệu transcript AnnotationHub được trích dẫn, gói genome GRCh38, transcript ENST00000378230 và deletion trên chr1. Bắt đầu với predictCoding và kiểm tra cách nó xử lý ranh giới exon 16 và intron liền kề. Được xem là hoàn thành khi các codon được báo cáo và vị trí CDS không còn kéo dài vào exon hoặc intron không đúng đối với biến thể này.

Using:

• VariantAnnotation_1.49.1
• txdb = AnnotationHub::AnnotationHub()[["AH100643"]] (UCSC seq styles)
• genome = BSgenome.Hsapiens.NCBI.GRCh38::BSgenome.Hsapiens.NCBI.GRCh38 (UCSC seq styles)

I've got the following variant:

chr1:3826659-3826710 (genomic)
"CGATTCTTTACACACCCCAGTTCTTTGTGCACCCCAATTCTTTACATACCCT" (ref) -> "C" (alt)
                                                ( ^ exon 16 end <--- )

This overlaps transcript ENST00000378230, which has the following exons:

seqnames          ranges strand |         exon_id exon_rank
   <Rle>       <IRanges>  <Rle> |     <character> <integer>
    chr1 3829266-3829373      - | ENSE00001690217        15
    chr1 3826708-3826744      - | ENSE00001730754        16
    chr1 3826370-3826436      - | ENSE00001769704        17

so the deletion only partially affects exon 16, and then the intron after. The predictCoding function reports:

REFCODON: CAGGACATTCAAGGAGGGAAAGCAGCCCCTGCTGAAGCTCTGGGAATCCCGGAT (= QDIQGGKAAPAEALGIPD)
VARCODON: CGT (= R; GT from extending into the intron?)
CDSLOC: 2186-2188

The sequence of each exon listed above is:

15: GCACGGAGAAAAGCGGCTACAGAAGAAGCAGAAAAACAAAAGAAAGAAGAAATAAAAGCTTTACAAGGGCAGCTGGCAGCACTGAAAGAAATTCAGGCTGAAGTTCAG
16: GAAAAAGAAAGTGATGCTGTGAAGCCAAAGAATCAGG [AGG = revcomp of variant end]
                                     ^ refcodon start
17: ACATTCAAGGAGGGAAAGCAGCCCCTGCTGAAGCTCTGGGAATCCCGGATGAGCACTATCTAGATAA
                                        refcodon end ^ (but this exon should not be altered)

So this looks like VariantAnnotation incorrectly extends the REFCODON from one exon to the next, while the VARCODON extens into the (deleted part of the) intron?