Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,534,963	A→T	G123G (GGA→GGT)	SAUSA300_RS07470 →	L‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,534,963	0	A	T	94.7%	38.2 / ‑4.3	19	G123G (GGA→GGT)	SAUSA300_RS07470	L‑asparaginase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (15/3);  minor base C (0/1);  total (15/4)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

TTAGGTATCGAGCAACCTGTTGTTATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGATTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAATGATGAAATTCATACGG  >  NC_007793/1534895‑1535057                                                                     |                                                                                               ctAGGTATCGAGCAACCTGTTGTTATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTg                                                                 >  1:3854415/2‑100 (MQ=255)  tAGGTATCGAGCAACCTGTTGTTATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGc                                                                >  1:2644972/1‑100 (MQ=255)   aGGTATCGAGCAACCTGTTGTTATTACTGGCGCAATGCGCTCGCCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGcc                                                               >  1:2238781/1‑100 (MQ=255)              ttccTGTTGTTATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGaaaa                                                    >  1:3453699/3‑100 (MQ=255)              caCCTGTTGTTATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGaaaa                                                    >  1:2552545/2‑100 (MQ=255)              aaCCTGTTGTTATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGaaaa                                                    >  1:2179050/1‑100 (MQ=255)                   gttgttATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGGTGCCTCTGATGAAAAGGccc                                               >  1:4141039/1‑100 (MQ=255)                          ttaCTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATaa                                        <  1:4498076/100‑1 (MQ=255)                           taCTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATaaa                                       >  1:3048082/1‑100 (MQ=255)                               ggCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCg                                   >  1:1532658/1‑100 (MQ=255)                               ggCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCg                                   >  1:4980865/1‑100 (MQ=255)                                gcgcAATGCGCTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGt                                  <  1:3888923/100‑1 (MQ=255)                                         gcTCGTCTAATGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGt                         <  1:1555343/100‑1 (MQ=255)                                                 aaTGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAatg                 >  1:3377338/1‑100 (MQ=255)                                                 aaTGAAATTGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAatg                 >  1:4941902/1‑100 (MQ=255)                                                    ctcttccgatCTGACGGCTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAatgatg              <  1:388858/91‑1 (MQ=255)                                                         tGGTTCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAATGATGAAAtt         >  1:679448/1‑100 (MQ=255)                                                          cgctCTGACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAATGATGAAATTc        >  1:2883731/4‑100 (MQ=255)                                                                gACGGTTTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAATGATGAAATTCATACgg  >  1:1439992/1‑100 (MQ=255)                                                                     |                                                                                               TTAGGTATCGAGCAACCTGTTGTTATTACTGGCGCAATGCGCTCGTCTAATGAAATTGGTTCTGACGGATTATATAATTATATTTCCGCTATTCGCGTTGCCTCTGATGAAAAGGCCCGTCATAAAGGCGTGATGGTTGTATTTAATGATGAAATTCATACGG  >  NC_007793/1534895‑1535057

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