CirSeq accurately identifies ultra-rare and low-frequency genetic variants in RNA viruses. The method uses a unique step for circularization of fragmented viral RNAs, followed by rolling-circle RT. , CirSeq corrects for mutations introduced during its amplification steps by aligning the tandem-repeat sequences with each other and excluding those reads using informatics tools. First, single-stranded RNAs are fragmented using Zn2+ and size-selected to no more than one-third of the sequencing read length. Next, they are circularized and reverse-transcribed using random primers. Rolling-circle RT is used to generate tandem-repeat cDNA strands. The first-strand cDNAs are amplified, generating double-stranded cDNAs, followed by end repair, poly(A) tailing, and adapter ligation. The cDNA libraries are ready for sequencing.