Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,210,233	1	.	G	77.3%	36.3 / 13.8	22	coding (268/648 nt)	SAUSA300_RS11250	phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base . (5/0);  new base G (17/0);  total (22/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-02

TTGCTTCATGATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AC‑CGA‑AGATC‑‑AAGAGCACATTTTATTTATTATTAATAATATGAGCTATCGCAATGGCAAAAATGATCATGTCACTTTATCTTTAGAAGGTCAAATT  >  NC_007793/2210134‑2210319                                                                                                        |                                                                                        ctGCTTCATGATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑AC‑CGA‑AGATc                                                                                          >  1:1233003/2‑100 (MQ=255)  tGCTTCATGATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCT‑AC‑CGG‑AGATcg                                                                                         >  1:4275472/1‑99 (MQ=255)       cATGATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑GC‑CGA‑AGATcggaaga                                                                                    >  1:186678/1‑94 (MQ=255)         tGATACCGTTGACGGAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑AC‑AGA‑AGATcggaagagc                                                                                  >  1:1959782/1‑92 (MQ=255)         tGATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑AC‑CGA‑AGATcggaagagc                                                                                  >  1:5003172/1‑92 (MQ=255)         tGATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑GA‑CGA‑AGATcggaaggca                                                                                  >  1:1508621/1‑92 (MQ=255)          taTACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑ACACGA‑AGATcggaagagc                                                                                  >  1:5116042/2‑92 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑AC‑TGA‑AGATcggaagagca                                                                                 >  1:2587309/1‑91 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑TC‑CGA‑AGATcggaagagca                                                                                 >  1:5093837/1‑91 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AT‑TGA‑AGATcggaagagca                                                                                 >  1:5652490/1‑91 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AC‑CGA‑AGATcggaagagca                                                                                 >  1:5707679/1‑91 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AC‑AGATAGATcggaagagc                                                                                  >  1:2097441/1‑92 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AA‑CGA‑AGATcggaagagca                                                                                 >  1:2978115/1‑91 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATC‑‑AC‑CGA‑AGATCGGAAGAGcac                                                                                >  1:2686738/1‑100 (MQ=255)          gATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTCTATCACTATCG‑AC‑TGA‑AAATcggaagagca                                                                                 >  1:2615283/1‑91 (MQ=255)          gATACCGTTGACG‑AAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCT‑AC‑CGA‑AGATCGGAAGAGcac                                                                                >  1:811001/1‑100 (MQ=255)           cgaCCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AC‑CGA‑AGATCGGAAGAGcac                                                                                >  1:2554258/3‑100 (MQ=255)           aTACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCT‑GC‑CGA‑AGATCGGAAGAGcac                                                                                >  1:5452071/1‑100 (MQ=255)           aTACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCT‑AC‑CGA‑AGATCGGAAGAGcac                                                                                >  1:877188/1‑100 (MQ=255)           aTACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑AC‑CGA‑AGATCGGAAGAGcac                                                                                >  1:3959958/1‑100 (MQ=255)           aTACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑TA‑CGA‑AGATCGGAAGAGcac                                                                                >  1:3703434/1‑100 (MQ=255)           aTACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AC‑CGA‑AGATC‑‑AAGAGCACAt                                                                              >  1:4431551/1‑100 (MQ=255)            tACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCATCC‑CGA‑AGATCGGAAGAGcac                                                                                >  1:659101/1‑100 (MQ=255)            tACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AT‑C‑A‑AGATCGGAAGAGcacac                                                                              >  1:793058/1‑99 (MQ=255)             aCCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCT‑AC‑CGA‑AGATCGGAAGAGcacac                                                                              >  1:627593/1‑99 (MQ=255)             aCCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑GC‑CGA‑AGATCGGAAGAGcacac                                                                              >  1:5870729/1‑99 (MQ=255)             aCCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AC‑CGA‑AGATC‑‑AAGAGCACAttt                                                                            >  1:4312722/1‑100 (MQ=255)             aCCGTTGACGAAAAAGTTGTAGATGCTAACAAAAAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AT‑AGA‑AGATCGGAAGAGcacac                                                                              >  1:2581172/1‑99 (MQ=255)              tccctGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATC‑‑GC‑CGG‑AGATCGGAAGAGcacacca                                                                            >  1:4974966/5‑97 (MQ=255)                  ttaCGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCG‑GC‑CGA‑AGATCGGAAGAGcacacgtctg                                                                         >  1:666858/3‑94 (MQ=255)                  tGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCATCC‑CGA‑AGATCGGAAGAGcacacgtct                                                                          >  1:2070677/1‑95 (MQ=255)                  tGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑GC‑CGA‑AGATCGGAAGAGcacacgtctg                                                                         >  1:4745725/1‑94 (MQ=255)                   gACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑ACGCGA‑AGATCGGAAGAGcacacgtctg                                                                         >  1:2852712/1‑94 (MQ=255)                                                                                   cACTATCA‑AC‑CGA‑AGATC‑‑AAGAGCACATTTTATTTATTATTAATAATATGAGCTATCGCAATGGCAAAAATGATCATGTCACTTTATCTTTAGAAGGTCa      >  1:5981419/1‑100 (MQ=255)                                                                                       aTCA‑AC‑CGA‑AGATC‑‑AAGAGCACATTTTATTTATTATTAATAATATGAGCTATCGCAATGGCAAAAATGATCATGTCACTTTATCTTTAGAAGGTCAAAtt  >  1:613638/1‑100 (MQ=255)                                                                                       aTCA‑AC‑CGA‑AGATC‑‑AAGAGCACATTTTATTTATTATTAATAATATGAGCTATCGCAATGGCAAAAATGATCATGTCACTTTATCTTTAGAAGGTCAAAtt  >  1:94772/1‑100 (MQ=255)                                                                                                        |                                                                                        TTGCTTCATGATACCGTTGACGAAAAAGTTGTAGATGCTAACAAACAATATGTTGAATTGAAGTCATTTTTATCTTCTTTATCACTATCA‑AC‑CGA‑AGATC‑‑AAGAGCACATTTTATTTATTATTAATAATATGAGCTATCGCAATGGCAAAAATGATCATGTCACTTTATCTTTAGAAGGTCAAATT  >  NC_007793/2210134‑2210319

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