Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,099	T→C	I474M (ATA→ATG)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,099	0	T	C	100.0%	79.0 / NA	29	I474M (ATA→ATG)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/27);  total (2/27)

GTGACATAGAATTGTTTTTTATGATTACCATCTTTTTTAATATCAATGCGATCAATCACTGATAGATACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAAAGATGTGTTGTAATAGTTTACTGATTTCTTTGGCATCAAATAAAGTCTTATCTTC  >  NC_007793/46007‑46173                                                                                             |                                                                           gtgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCt                                                                     >  1:4739178/1‑100 (MQ=255)  tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGGGATTTTTCCATGTGCtt                                                                    <  1:695619/100‑1 (MQ=255)  tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt                                                                    <  1:1522286/100‑1 (MQ=255)  tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt                                                                    <  1:5743295/100‑1 (MQ=255)  tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt                                                                    <  1:466737/100‑1 (MQ=255)  tgACATAGAATTGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt                                                                    <  1:239737/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:1921013/90‑1 (MQ=255)  tcttccgatcttGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt                                                                    <  1:4277142/90‑1 (MQ=255)  tcttccgatcttGTTTTTTATGATTACCATCTTTTCTAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCtt                                                                    <  1:3937289/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)                                    tctAATATTAATACGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTg                                  <  1:5781827/98‑1 (MQ=255)                                              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:2524489/100‑1 (MQ=255)                                                 cGATCAATCACTGATAGGTACAATGCTTTGAGTCGAGATTTTTCCATGTGCTTAATATCTTGAAAGATGTGTTGTAATAGTTTACTGATTTGTTTGGCAt                     <  1:2230666/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)                                                                                             |                                                                           GTGACATAGAATTGTTTTTTATGATTACCATCTTTTTTAATATCAATGCGATCAATCACTGATAGATACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAAAGATGTGTTGTAATAGTTTACTGATTTCTTTGGCATCAAATAAAGTCTTATCTTC  >  NC_007793/46007‑46173

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