Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,632,325	C→T	E35K (GAA→AAA)	SAUSA300_RS08070 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,632,325	0	C	T	92.9%	27.9 / ‑4.4	14	E35K (GAA→AAA)	SAUSA300_RS08070	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (4/9);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCA  >  NC_007793/1632280‑1632418                                              |                                                                                              cAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGc                                         <  1:4721570/100‑1 (MQ=255)                                 tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa         >  1:2390683/7‑100 (MQ=255)                                 tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa         >  1:2946612/7‑100 (MQ=255)                                  cactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATAc        >  1:4568314/6‑100 (MQ=255)                                   actcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACa       <  1:2683281/96‑1 (MQ=255)                                    ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa      >  1:2131754/4‑100 (MQ=255)                                    ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa      <  1:5736548/97‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc   <  1:1147150/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc   <  1:2204886/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc   <  1:2720558/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc   <  1:2887236/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc   <  1:484466/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc   <  1:521497/100‑1 (MQ=255)                                        tACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCa  <  1:4389652/100‑1 (MQ=255)                                              |                                                                                              CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCA  >  NC_007793/1632280‑1632418

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