Question: meaning of XM:i:N in Bowtie output
1
gravatar for Don
3.2 years ago by
Don10
Florida.USA
Don10 wrote:

Hi, I am a newbie for NGS, I was experimenting with chip-seq analysis with galaxy-bowtie. Question is: what does the XM:i:N mean in the bowtie file? Reason I ask is: My research shows it means (from bowtie manual) : "For a read with no reported alignments, <N> is 0 if the read had no alignments. If [-m] was specified and the read's alignments were supressed because the [-m] ceiling was exceeded, <N> equals the [-m] ceiling + 1, to indicate that there were at least that many valid alignments (but all were suppressed). In [-M] mode, if the alignment was randomly selected because the [-M] ceiling was exceeded, <N> equals the [-M] ceiling + 1, to indicate that there were at least that many valid alignments (of which one was reported at random)." I've pasted few lines from the bowtie output. I noticed when there is a unique match, the XM:I field is not shown. When XM:i:1 or XM:i:0, there is no match. Problem is shouldn't XM:i:1 mean there is a match???? Here is the bowtie parameters: ID:Bowtie VN:0.12.7 CL:"bowtie -q -p 6 -S -n 2 -e 70 -l 28 --maxbts 125 -k 1 -m 1 --phred33-quals /galaxy/data/hg19/hg19full/bowtie_index/hg19full /galaxy-repl/main/files/019/148/dataset_19148014.dat"

QNAME   FLAG    RNAME   POS MAPQ    CIGAR   MRNM    MPOS    ISIZE   SEQ QUAL    OPT
SRR4237722.115 HS38_18889:1:1101:10918:1975 length=50   4   *   0   0   *   *   0   0   GGGTCNCACTCTGTCACCCAGGCTGGAGTGCAGTGGCCGTGATCTCGGCT  BCCDF!2AAFFHHIJJIJIJFHIIDGIGFCEGGHIIJGIDHIIJIEEHHH  XM:i:1
SRR4237722.108 HS38_18889:1:1101:10120:1993 length=50   4   *   0   0   *   *   0   0   ATTCCATTCCACTCCATTCCATTCCAATCCTTTCCGCTCGGGTTGATTCC  CCCFFEFFGHHHGJIJJIIIIIJJJJJJIJJJIJJJJJJJII?FHIIJJJ  XM:i:0
SRR4237722.168  0   chr13   75065209    255 50M *   0   0   CACTGGAGCCTCTGATATCTTTTTTCTTTTTCCTGCTGACACAGGACTGC  CCCFFFDEHHHHHIIJJJJJJJJJJIJIGGJIJJIIJJJJJIJJJIJIJJ  XA:i:0 MD:Z:50 NM:i:0
SRR4237722.94   0   chr3    121413904   255 50M *   0   0   TGATGNAACCGATTCTGAACATGTAGGTCTTGTGCTCATACTCAGAGAGT  @@?DD!2=CFHHHFHIIF>GHIIHCGI*@FGHFGIIIIHIIHGIIIGIID  XA:i:1 MD:Z:5G44 NM:i:1
chip-seq • 3.4k views
ADD COMMENTlink modified 3.2 years ago by Sej Modha4.7k • written 3.2 years ago by Don10

I will really appreciate any help.

ADD REPLYlink written 3.2 years ago by Don10
1
gravatar for Sej Modha
3.2 years ago by
Sej Modha4.7k
Glasgow, UK
Sej Modha4.7k wrote:

You can check details about these tags here:

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by Sej Modha4.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1802 users visited in the last hour