Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,270,478	1	.	G	22.7%	50.4 / 11.2	22	noncoding (115/377 nt)	rnpB	RNase P catalytic RNA component
Reads supporting (aligned to +/- strand):  ref base . (5/12);  new base G (5/0);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 9.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCT  >  NC_000913/3270228‑3270727                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                             ccGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGTTCTC‑TGTTGCTCTCGTCGTggg                                                                                                                                                                                                                                                               >  1:79592/1‑251 (MQ=255)                      gggACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCt                                                                                                                                                                                                                                         <  1:80920/250‑1 (MQ=255)                       ggACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTTTC‑TGGTGCCCCTGTCGGGGGGTTTcccccccggcggtatctg                                                                                                                                                                                                                                         >  1:236962/1‑240 (MQ=255)                            gTCATTCATCTAGGCCAGCACTCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTTCTC‑TTTTTCACTGGTCGTGGGGTTTCCCCCCC‑GGCGTTACCttgggcc                                                                                                                                                                                                                                   >  1:157875/1‑245 (MQ=255)                                                        cTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCCTGTTGCCCTG‑TCCTGGG‑GTTCCCCCCCCAGCGCTACCTGGCTACCC‑CCCCTCTGGGGCCC‑GGGCTTtcctcc                                                                                                                                                                                                     >  1:144622/1‑251 (MQ=255)                                                          ggCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCTGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTtcctcc                                                                                                                                                                                                     <  1:249163/249‑1 (MQ=255)                                                                     ccTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCTTTCTCTC‑TGTTTCCCCGGTCGGGGGGTTCCCCCCCCGGGCGTTCCCTGTC‑‑CCCTTCCCTATTGGCCCC‑GGGCTTTtctcccccccccccctcc                                                                                                                                                                                        >  1:174384/1‑240 (MQ=255)                                                                                  aGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTTTC‑TTTTTCACCGGTCGTGGG‑GTTCTCCCCCCCGCGGTTCCCGGC‑GCCCCGCCCCATTGTGCCCCGGCCTTTCC‑CCCCTCCCCCCGCCTccccccccagggg                                                                                                                                                                           >  1:262808/1‑243 (MQ=255)                                                                                       gggCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTC‑TGTTGCCCCGGTCGGGGGGTTTCCCCCCC‑GGCGTTTCCCGGC‑CCCCCGCCCTCTTGGGCCC‑GGGCTTTCCTCCCCCCCCCCCCCCCCCCCCGCAGA‑Ggacgacc                                                                                                                                                                     >  1:274486/1‑250 (MQ=255)                                                                                              cccTTATGAACCCCTATTTGGCCTTGCTCCGGTTGGAGTTTACCGTGCCACGGACTGTTCCCACCCGCGCGGTGCGCTCTTACCGCACCCTTTCCCCCTGTCCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCCCCCCACCGCCCGTCTCCCCCGACGA‑GGACGACGACGGAgc                                                                                                                                                              <  2:300365/250‑1 (MQ=255)                                                                                                                                ggAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCTCCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGtata                                                                                                                            <  2:111833/251‑1 (MQ=255)                                                                                                                                     ttACCGTGCCACGGACTGTTACCCGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTAGGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAgg                                                                                                                        <  2:79592/250‑1 (MQ=255)                                                                                                                                                           cAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑Ggacg                                                                                                                                                                        <  1:86297/180‑1 (MQ=255)                                                                                                                                                              ccGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTTCTCTC‑TGTTGCACCGGTCGTGGGGTTTCCCCCCC‑GACGTTAACTGGC‑ACCCTGCCCTATGGAGCCC‑CGACCTTTCTCCCCCCCGCCCGGCTCCCCCGCAGAAGGA‑GAGGACGAAACGGCGCGCGGCTGGTCGGCTTCGGCGCGCAGATTAGAGGGTTTTGGCGCCCTTTTCAccccccgc                                                                                              >  1:385603/1‑247 (MQ=255)                                                                                                                                                                                          ccccccTTCCCTGATCCCGCTTCCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCgcggc                                                                  <  2:262808/248‑1 (MQ=255)                                                                                                                                                                                           cACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCgcggcg                                                                 <  1:375627/251‑1 (MQ=255)                                                                                                                                                                                                               ttcgcgGGCCCTCGGCGGTTTGCTCTC‑TGTTGCCCGGGTCGTGGG‑GTTCCCCCCCAGGCGTTACCTGGC‑CCCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACgg                                               <  2:88198/247‑1 (MQ=255)                                                                                                                                                                                                                      cccaccggcggTTTGCTCCC‑TTTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACGGCAGATATTg                                      <  2:274486/246‑1 (MQ=255)                                                                                                                                                                                                                           cggcggTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGGTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACGGCAGATATTGTCAt                                  >  1:337667/1‑250 (MQ=255)                                                                                                                                                                                                                                   tGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACGGCAGATATTGTCATAAGCCCcg                          <  1:179750/250‑1 (MQ=255)                                                                                                                                                                                                                                              tGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGatat                                                             <  1:396458/205‑1 (MQ=255)                                                                                                                                                                                                                                              tGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGatat                                                             >  2:396458/1‑205 (MQ=255)                                                                                                                                                                                                                                                cACTGGCCGTGGG‑TTTCCCACCCGGGCGTTACCTGCC‑ACCCTGCCCTATGCAGCCC‑GGACTTTCATCCCCTCCGCCCGTCTCCACCGCAGA‑GGACCACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATCGAGGATTTGCGCGCCCTTGTCCCACCGCGTTGCGCATTCCCATCGCCAGTGTCGCGGCGATATTACGTGACGCACGGCCGATCTTGTAATACACCCCGCGGAAtacctcgg             >  2:301126/1‑244 (MQ=255)                                                                                                                                                                                                                                                          ggg‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCGCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATGACGGGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTAc    <  2:254648/250‑1 (MQ=255)                                                                                                                                                                                                                                                           gg‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCt  >  2:11379/1‑251 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                             CCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTC‑TGTTGCACTGGTCGTGGG‑TTTCCCCCCCAGGCGTTACCTGGC‑ACCCTGCCCTATGGAGCCC‑GGACTTTCCTCCCCTCCGCCCGTCTCCCCCGAAGA‑GGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGCAGTATAGAGGGTTTGCGCGCCCTTGTCACCCCGCGTTGCGCATTCCAATCGCCAGTGTCGCGGCGATATTACGTGAAGCACGGCAGATATTGTCATAAGCCCCGCGGAATGCTTCGTCCAACGTACCT  >  NC_000913/3270228‑3270727

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑