Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,398	T→C	T375A (ACA→GCA)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,398	0	T	C	95.2%	51.3 / ‑5.5	21	T375A (ACA→GCA)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (1/19);  minor base A (0/1);  total (1/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTATTTTGATTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTACTACGCGTTGAATGACTTTATCACTTTTGACTATTTCAAGTATTTGCTTCATCACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGA  >  NC_007793/46302‑46494                                                                                                 |                                                                                                 ttagttTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCtac                                                                                               <  1:963282/100‑1 (MQ=255) ttagttTGGTTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTTATTTTCTTGATTAACACGTGCtac                                                                                               <  1:215338/100‑1 (MQ=255)         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:4616587/99‑1 (MQ=255)                                                                                             cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg   <  1:509427/99‑1 (MQ=255)                                                                                             cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg   <  1:2233441/99‑1 (MQ=255)                                                                                             cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg   <  1:5235614/99‑1 (MQ=255)                                                                                             cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg   <  1:1439757/99‑1 (MQ=255)                                                                                             cgtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCg   <  1:5727299/99‑1 (MQ=255)                                                                                              gtgcTACAACGCGTTGAATGACTTTATCACTTTTGACTATTTCGAGTATTTGCTTCATTACATAATCTTCAATCACATCAGCTCTTACACTGTTTGccag  <  1:1853455/100‑3 (MQ=255)                                                                                                 |                                                                                                 TTATTTTGATTTGTACTTCATCATATTGTTGTTGCTTATAAGCAATATCGTGATGAAGTGCAGCGCCATCAACTTGATTTTCTTGATTGACGTGTGTTACTACGCGTTGAATGACTTTATCACTTTTGACTATTTCAAGTATTTGCTTCATCACATAATCTTCAATCACATCAGCTCTTACACTGTTTGCCGA  >  NC_007793/46302‑46494

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