Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,586,131	C→T	E149K (GAA→AAA)	SAUSA300_RS13290 ←	peptidase M28

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,586,131	0	C	T	100.0%	48.4 / NA	19	E149K (GAA→AAA)	SAUSA300_RS13290	peptidase M28
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (1/18);  total (1/18)

GCTTTCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTCAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGCGGTTT  >  NC_007793/2586043‑2586227                                                                                         |                                                                                                 gCTTTCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCt                                                                                       <  1:6250641/100‑1 (MQ=255) gCTTTCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCt                                                                                       <  1:5508404/100‑1 (MQ=255)     tCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGct                                                                                   <  1:1108974/100‑1 (MQ=255)     tCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGca                                                                                   <  1:337262/100‑2 (MQ=255)     tCGCAATTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGct                                                                                   <  1:585876/100‑1 (MQ=255)      cGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCCTCTTCCTTGcca                                                                                  <  1:6028194/100‑3 (MQ=255)       gCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCtggtagt                                                                                 <  1:2199556/100‑7 (MQ=255)        cAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAGCCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctact                                                                                <  1:3036982/100‑1 (MQ=255)        cAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCtccagcgc                                                                                <  1:2851181/100‑8 (MQ=255)         aGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctagcg                                                                               <  1:513066/100‑4 (MQ=255)         aGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctacta                                                                               <  1:3482526/100‑1 (MQ=255)         aGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctacta                                                                               <  1:4149484/100‑1 (MQ=255)              aaTATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACAcc                                                                          <  1:6346299/100‑1 (MQ=255)                      tCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTCCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCta                                                                  <  1:4110957/100‑1 (MQ=255)                                     tCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGa                                                   >  1:5974001/1‑100 (MQ=255)                                               tataTGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTgg                                         <  1:1630556/100‑1 (MQ=255)                                                              tatTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACgggttgt                          <  1:464579/100‑7 (MQ=255)                                                                       ccgatctgCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTat                 <  1:4363646/93‑2 (MQ=255)                                                                              gCCAGCTTCTTTAGCTTCTTCCTTACTACTATCACCAATAGCTATAAACATATTTTTGAATTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACagag          <  1:296118/100‑4 (MQ=255)                                                                                    ttctttagcttcttcCTTGCTACTAACACCAATATCTAAAAACAGATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGCGGt    <  1:5234/100‑1 (MQ=255)                                                                                      ctttagcttcttcCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGCGGttt  <  1:195210/100‑1 (MQ=255)                                                                                      ctttagcttcttcCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTCTACGTTCTTCAGGCGTTAAGACATGTGGCGGttt  <  1:5018372/100‑1 (MQ=255)                                                                                         |                                                                                                 GCTTTCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTCAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGCGGTTT  >  NC_007793/2586043‑2586227

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