Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	47,290	G→A	D77D (GAC→GAT)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	47,290	0	G	A	100.0%	54.1 / NA	20	D77D (GAC→GAT)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/5);  total (15/5)

TTTGTTTGTGTAAATCTTCGACAATTTTGAGAAGATCAGATGTATTACGAGCTAATCGGTTTGTTTTGTATACCATAACACAGTCGATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGGTCAGCGTATATATC  >  NC_007793/47208‑47382                                                                                   |                                                                                             tttgtttaTGCAAATCTTCGACGATTTTCAGAAGATCCGATGTATTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCttt                                                                             >  1:1696655/1‑100 (MQ=255)        atgcaaaTCTTCGACGATTTTCAGAAGATCCGATGTATTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCAtac                                                                      <  1:118557/96‑3 (MQ=255)                                     ccgATGTATTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGtt                                         >  1:2573188/3‑100 (MQ=255)                                      ccaTGTATTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTTCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTc                                        >  1:234237/3‑100 (MQ=255)                                        ttGTATTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCAt                                      >  1:3759972/2‑100 (MQ=255)                                        aTGTATTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCAt                                      <  1:4991078/100‑1 (MQ=255)                                        aTGTATTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCAt                                      <  1:2785869/100‑1 (MQ=255)                                            aTTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGAt                                  >  1:4953004/1‑100 (MQ=255)                                            aTTACGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGAt                                  >  1:4818164/1‑100 (MQ=255)                                             ccaCGAGCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGAtt                                 >  1:1142688/3‑100 (MQ=255)                                                   gCTAATCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTg                           >  1:5069427/1‑100 (MQ=255)                                                      aaTCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGaaa                        >  1:2311607/1‑100 (MQ=255)                                                      aaTCGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGaaa                        >  1:4006212/1‑100 (MQ=255)                                                         cGGTTTGTTTTGTAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATAc                     >  1:530363/1‑100 (MQ=255)                                                             ttgttttgtAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACg                 >  1:1444927/1‑100 (MQ=255)                                                              tgttttgtAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGCTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCAcca                <  1:1056822/100‑3 (MQ=255)                                                                   tgtAGATCATAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTTCCTGAAATACCACGATCAGc           <  1:5582081/100‑1 (MQ=255)                                                                            tAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATc  >  1:869067/1‑100 (MQ=255)                                                                            tAACACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATc  >  1:1992261/1‑100 (MQ=255)                                                                            aaaCACAATCTATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGATCAGCGTATATATc  >  1:3732217/2‑100 (MQ=255)                                                                                   |                                                                                             TTTGTTTGTGTAAATCTTCGACAATTTTGAGAAGATCAGATGTATTACGAGCTAATCGGTTTGTTTTGTATACCATAACACAGTCGATATTGCCTTCTTTTGCATCTTTCAACATACGTTGGAGCTCAGGTCTGTTCATAGATTTACCTGAAATACCACGGTCAGCGTATATATC  >  NC_007793/47208‑47382

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