Entering edit mode
10.2 years ago
Abrari
•
0
I'm using SOAP2 aligner to map some short reads to Glycine max v1.01 reference genome. The reads are paired-end soybean reads from Lam et al. 2010. SOAP2 produces two output files: the paired-end and single-end mapping. The paired-end one is significantly larger in size.
My question is, for the purpose of SNP calling (I will try some variant caller tools), which mapping result should I use? The paired-end, single-end, or both? If I have to use both, could I simply concatenate those two output files?
Any suggestion is appreciated, since I'm new to bioinformatics. Thanks.