The 5’end of a eukaryotic mRNA generally has a methyl guanosine cap (m7G cap) that not only protects the mRNA from decay by 5'-3’exonucleases, but also plays an essential role in almost all aspects of gene expression. Some RNAs in E. coli, yeast, and mammals were recently found to contain an NAD+ cap at their 5’ends. We have developed a new method – NAD tagSeq – for transcriptome-wide identification and quantification of NAD+-capped RNAs (NAD-RNAs). The method uses an enzymatic reaction and a click chemistry reaction to label NAD-RNAs with a synthetic RNA tag. The tagged RNA molecules can be enriched and directly sequenced using the Oxford Nanopore sequencing technology. NAD tagSeq not only allows more accurate identification and quantification of NAD-RNAs but can also reveal sequences of whole NAD-RNA transcripts. Using NAD tagSeq, we found that NAD-RNAs in Arabidopsis are mostly produced from a few thousand protein-coding genes. The top 2,000 genes that were found to produce the highest numbers of NAD-RNAs were enriched in the gene ontology terms of responses to stresses, photosynthesis, and protein synthesis. For some Arabidopsis genes, over 10% of their transcripts could be NAD-capped. The NAD-RNAs in Arabidopsis have similar overall sequence structures to their canonical m7G-capped mRNAs. The identification and quantification of NAD-RNAs and revealing their sequence features provide essential steps toward understanding functions of NAD-RNAs.