Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	3,946,352	0	C	T	71.4%	4.1 / 11.3	21	noncoding (2649/2904 nt)	rrlC	23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (5/10);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTG  >  NC_000913_3_hsa_tpiA/3946203‑3946376                                                                                                                                                      |                         tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgctgc                          <  2:555306/150‑1 (MQ=32)    aTGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTcc                       <  4:770500/150‑1 (MQ=32)         ggCTCATCACATCCTGGAGCTGAAGTAGTTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTa                  <  4:682778/150‑1 (MQ=17)            tcatcgCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACga               <  1:555409/150‑1 (MQ=2)            tcatcgCATCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACga               <  1:93627/150‑1 (MQ=2)                cgcaTCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAgg           >  1:936524/3‑150 (MQ=2)                cacaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAgg           >  2:1095702/1‑150 (MQ=32)                 gcaTCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGa          <  2:744773/149‑1 (MQ=2)                 acaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGa          >  4:813992/1‑150 (MQ=32)                   aTCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGAcc        >  4:103350/1‑150 (MQ=2)                    tCCTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCg       <  1:521718/150‑1 (MQ=2)                      cTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGa     <  2:62414/150‑1 (MQ=2)                      cTGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGa     <  2:121192/150‑1 (MQ=2)                      cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa     <  3:278972/150‑1 (MQ=32)                       tGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGaa    >  3:785990/1‑149 (MQ=2)                       tGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGaa    >  4:315265/1‑149 (MQ=2)                       tGGAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGaa    >  1:781569/1‑149 (MQ=2)                        ggAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGaat   <  2:270137/150‑3 (MQ=2)                        ggAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGaat   <  2:239906/150‑3 (MQ=2)                        ggAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGaat   <  3:396870/150‑3 (MQ=2)                        ggAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGt   >  4:497820/1‑150 (MQ=32)                         gAGCTGTAGTCGGTTCCAAGGGTTGGGCTGTTCGCCCATTAAAGCGGTACGCGAGCTGGGTTCAGAACGTCGTGAGACAGTTTGGTCCCTATCTGTTGTGGGCGTTGGAAAATTGAAAAGAGCTGTTCCTAGTACGAGAGGACCGGAATg  <  3:391104/150‑1 (MQ=2)                                                                                                                                                      |                         TCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTG  >  NC_000913_3_hsa_tpiA/3946203‑3946376

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑