This pipeline is created to analyze allelic imbalance for F1 crosses of inbred mouse lines. It constructs individual paternal and maternal genomes, then maps the reads from RNA-seq experiments to these genomes and counts the number of reads which map to either the reference or alternate allele at each heterozygous SNP. Next it estimates allelic imbalance for individual genes summarizing information from SNPs and constructs confidence intervals if technical replicates are available.
You will need to run prepare_reference.py once for every F1 genome and also download a bunch of reference files (for instructions, please see GenomePreparation.md).
For each RNA-seq experiment, you need to get allelic counts for individual SNPs or genes. Please refer to Allelic_counts.md) for more information.
As a result of completing the previous step, you should have a file "*_processed_gene_v3.txt" containing information about number of maternal and paternal counts per gene. Next step is to estimate allelic imbalances for each gene and condition, calculate confidence intervals using technical replicates and then finally detect genes demonstrating differential allelic imbalance. To do this, please refer to the manual.
TODO: Write credits
TODO: Write license