Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	48,754	A→T	G46G (GGT→GGA)	SAUSA300_RS00205 ←	recombinase RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	48,754	0	A	T	100.0%	106.4 / NA	38	G46G (GGT→GGA)	SAUSA300_RS00205	recombinase RecA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (16/22);  total (16/22)

TATTAAGACGATAACAGCATAATATGTCACACTGTCCTTGTTGAATTAATTGTGTTATTTGTCGATACCCACTACGATTATCTTTGCGACCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAAGAGATTGTTGT  >  NC_007793/48666‑48850                                                                                         |                                                                                                 tATTAAGACGATAACAGCATAATATGTCACATTGTCCTTGTTGAATTAATTGTGTCATTTGTCGATACCCACTACGATTATCTTTGCGTCCTGAttgttt                                                                                       <  1:4812794/100‑1 (MQ=255) aacttctaCGATAACAGCATAATATGTCACATTGTCCTTGTTGAATTAATTGTGTCATTTGTCGATACCCACTACGATTATCTTTGCGTCCTGAttgttt                                                                                       >  1:5305/8‑100 (MQ=255)           aTAACAGCATAATATGTCACATTGTCCTTGTTGAATTAATTGTGTCATTTGTCGATACCCACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATaa                                                                             >  1:440428/1‑100 (MQ=255)                                     cTTGTTGAATTAATTGTGTCATTTGTCGATACCCACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGttt                                                   >  1:3920225/1‑100 (MQ=255)                                                      gtcaTTTGTCGATACCCACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGa                                  >  1:3992636/4‑100 (MQ=255)                                                                  tACCCACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCt                      >  1:782005/1‑100 (MQ=255)                                                                  tACCCACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCt                      <  1:1743894/100‑1 (MQ=255)                                                                   aCCCACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCtg                     >  1:5570020/1‑100 (MQ=255)                                                                   aCCCACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCtg                     >  1:4449381/1‑100 (MQ=255)                                                                    cccACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCtgt                    >  1:5111010/1‑100 (MQ=255)                                                                     ccACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCtgtg                   >  1:5495544/1‑100 (MQ=255)                                                                     ccACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCtgtg                   <  1:5576393/100‑1 (MQ=255)                                                                     ccACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCtgtg                   <  1:1886637/100‑1 (MQ=255)                                                                      cACTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGtgc                  >  1:753185/1‑100 (MQ=255)                                                                       aCTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGtgct                 <  1:4444778/100‑1 (MQ=255)                                                                        cTACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGtgcat                <  1:2444775/100‑3 (MQ=255)                                                                         tACGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGtgctgc               <  1:1735576/100‑1 (MQ=255)                                                                          tccgatctttttGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:3624434/91‑1 (MQ=255)                                                                          cttccgaTCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGtgctgct              <  1:5125559/94‑2 (MQ=255)                                                                          cttccgaTCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:2824291/94‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:3516568/100‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              >  1:3411903/1‑100 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:3029209/100‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:2979521/100‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:1269310/100‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:1845113/100‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              >  1:2585767/1‑100 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:609483/100‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:2461974/100‑1 (MQ=255)                                                                          aCGATTATCTTTGCGTCCTGATTGTTTGTAGCGATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCa              <  1:5149365/100‑1 (MQ=255)                                                                            gATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAag            <  1:2634179/100‑1 (MQ=255)                                                                            gATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAag            <  1:2086765/100‑1 (MQ=255)                                                                            gATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAag            <  1:5175191/100‑1 (MQ=255)                                                                            gATTATCTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAag            <  1:2708320/100‑1 (MQ=255)                                                                                  cTTTGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGC‑AGAGAtcgg     >  1:5530097/1‑97 (MQ=255)                                                                                    ttGCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAAGAGAttgtt    >  1:867685/1‑100 (MQ=255)                                                                                    ccgCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTTTTTTTGTGCTGCAagagctcttc    >  1:630728/3‑94 (MQ=255)                                                                                      gCGTCCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAagatcggaagag  >  1:7935/1‑91 (MQ=255)                                                                                         |                                                                                                 TATTAAGACGATAACAGCATAATATGTCACACTGTCCTTGTTGAATTAATTGTGTTATTTGTCGATACCCACTACGATTATCTTTGCGACCTGATTGTTTGTCGCTATAAAAGTTGATATGTTGAATATGATGTTTTTCGGCTATTGCTTCGATAGCTTGTTTCTGTGCTGCAAGAGATTGTTGT  >  NC_007793/48666‑48850

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