Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,980,852	T→A	intergenic (‑200/+156)	SAUSA300_RS09835 ← / ← SAUSA300_RS09840	transcriptional regulator/DNA‑binding response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,980,852	0	T	A	96.4%	77.9 / ‑3.5	28	intergenic (‑200/+156)	SAUSA300_RS09835/SAUSA300_RS09840	transcriptional regulator/DNA‑binding response regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (1/26);  minor base . (0/1);  total (1/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAGAGTCAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAATTTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATCAAATTTACAAAGAAAACTGTATAATAAACCCTATCGCTATA  >  NC_007793/1980763‑1980945                                                                                          |                                                                                              aaGAGGAAATTTCACATTTGTATTAAATTTTGATCAAGTCAAAAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCAt                                                                                     <  1:2276643/100‑1 (MQ=255)  agagTCAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCAtt                                                                                    <  1:23127/100‑2 (MQ=255)  agagTCAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATa                                                                                    <  1:1352290/100‑1 (MQ=255)   acgccAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATaa                                                                                   >  1:5302230/5‑100 (MQ=255)     gTCAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATc                                                                                 <  1:395966/100‑1 (MQ=255)      tCAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAAtgg                                                                                <  1:4312273/100‑3 (MQ=255)       cAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAA‑TTTTATGCATAATCACt                                                                              <  1:1492153/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:2161013/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:3758899/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:2834802/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:2805418/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:2771082/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:4298961/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:2413967/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:1707593/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:5096143/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:1435952/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:1425030/100‑1 (MQ=255)        aaTTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACt                                                                              <  1:1290186/100‑1 (MQ=255)          tttCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACTAg                                                                            <  1:2621306/100‑1 (MQ=255)               caTTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACTAGATTTg                                                                       <  1:1602055/100‑1 (MQ=255)                                          acaaGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGatac                                            <  1:991643/100‑2 (MQ=255)                                           caaGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATc                                           <  1:1554225/100‑1 (MQ=255)                                                               taATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATCAAATTTACAAAGAAAACTGt                       <  1:2933075/100‑1 (MQ=255)                                                               taATGAAGCGAAAGTAAATGAAGTAAAATTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATCAAATTTACAAAGAAAACTGt                       <  1:4040430/100‑1 (MQ=255)                                                                                  gAAGTAAAATTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATCAAATTTACAAAGAAAACTGTATAATAAACCCTATCGCtg    <  1:5102060/100‑2 (MQ=255)                                                                                  gAAGTAAAATTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATCAAATTTACAAAGAAAACTGTATAATAAACCCTATCGCta    <  1:1401962/100‑1 (MQ=255)                                                                                    aGTAAAATTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATCAAATTTACAAAGAAAACTGTATAATAAACCCTATCGccaca  <  1:1131154/100‑5 (MQ=255)                                                                                          |                                                                                              AAGAGTCAATTTCACATTTGTATTAAATTTTGATCAAGTCAACAAGAGTTAAGTTTACTTTATAATGAAGCGAAAGTAAATGAAGTAAATTTTTATGCATAATCACTAGATTTGATAAAACTTACACTTTTATAATGATATCAAATTTACAAAGAAAACTGTATAATAAACCCTATCGCTATA  >  NC_007793/1980763‑1980945

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