Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,164	T→A	T453S (ACT→TCT)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,164	0	T	A	100.0%	92.2 / NA	33	T453S (ACT→TCT)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (3/30);  total (3/30)

ACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAAAGATGTGTTGTAATAGTTTACTGATTTCTTTGGCATCAAATAAAGTCTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTG  >  NC_007793/46074‑46259                                                                                           |                                                                                                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)                                                                         gtttGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACt                >  1:3956905/2‑100 (MQ=255)                                                                          tttGGCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACtt               <  1:3246502/100‑1 (MQ=255)                                                                              gCATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTc           <  1:5057914/100‑1 (MQ=255)                                                                               cATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTca          <  1:2153988/100‑2 (MQ=255)                                                                               cATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCt          <  1:4516987/100‑1 (MQ=255)                                                                               cATCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCt          <  1:5328213/100‑1 (MQ=255)                                                                                aTCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTg         <  1:2146962/100‑1 (MQ=255)                                                                                aTCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTg         <  1:4796756/100‑1 (MQ=255)                                                                                 tCAAATAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGa        <  1:2605858/100‑1 (MQ=255)                                                                                     aTAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTAc    <  1:2603118/100‑1 (MQ=255)                                                                                     aTAAAGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTAc    <  1:1176234/100‑1 (MQ=255)                                                                                       aaaGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTg  <  1:4756219/100‑1 (MQ=255)                                                                                       aaaGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTg  <  1:2129927/100‑1 (MQ=255)                                                                                       aaaGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTg  <  1:2571779/100‑1 (MQ=255)                                                                                       aaaGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTg  <  1:1267407/100‑1 (MQ=255)                                                                                       aaaGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTg  <  1:3932599/100‑1 (MQ=255)                                                                                       aaaGACTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTg  <  1:71478/100‑1 (MQ=255)                                                                                           |                                                                                                ACAATGCTTTGAGTCGAGATTTTTCTATGTGCTTAATATCATGAAAGATGTGTTGTAATAGTTTACTGATTTCTTTGGCATCAAATAAAGTCTTATCTTCATTTTGTTGATTTTTGAGTTGGTTGATTTGATTCGTAATGTCATTGAGTTGCTTTTCATATTTTTGAATACTTGGTCTGATTACTG  >  NC_007793/46074‑46259

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: ‑