Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,753,451	T→G	S148A (TCT→GCT)	SAUSA300_RS14145 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,753,451	0	T	G	100.0%	56.9 / NA	23	S148A (TCT→GCT)	SAUSA300_RS14145	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (23/0);  total (23/0)

CAATTCATTACTTGAAAGCAGTGATTTAACTGAAGCTGAAACACAAGAGGCAC‑AAAAAATTAAATCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAAATTTTCAAATAT  >  NC_007793/2753387‑2753534                                                                  |                                                                                    cAATTCATTACTTGAAAGCAGTGATTTAACTGAAGCTGAAACACAAGAGGCACAAAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATg                                                   >  1:3643575/1‑100 (MQ=255)             tGAAAGCAGTGATTTAACTGAAGCTGAAACACAAGAGGCAC‑AAAAAATTATAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTa                                      >  1:1780064/1‑100 (MQ=255)                taGCAGTGATTTAACTGAAGCTGAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAAcc                                   >  1:116570/2‑100 (MQ=255)                    tgTGATTTAACTGAAGCTGAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAg                               >  1:3968989/2‑100 (MQ=255)                            aaCTGAAGCTGAAACACAAGAGGCAC‑AAAGAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGaa                       >  1:1106579/1‑100 (MQ=255)                                     tGAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGaa              >  1:3747830/1‑100 (MQ=255)                                     accgacacAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCggacag              >  1:5534247/5‑96 (MQ=255)                                      gAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTGCGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGaaa             >  1:6391744/1‑100 (MQ=255)                                      gAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATTGTTAGAACAGTTAACCGAAGCGATTGAAAGCGgagatc             >  1:4109173/1‑98 (MQ=255)                                      gAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGgagatc             >  1:6004074/1‑98 (MQ=255)                                      gAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGtagat             >  1:3792036/1‑96 (MQ=255)                                      gAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGaaa             >  1:5569444/1‑100 (MQ=255)                                      gAAACACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGaaa             >  1:3566892/1‑100 (MQ=255)                                        cccacaAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAAAtt           >  1:5917547/3‑100 (MQ=255)                                        cacacAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAGAtc           >  1:2692169/2‑99 (MQ=255)                                        aaCACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGgagatcgg           >  1:3648637/1‑96 (MQ=255)                                        aaCACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAGAtc           >  1:151057/1‑99 (MQ=255)                                        aaCACAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAAAtt           >  1:67176/1‑100 (MQ=255)                                         acacAAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGtagatcgg          >  1:319236/1‑93 (MQ=255)                                            caAGAGGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAGAtcggaa       >  1:3084585/1‑95 (MQ=255)                                               cagGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAAATTTTCAAat    >  1:344078/2‑100 (MQ=255)                                                agGCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAAATTTTCAAata   >  1:1726217/1‑100 (MQ=255)                                                 ggCAC‑AAAAAATTAAAGCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAAATTTTCAAatat  >  1:3356894/1‑100 (MQ=255)                                                                  |                                                                                    CAATTCATTACTTGAAAGCAGTGATTTAACTGAAGCTGAAACACAAGAGGCAC‑AAAAAATTAAATCTCAGCTTGATCATACACATGTTTATTACGAATGGTTAGAACAGTTAACCGAAGCGATTGAAAGCGGTGAAATTTTCAAATAT  >  NC_007793/2753387‑2753534

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