Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,402	T→A	V373V (GTA→GTT)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,402	0	T	A	100.0%	84.1 / NA	29	V373V (GTA→GTT)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (1/28);  total (1/28)

ATTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTACTACGCGTTGAATGACTTTATCACTTTTGACTATTTCAAGTATTTGCTTCATCACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATAC  >  NC_007793/46310‑46500                                                                                             |                                                                                                   gtttGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGtt                                                                                             <  1:1758007/99‑1 (MQ=255)  tttGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTg                                                                                            <  1:5337903/100‑1 (MQ=255)  tttGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTg                                                                                            <  1:2318032/100‑1 (MQ=255)   ttGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTgg                                                                                           <  1:4885458/100‑2 (MQ=255)   ttGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGa                                                                                           <  1:1280373/100‑1 (MQ=255)      tACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATg                                                                                        <  1:63052/100‑1 (MQ=255)        cTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGATACAACGCGTTGAATgca                                                                                      <  1:5226855/100‑3 (MQ=255)         ttCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATggtg                                                                                     <  1:3457821/100‑4 (MQ=255)          tcatcaTATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGaggt                                                                                    <  1:1265480/100‑4 (MQ=255)                   tgttgttgCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACtta                                                                           <  1:1492886/100‑2 (MQ=255)                        ttgCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACt                                                                      <  1:208586/100‑1 (MQ=255)                             tataAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTc                                                                 >  1:961419/1‑100 (MQ=255)                                               gaAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATtac                                               <  1:5461095/100‑1 (MQ=255)                                                              cAACTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATcac                                <  1:4176767/100‑1 (MQ=255)                                                                aCTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATcacca                              <  1:210142/100‑3 (MQ=255)                                                                 cTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATc                             <  1:3453916/100‑1 (MQ=255)                                                                 cTTTATTTTCTTGATTAACACGTGCTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATc                             <  1:3692448/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:1439757/99‑1 (MQ=255)                                                                                     cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg         <  1:5727299/99‑1 (MQ=255)                                                                                     cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg         <  1:5235614/99‑1 (MQ=255)                                                                                     cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg         <  1:509427/99‑1 (MQ=255)                                                                                     cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg         <  1:2233441/99‑1 (MQ=255)                                                                                     cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg         <  1:4616587/99‑1 (MQ=255)                                                                                      gtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGccag        <  1:1853455/100‑3 (MQ=255)                                                                                            ctacGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATAc  >  1:2572273/1‑100 (MQ=255)                                                                                             |                                                                                                   ATTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTACTACGCGTTGAATGACTTTATCACTTTTGACTATTTCAAGTATTTGCTTCATCACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGAACATAC  >  NC_007793/46310‑46500

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