Mopepgen

Fixed

- Fixed that reference source are not recognized. Switched to use upper case values. 758

- Fixed `generateIndex` that symlink was not created properly for GTF with `--gtf-symlink`.

- Fixed matplotlib warning message 742

- Fixed `parseVEP` that insertions not parsed successfully if the location is a single base. 766

- Fixed `callVariant` that peptides with upstream cleavage altering mutations were not called. 670

- Fixed fuzzTest to take multiple CPUs and use temporary directory.

- Fixed issue in callVariant of circRNA with a SNV being silent in the first loop but not in the second.

- Fixed issue that circRNA peptide nodes with and without a silent mutation were collapsed. 778

- Fixed that circRNA peptides that carry indels incompatible with the orf start node should not be called. 780

- Fixed that hybrid node was not identified if the variant is in the end of an exon. 782

- Fixed that in circRNA, cleavage gain from upstream node is added to the the wrong ORF. 783

- Fixed callVariant that `fit_into_codon` terminated early in fusion transcripts when there donor has a frameshift and accepter has a variant right after the breakpoint. 786

Fixed

- Reduced memory usage by mapping the GTF files into memory instead of reading it all at once. 371

Added

- Added the support for calling peptides of Selenocysteine terminated. 684

- Added the support for calling peptides of W > F codon reassignments. 484

- Added the support for calling peptides with adjacent SNP/INDEL. 691

- Updated fake.py and bruteForce to handle selenocysteine, W2F and MNV. 689

- `callAltTranslation` added to call peptides with alternative translation without any genomic or transcriptomic variations.

- Enabled `summarizeFasta` to create bar plot of the summary results.

Fixed

- Fixed `fake` that simulated selenocysteine positions could be in introns.

- Fixed `fake` that the last exon was picked for A3SS or first exon for A5SS.

- Fixed fusion with very small intronic insertion. 707

- In ThreeFrameTVG when aligning variant bubbles and when nodes are merged, variants were not merged correctly.

- Fixed TVG that indel merged with downstream fusion treated as subgraph out. 708

- Fixed `parseRMATS` to handle more complex situations such as exons interjacent between splicing sites and exons spanning over the splicing site. 715, 716, 717, and PR 720

- Fixed `callVariant` that failed when there is a SNV very close to the end on a AltSplice insertion. 723

- Fixed `TranscriptAnnotationModel` for not recognizing transcripts with `mRNA_end_NF` correctly. 724

- Fixed `callVariant` issue of altSplice insertion carries an intronic indel that goes back to the original reading frame. 726

- Fixed `callVariant` to handle deletion that spans over an entire intron. 732

- Fixed `callVariant` to skip peptides earlier if they are either too long or too short to significantly improve efficiency. 736

- Fixed `callVariant` to handle hypermutated region with a dynamic cutoff. 738

- Fixed `decoyFasta` to make it as default to keep cleavage site amino acid residues unmodified. 750

- Fixed `SeqFeature` and `GTFSeqFeature` to remove the definition of strand and use `location.strand`. 616

- Refactored `util` so all functions are accessible. 749

Fixed

- Fixed `callVariant` that the command line argument `--max-variants-per-node` and `--additional-variants-per-misc` not passed to it and the default value was always used.

Fixed

- Fixed `callVariant` that variant peptides may get redundant labels with same information (transcript ID and variants). 679

Fixed

- `filterFasta` failed with the new noncoding peptide FASTA header. 675