Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,374	C→T	Q151Q (CAG→CAA)	SAUSA300_RS07730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,374	0	C	T	90.6%	68.3 / ‑2.7	32	Q151Q (CAG→CAA)	SAUSA300_RS07730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (6/23);  total (9/23)
Fisher's exact test for biased strand distribution p-value = 1.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTGGTCAGGATAAAATTCTTTGAATATTATCTGTGGTGTGC  >  NC_007793/1579281‑1579468                                                                                               |                                                                                               aataTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCAGATCg                                                                                          >  1:132888/4‑100 (MQ=255) aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTct                                                                                          <  1:6616270/97‑2 (MQ=255) aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                                                          >  1:722819/4‑100 (MQ=255) aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                                                          >  1:1652475/4‑100 (MQ=255) aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                                                          >  1:3111560/4‑100 (MQ=255) aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                                                          >  1:3280892/4‑100 (MQ=255) aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATATTGTCTTTCGTGGTTCCTTTTCATTTTCCTTTTGTATGCGTCAATGAGTTGGTCg                                                                                          >  1:6650409/4‑100 (MQ=255) aaGATGTCTTTACATGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCAGATCg                                                                                          >  1:5540304/1‑100 (MQ=255) aaGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCtagatc                                                                                          >  1:2069662/1‑95 (MQ=255)        cTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGATCGATAgtgt                                                                                   <  1:2171842/100‑4 (MQ=255)          ttACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAatat                                                                                 <  1:1606397/100‑1 (MQ=255)                 tATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTgt                                                                          <  1:5644890/100‑2 (MQ=255)                 tATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTgt                                                                          <  1:4394623/100‑2 (MQ=255)                   tCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTgtag                                                                        <  1:4759293/100‑4 (MQ=255)                        ttCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTAt                                                                   >  1:3056748/1‑100 (MQ=255)                          ccTGCGTCGGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGt                                                                 <  1:963394/100‑1 (MQ=255)                          cagacGTGTGCTCTTCCGATC‑TGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGt                                                                 <  1:3852389/96‑1 (MQ=255)                                          c‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTagctgggc                                                 <  1:3300474/100‑8 (MQ=255)                                                  tCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGcttgt                                          <  1:4233818/100‑5 (MQ=255)                                                  ggcttCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGcttgt                                          <  1:3104356/97‑5 (MQ=255)                                                  ccTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGcttgt                                          <  1:3510153/99‑5 (MQ=255)                                                          ggttCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCggg                                  <  1:487095/98‑1 (MQ=255)                                                             tCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCggggat                               <  1:1151692/100‑4 (MQ=255)                                                                                       tCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTgtggtg     <  1:3404193/100‑1 (MQ=255)                                                                                       tCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTgtggtg     <  1:6186487/100‑1 (MQ=255)                                                                                       tCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTgtggtg     <  1:1467623/100‑1 (MQ=255)                                                                                        cAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTgtggtgc    <  1:3122245/100‑2 (MQ=255)                                                                                        cAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgt    <  1:339699/100‑1 (MQ=255)                                                                                        cAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgt    <  1:4976794/100‑1 (MQ=255)                                                                                        cAATGAGTTGGTCGATAGAATATAAGGTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgt    <  1:300032/100‑1 (MQ=255)                                                                                        cAATGAGTTGGTCGATAGAATATAAGATGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgt    <  1:2504517/100‑1 (MQ=255)                                                                                          aTGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGCTTGTTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgtat  <  1:4380195/100‑3 (MQ=255)                                                                                               |                                                                                               AAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTGGTCAGGATAAAATTCTTTGAATATTATCTGTGGTGTGC  >  NC_007793/1579281‑1579468

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