Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,438	A→G	L361L (CTT→CTC)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,438	0	A	G	100.0%	71.8 / NA	28	L361L (CTT→CTC)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/22);  total (6/22)

GAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTACTACGCGTTGAATGACTTTATCACTTTTGACTATTTCAAGTATTTGCTTCATCACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCTTGTTCCGAAAATTACTACATGAA  >  NC_007793/46356‑46534                                                                                   |                                                                                                 gaAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATtac                                                                                 <  1:5461095/100‑1 (MQ=255)                cAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATcac                                                                  <  1:4176767/100‑1 (MQ=255)                  aCTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATcacca                                                                <  1:210142/100‑3 (MQ=255)                   cTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATc                                                               <  1:3692448/100‑1 (MQ=255)                   cTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATc                                                               <  1:3453916/100‑1 (MQ=255)                   cTTTATTCTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATc                                                               <  1:5059780/100‑1 (MQ=255)                               gattgaCACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTacaca                                                   <  1:5671130/100‑2 (MQ=255)                                   aacacgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGttt                                               <  1:4115309/95‑1 (MQ=255)                                    acacgtgcTACAACGCGTCGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGttgc                                              <  1:302236/96‑3 (MQ=255)                                     cacgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGc                                             <  1:1601367/97‑1 (MQ=255)                                       cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg                                           <  1:2233441/99‑1 (MQ=255)                                       cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg                                           <  1:1439757/99‑1 (MQ=255)                                       cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg                                           <  1:4616587/99‑1 (MQ=255)                                       cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg                                           <  1:5727299/99‑1 (MQ=255)                                       cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg                                           <  1:509427/99‑1 (MQ=255)                                       cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg                                           <  1:5235614/99‑1 (MQ=255)                                        gtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGccag                                          <  1:1853455/100‑3 (MQ=255)                                             acaacGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATa                                     >  1:2460344/4‑100 (MQ=255)                                              ctacGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATAc                                    >  1:2572273/1‑100 (MQ=255)                                               aacGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACt                                   >  1:4669932/2‑100 (MQ=255)                                                acGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACtt                                  <  1:4861763/100‑1 (MQ=255)                                                   cgTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTgg                               <  1:4742209/100‑1 (MQ=255)                                                            aCTTTATCTCTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCTTgg                      <  1:2834644/100‑2 (MQ=255)                                                              atctacACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCtgcatg                    <  1:4499801/95‑6 (MQ=255)                                                                         tGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCTTGTTCCGAAAATtact         >  1:4322418/1‑100 (MQ=255)                                                                         tGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCTTGTTCCGAAAATtact         >  1:5384136/1‑100 (MQ=255)                                                                              aTTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCTTGTTCCGAAAATTACTACATg    >  1:1875233/1‑100 (MQ=255)                                                                                ttCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCTTGTTCCGAAAATTACTACagggg  <  1:727912/100‑5 (MQ=255)                                                                                   |                                                                                                 GAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTACTACGCGTTGAATGACTTTATCACTTTTGACTATTTCAAGTATTTGCTTCATCACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATACTTTGGAACCCTTGTTCCGAAAATTACTACATGAA  >  NC_007793/46356‑46534

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