Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,043	T→C	K493R (AAA→AGA)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,043	0	T	C	100.0%	68.7 / NA	25	K493R (AAA→AGA)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/23);  total (2/23)

GATTTCGTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTTTAATATCAATGCGATCAATCACTGATAGATACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAA  >  NC_007793/45985‑46118                                                           |                                                                            gATTTCGTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTg                                    <  1:4011933/100‑1 (MQ=255)     ttgTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAgat                                >  1:499406/3‑98 (MQ=255)      cGTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCg                               <  1:5712385/100‑1 (MQ=255)      cGTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCg                               <  1:2008559/100‑1 (MQ=255)      cGTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCg                               <  1:2218565/100‑1 (MQ=255)      cGTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCg                               <  1:4831558/100‑1 (MQ=255)                       gtgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCt              >  1:4739178/1‑100 (MQ=255)                        tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGGGATTTTTCCATGTGCtt             <  1:695619/100‑1 (MQ=255)                        tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:5743295/100‑1 (MQ=255)                        tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:1522286/100‑1 (MQ=255)                        tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:239737/100‑1 (MQ=255)                        tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:466737/100‑1 (MQ=255)                        tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtg             <  1:808381/100‑2 (MQ=255)                        tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCta             <  1:102945/100‑2 (MQ=255)                        tgACATAGAATTGTTTCTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:1973783/100‑1 (MQ=255)                        tcttccgatcttGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:1410352/90‑1 (MQ=255)                        tcttccgatcttGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:4277142/90‑1 (MQ=255)                        tcttccgatcttGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:3937289/90‑1 (MQ=255)                        tcttccgatcttGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:1921013/90‑1 (MQ=255)                        tcttccgatcttGTTCTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:2919369/90‑1 (MQ=255)                        tcttccgatcttGTCTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt             <  1:3491555/90‑1 (MQ=255)                         gACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCttg            <  1:3533550/100‑2 (MQ=255)                         gACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTa            <  1:1890527/100‑1 (MQ=255)                                 cttGTTTCTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATcttt    <  1:1906206/99‑4 (MQ=255)                                   tGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGaa  <  1:897696/100‑1 (MQ=255)                                                           |                                                                            GATTTCGTTATTAAGTTTGAGTGTGACATAGAATTGTTTTTTATGATTACCATCTTTTTTAATATCAATGCGATCAATCACTGATAGATACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAA  >  NC_007793/45985‑46118

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