Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,503,835	T→C	D111D (GAT→GAC)	SAUSA300_RS12860 →	general stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,503,835	0	T	C	100.0%	73.0 / NA	30	D111D (GAT→GAC)	SAUSA300_RS12860	general stress protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/2);  total (28/2)

GAGATGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGATCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGATTTATAACACAA  >  NC_007793/2503740‑2503930                                                                                                |                                                                                                gagaTGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCag                                                                                             >  1:1111475/1‑100 (MQ=255)   cctGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCagag                                                                                           >  1:5280161/3‑100 (MQ=255)     tGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTCTTCAGTTCAAAAGAAGACCCAGAGtt                                                                                         <  1:3845169/100‑1 (MQ=255)         cgCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCagatcggaag                                                                                     >  1:2455744/2‑93 (MQ=255)                  cGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAaga                                                                            >  1:4008523/1‑100 (MQ=255)                       ctgctACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGagtagt                                                                       >  1:944972/5‑100 (MQ=255)                        cGTTACAGAACAGGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGtt                                                                      >  1:2461443/1‑100 (MQ=255)                         cctACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTc                                                                     >  1:4224965/3‑100 (MQ=255)                          ttACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTcc                                                                    >  1:522260/1‑100 (MQ=255)                            aCAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGc                                                                  >  1:1768299/1‑100 (MQ=255)                            aCAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGc                                                                  >  1:1234637/1‑100 (MQ=255)                                  cAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCa                                                            >  1:1532667/1‑100 (MQ=255)                                              ttttgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatga                                                >  1:2257015/3‑100 (MQ=255)                                               accggCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTTCGCAATCCATTAAGCTaatgaa                                               >  1:1024947/4‑100 (MQ=255)                                                ttgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaat                                              >  1:4614421/1‑100 (MQ=255)                                                 tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGCCCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaatg                                             >  1:1416510/1‑100 (MQ=255)                                                 tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaatg                                             >  1:3228147/1‑100 (MQ=255)                                                 tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaatg                                             >  1:4479937/1‑100 (MQ=255)                                                 tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGATCCGCAATCCATTAAGCTaatgaatg                                             >  1:5543266/1‑100 (MQ=255)                                                  ggCTATGGGAGACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAatgaatga                                            >  1:3525052/1‑100 (MQ=255)                                                  ggCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAatgaatga                                            >  1:900739/1‑100 (MQ=255)                                                        gggAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATc                                      >  1:502186/1‑100 (MQ=255)                                                          gAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCAt                                    >  1:6452621/1‑100 (MQ=255)                                                              cacaAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGa                                >  1:1425851/1‑100 (MQ=255)                                                                 caGACAAAAGCTTTTTCAGTTCAAAAGCAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATac                             >  1:112217/2‑100 (MQ=255)                                                                      taaaGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTa                        >  1:2904233/2‑100 (MQ=255)                                                                        aaGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATc                      <  1:5360679/100‑1 (MQ=255)                                                                                          tgaagaCCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGATTTATAacac    >  1:1103104/2‑100 (MQ=255)                                                                                           caagaCCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGATTTATAAcaca   >  1:2442124/2‑100 (MQ=255)                                                                                            cagaCCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGATTTATAACACaa  >  1:4588873/2‑100 (MQ=255)                                                                                                |                                                                                                GAGATGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGATCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGATTTATAACACAA  >  NC_007793/2503740‑2503930

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