Hey guys, I am using the command line "samtools pileup -cv -f ref_genome snp.sorted.bam" trying to get SNP information between RNA-Seq data and reference genome. The output file has ten columns, which I've no idea for most of them. Here are some of the output:
Group1.1 225532 T G 0 3 60 1 g B
Group1.2 1381414 C Y 0 3 60 1 T$ B
Group1.3 589260 A R 0 3 60 1 ^~g B
From what I can tell so far, the 1st column is the scaffold id and the 2nd is the scaffold coordinate where SNP occurs between the scaffold and the RNA-Seq read. Anyone has previous experiences? Thanks so much!
Group1.1 225532 T G 0 3 60 1 g B
Group1.2 1381414 C Y 0 3 60 1 T$ B
Group1.3 589260 A R 0 3 60 1 ^~g B
From what I can tell so far, the 1st column is the scaffold id and the 2nd is the scaffold coordinate where SNP occurs between the scaffold and the RNA-Seq read. Anyone has previous experiences? Thanks so much!
Comment