Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	503,675	0	A	T	77.8%	11.0 / ‑1.0	9	noncoding (25/93 nt)	SAUSA300_RS02390	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (7/0);  total (8/1)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TTGTGTTTAGGAGGAATACCCAAGTCCGGCTGAAGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGTTCGAATCCCTCTTCCTCCGCCATCAATATTTATATTAAATTCTATA  >  NC_007793/503642‑503767                                  |                                                                                             ttGTGTTTAGGAGGAATACCCAAGTCCGGCTGAAGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGTTCGAAACCCTCTtcctcccc                            >  1:4360514/1‑98 (MQ=255) cagtgtTTAGGAGGAATACCCAAGTCCGGCTGATGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGGTCGAATCCCTCTTCCTCCGc                            >  1:3768504/3‑100 (MQ=255)  ggtgtTTAGGAGGAATACCCAAGTCCGGCTGAAGGGACTGGTCTTGAAAACGGACAGGGGCTTAACGGTTCGCGGGGGGTCGAATCCCACCTCCTTCGcc                           >  1:4662/2‑100 (MQ=14)   gtgtTTAGGAGGAATACCCAAGTCCGGCTGATGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGGTCGAAACCCTCCTCCTCCGCCa                          >  1:4477255/1‑100 (MQ=14)       atAGGAGGAATACCCAAGTCCGGCTGATGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGGTCGAATCCCTCTTCCTCCGCCATCaa                      >  1:4586561/2‑100 (MQ=255)          ggaggaATACCCAAGTCCGGCTGATGGGATCGGTCTTGAGAACCGACAGGGGCTTAACGGCTCGGGGGGGTTTGAATCCCTCTTCCTCCGCCATCAattt                   >  1:2082114/1‑98 (MQ=14)                  aCCCAAGTCCGGCTGATGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGTTCGAATCCCTCTTCCTCCGCCATCAATATTtaacaga           >  1:578687/1‑95 (MQ=255)                  aCCCAAGTCCGGCTGATGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGGTCGAATCCCTCTTCCTCCGCCATCAATATTTATATTa           >  1:2198886/1‑100 (MQ=255)                   cccAAGTCCGGCTGAAGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGTTCGAATCCCTCTTCCTCCGCCATCAATATTtaacgtac          <  1:226180/100‑7 (MQ=255)                           cGGCTGATGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGGTCGAATCCCTCTTCCTCCGCCCTCAATATTTATATTAAATTCtata  >  1:821151/1‑100 (MQ=255)                                  |                                                                                             TTGTGTTTAGGAGGAATACCCAAGTCCGGCTGAAGGGATCGGTCTTGAAAACCGACAGGGGCTTAACGGCTCGCGGGGGTTCGAATCCCTCTTCCTCCGCCATCAATATTTATATTAAATTCTATA  >  NC_007793/503642‑503767

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