Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,411,453	0	T	G	62.5%	9.0 / 12.9	16	Y37D (TAT→GAT)	SAUSA300_RS12395	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base G (0/10);  total (1/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

AAAATATAAAAAAAGAAGAGCTTTAAATCATAACGCATGCCATTTATCTTTAAATGCAATGAGTTATCATTCAAACTCTTCTTATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTGGGACATT  >  NC_007793/2411371‑2411552                                                                                   |                                                                                                    aaaatataaaaAAAGAAGAGCTTTAAATCATAACGCATGCCATTTATCTTTAAATGCAATGAGTTATCATTCAAACTCTTCTTATCAACCTACATACAtt                                                                                    <  1:910056/100‑1 (MQ=255)                                       gCCATTTATCTTTAAATGCAATGAGTTATCATTCAAACTCTTCTTATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTcaaca                                              >  1:3507342/1‑100 (MQ=255)                                                           aTGAGTTATCATTCAAACTCTTCTTATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACt                          <  1:4869753/100‑1 (MQ=255)                                                                    cATTCAAACTCTTCTTATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATc                 <  1:3425244/100‑1 (MQ=255)                                                                       gtgtgctctTCCGATCTATCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCagg              <  1:150667/95‑3 (MQ=255)                                                                        tgtgctctTCCGATCTACCTATATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGc             <  1:189113/96‑1 (MQ=255)                                                                        tccgatcttctTATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGc             <  1:3598186/95‑1 (MQ=255)                                                                         gtgctctTCCGATCTACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCa            <  1:141417/97‑1 (MQ=255)                                                                          tgctctTCCGATCTACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCAt           <  1:2404160/98‑1 (MQ=255)                                                                            ctctTCTTATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTg         <  1:1353905/100‑1 (MQ=255)                                                                            ctctTCCGATCTTCCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTg         <  1:2777396/100‑1 (MQ=255)                                                                            ctctTCCGATCTTCCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTg         <  1:3650614/100‑1 (MQ=255)                                                                            ctctTCCGATCTACCTATATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTg         <  1:1210471/100‑1 (MQ=255)                                                                            ctctTCCGATCTACCTATATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTg         <  1:3939721/100‑1 (MQ=255)                                                                            ctctTCCGATCTACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTg         <  1:5098239/100‑1 (MQ=255)                                                                            ctctTCCGATCTACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTg         <  1:5332421/100‑1 (MQ=255)                                                                              cttcttATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTggg       <  1:5473990/100‑1 (MQ=255)                                                                                   tATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTGGGACAtt  >  1:3919591/1‑100 (MQ=255)                                                                                   |                                                                                                    AAAATATAAAAAAAGAAGAGCTTTAAATCATAACGCATGCCATTTATCTTTAAATGCAATGAGTTATCATTCAAACTCTTCTTATCAACCTACATACATTGACATTATACTATCATATGTTGTTTTACATATTCAACAAACATTTCAACGTTAGCACTTTTAATATCATGCATTGGGACATT  >  NC_007793/2411371‑2411552

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