Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	613,657	0	T	C	60.9%	‑0.1 / 17.1	23	D837D (GAT→GAC)	SAUSA300_RS02915	serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base T (8/1);  new base C (14/0);  total (22/1)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAG  >  NC_007793/613615‑613751                                           |                                                                                               ctcagactcagaCAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCg                                       >  1:3540752/1‑100 (MQ=255)       cgcagacAGCGACTCAGACTCAGACAGCGACCCAGACCCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGACAGCGCCTCAGA‑TCGGaagagc                                >  1:3590098/3‑95 (MQ=11)         caggcAGCGACTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGTAACTCAGA‑TCGGaagagcac                              >  1:5832033/1‑93 (MQ=11)         cagatAGCGACTCAGATTCAGACAGCGACTCAGATTCAGACAGCGACTCAGACTCAGATAGTGATTCAGACTCAGATAGCGACCCAGA‑TCGGaagagcac                              >  1:4945394/1‑93 (MQ=14)         cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCTAATCAGA‑TCGGaagagcac                              >  1:3944386/1‑93 (MQ=11)         cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                              >  1:4953360/1‑93 (MQ=11)         cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                              >  1:2598541/1‑93 (MQ=11)         cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                              >  1:324243/1‑93 (MQ=14)         cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCAACCCAGA‑TCGGaagagcac                              >  1:4754719/1‑93 (MQ=14)         cagacagcgactcagactcagacagcgactcagactCAGATAGCGACTCAGACTCAGACAGCGATTCAGATTCAGACAGCGACTCAGA‑TCGGaagagcac                              >  1:81406/1‑93 (MQ=11)         cagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGA‑TCGGaagagcac                              >  1:4584388/1‑93 (MQ=14)         cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGTGACTCAGA‑TCGGaagagcac                              >  1:3926645/1‑93 (MQ=14)         cagacAGTGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGACTCAGACTCAGACAGCGATTCAGACTCAGATAGCCCCTCAGA‑TCGGaagagcac                              >  1:2331947/1‑93 (MQ=14)         cagacAGCGATTCAGATTCAGACAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGTGATTCAGA‑TCGGaagagcac                              >  1:437002/1‑93 (MQ=14)         cagacAGCGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGA‑TCGGaagagcac                              >  1:513976/1‑93 (MQ=17)         cagacAGCGATTCAGATTCAGACAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGA‑TCGGaagagcac                              >  1:4727300/1‑93 (MQ=17)         cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGCCTCAGACAGCGACTCAGACTCAGATAGAGACTCAGA‑TCGGaagagcac                              >  1:1120207/1‑93 (MQ=14)         cagacAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGACCCAGA‑TCGGaagagaca                              >  1:5687576/1‑93 (MQ=14)         cagacAGCGACTCAGATTCAGATAGCGATTCAGACACAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGC‑ACTCAGA‑TCGGaagagcaca                             >  1:2755015/1‑92 (MQ=14)         cagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                              >  1:2124420/1‑93 (MQ=11)         cagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                              >  1:3854234/1‑93 (MQ=11)         cagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACCCAGA‑TCGGaagagcac                              >  1:4036641/1‑93 (MQ=11)         cagacAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGAGACTCAGA‑TCGGaagagcac                              >  1:2184887/1‑93 (MQ=11)         cagacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGCAACTCAGA‑TCGGaagagcac                              >  1:2171200/1‑93 (MQ=14)         cagacAGCGACTCAGACTCAGACAGCGACTCAGATTTAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGA‑TCGGaagagcac                              >  1:2424658/1‑93 (MQ=14)         cagacAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAACGACTCAGA‑TCGGaagagcac                              >  1:5859409/1‑93 (MQ=17)          tgacAGCGACTCAGACTCGGATAGCGATTCAGATTCAGATAGCGATTCGGATTCAGACAGTGATTCAGATTCAGACAGCGACTCAGA‑TCGGaagagcaca                             >  1:2968376/2‑92 (MQ=14)          agacagcgactcagactcagacagcgactcagactcAGATAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACCCAGA‑TCGGaagagcaca                             >  1:4996367/1‑92 (MQ=14)            acAGCGACTCAGATTCAGATAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGACTCAGA‑TCGGaagagcacacg                           >  1:2293262/1‑90 (MQ=14)                                    aCTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcagattc    <  1:3530089/100‑1 (MQ=255)                                      tcagattcagacAGCGACTCAGACTCGGATAGCGACTCAGACTCAGATAGCGATTCAGACTCAGACAGCGATTCAGATTCAGACAGCGACtcagattcag  >  1:4671529/1‑100 (MQ=11)                                           |                                                                                               CTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAG  >  NC_007793/613615‑613751

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