Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,146	C→G	E459Q (GAA→CAA)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,146	0	C	G	100.0%	60.3 / NA	24	E459Q (GAA→CAA)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (2/22);  total (2/22)

ATGCGATCAATCACTGATAGATACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAAAGATGTGTTGTAATAGTTTACTGATTTCTTTGGCATCAAATAAAGTCTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCAT  >  NC_007793/46052‑46232                                                                                               |                                                                                       atacGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTgg                                                                                   <  1:2085558/97‑1 (MQ=255) atacGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTCAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTgg                                                                                   <  1:5643970/97‑1 (MQ=255)  tccGATCTATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGc                                                                                  <  1:279526/98‑1 (MQ=255)  tacGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTGTGGc                                                                                  <  1:4798127/98‑1 (MQ=255)    cGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCAt                                                                                <  1:2230666/100‑1 (MQ=255)    cGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCAt                                                                                <  1:2524489/100‑1 (MQ=255)     gATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCAta                                                                               <  1:2877909/100‑2 (MQ=255)       tctatcaCTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCaa                                                                             <  1:1816485/97‑1 (MQ=255)                       aCAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTc                                                             >  1:4754349/1‑100 (MQ=255)                                               tcttgtgCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAgttggtt                                     <  1:3920847/97‑1 (MQ=255)                                                caTGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAgttggttg                                    <  1:3097049/99‑1 (MQ=255)                                                           gatcttgcacgATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGt                         <  1:3760819/90‑1 (MQ=255)                                                                tGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGt                    >  1:3281294/1‑100 (MQ=255)                                                                 gAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTc                   <  1:2553624/100‑1 (MQ=255)                                                                   aaGATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCAt                 <  1:5068852/100‑1 (MQ=255)                                                                     gATGTGTTGTAATAGTTTACTGATTTGTTTGGTATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTg               <  1:1412076/100‑1 (MQ=255)                                                                     gATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCAtat               <  1:5307493/100‑3 (MQ=255)                                                                     gATGTGTTGTAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTg               <  1:940421/100‑1 (MQ=255)                                                                         tgttgtAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGtt           <  1:5467727/100‑1 (MQ=255)                                                                          gttgtAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGtgc          <  1:2731335/100‑3 (MQ=255)                                                                             gtAATAGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCtg       <  1:2053238/100‑2 (MQ=255)                                                                             gtAATAGTTTACTGATTAGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCta       <  1:3716569/100‑2 (MQ=255)                                                                                  aGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCat  <  1:4541121/100‑1 (MQ=255)                                                                                  aGTTTACTGATTTGTTTGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCat  <  1:5336898/100‑1 (MQ=255)                                                                                               |                                                                                       ATGCGATCAATCACTGATAGATACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAAAGATGTGTTGTAATAGTTTACTGATTTCTTTGGCATCAAATAAAGTCTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCAT  >  NC_007793/46052‑46232

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