Workflow from Scientific Research

Alignment file (BAM) transformations to make spliced lrRNA-seq alignments suitable for variant calling. First, GATK’s SNCR function is used to split the reads at Ns in their cigar string, such that exons become distinct reads. Second, the flagCorrection function attributes the flag of the original read to all corresponding fragment reads.

DOI