Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,684	C→T	K279K (AAG→AAA)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,684	0	C	T	100.0%	90.1 / NA	33	K279K (AAG→AAA)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/25);  total (8/25)

CTGCCATAGGTGCGCCACATTGGGGACAGTGAATAATGCCTGTAAGCAGATTCGTTCCTTTGCCATGGACTTGGGGTTTTTGACTGACTTGTTTCTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTATCATTCAGCCC  >  NC_007793/46590‑46770                                                                                               |                                                                                       ctgcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACg                                                                                   >  1:3294330/1‑100 (MQ=255)  tgcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGc                                                                                  <  1:156571/100‑1 (MQ=255)   gcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCa                                                                                 >  1:446024/1‑100 (MQ=255)   gcCATAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCa                                                                                 >  1:5728413/1‑100 (MQ=255)    ccgtAGGTGCGCCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCAt                                                                                >  1:2791684/4‑100 (MQ=255)            gcgcCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTTCCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACta                                                                        <  1:404656/100‑2 (MQ=255)            gcgcCACATTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACtt                                                                        >  1:2553641/1‑100 (MQ=255)                 acaTTGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATcc                                                                   <  1:3217511/100‑1 (MQ=255)                    ttGGGGACAATTGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGAGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCAt                                                                <  1:3718240/100‑1 (MQ=255)                    ttGGGGACAATGGATAATGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCAt                                                                <  1:3788460/100‑1 (MQ=255)                           caatggataatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCtt                                                         <  1:3758372/94‑1 (MQ=255)                            aatggataatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTg                                                        <  1:1912494/95‑1 (MQ=255)                                  caatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATtaat                                                  >  1:4542211/2‑100 (MQ=255)                                  agatGCCTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATtaat                                                  >  1:829127/3‑100 (MQ=255)                                       ccTGTAAGCAGATTCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATgg                                             <  1:3902717/100‑1 (MQ=255)                                                   ttCGTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTAc                                 <  1:1781260/100‑1 (MQ=255)                                                      gTTCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCtt                              <  1:3947526/100‑1 (MQ=255)                                                       ttCCTTTACCATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTAccctg                             <  1:1625192/100‑4 (MQ=255)                                                             taccATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTAt                       <  1:2171422/98‑1 (MQ=255)                                                              accATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTAtg                      <  1:359361/99‑2 (MQ=255)                                                               ccATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTacta                     <  1:5071557/100‑4 (MQ=255)                                                               ccATGAACTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCa                     <  1:1127563/100‑1 (MQ=255)                                                                     atctGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGta               <  1:2247327/97‑2 (MQ=255)                                                                      cTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGttt              <  1:3219722/100‑1 (MQ=255)                                                                      cTTGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGttt              <  1:3820446/100‑1 (MQ=255)                                                                       ttGGGGTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTa             <  1:4795915/100‑1 (MQ=255)                                                                         ggggTTTTTGACTGACTTGTTTTTTGCGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGtttcac           <  1:4463095/100‑4 (MQ=255)                                                                         ggggTTTTTGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAACTATCACTGGTTTATc           <  1:197786/100‑1 (MQ=255)                                                                             tttttGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTATCATTc       <  1:1230287/100‑1 (MQ=255)                                                                             tttttGACTGACTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTATCATTc       <  1:341768/100‑1 (MQ=255)                                                                                 tgactgacTTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTATCATTCAGct   <  1:1958438/100‑2 (MQ=255)                                                                                  gactgactTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTATCATTCAggca  <  1:862105/100‑4 (MQ=255)                                                                                  gactgactTGTTTTTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTATCATTCAGccc  >  1:130853/1‑100 (MQ=255)                                                                                               |                                                                                       CTGCCATAGGTGCGCCACATTGGGGACAGTGAATAATGCCTGTAAGCAGATTCGTTCCTTTGCCATGGACTTGGGGTTTTTGACTGACTTGTTTCTTACGCATTTGTACTTTATCCCATAAATCTTGATTAATAATGGGGGAATGCTTACCTTCAGCTATCACTGGTTTATCATTCAGCCC  >  NC_007793/46590‑46770

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