Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,127,723	(G)5→4	coding (18/825 nt)	SAUSA300_RS10840 →	phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,127,719	0	G	.	100.0%	29.5 / NA	9	coding (14/825 nt)	SAUSA300_RS10840	phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/0);  total (9/0)

AAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGC  >  NC_007793/2127620‑2127812                                                                                                    |                                                                                              gctacaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGAATCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTg                                                                                               >  1:3614069/4‑100 (MQ=255)    acaacTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGAATCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTgggg                                                                                            <  1:4203653/100‑1 (MQ=255)                        tATGTAATTTAGGGACCCATTAGGGAATCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACT‑GGGGCAATATAAACGCGCAGAtcg                                                                       >  1:460727/1‑98 (MQ=255)                                 caGGGACCCATTAGGGAATCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACT‑GGGGCAATATAAACGCGCTGATTTAATTGGACa                                                              >  1:3993357/2‑100 (MQ=255)                                                    ccAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACT‑TGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAaat                                           >  1:4499913/1‑100 (MQ=255)                                                    ccAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACT‑GGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAaat                                           >  1:2525830/1‑100 (MQ=255)                                                                         gttACAAGGTTTCTATGTATCCAAACT‑GGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAAt                      >  1:2906125/1‑100 (MQ=255)                                                                                    tCTATGTATCCAAACT‑GGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGa           >  1:1982455/1‑100 (MQ=255)                                                                                       aTGTATCCAAACT‑GGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATaa        >  1:3349536/1‑100 (MQ=255)                                                                                           aTCCAAACT‑GGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGt    >  1:2718612/1‑100 (MQ=255)                                                                                             ccAAACT‑GGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGc  >  1:4042438/1‑100 (MQ=255)                                                                                                    |                                                                                              AAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGACTCCAAACCCAATAAATACTGTTGTTACAAGGTTTCTATGTATCCAAACTGGGGGCAATATAAACGCGCTGATTTAATCGGACAATCTTCTTATATTAAAAATAATGATGTCGTAATATTCAATGAAGCATTTGATAATGGTGC  >  NC_007793/2127620‑2127812

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