Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,467	G→A	S120S (AGC→AGT)	SAUSA300_RS07730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,467	0	G	A	100.0%	47.4 / NA	18	S120S (AGC→AGT)	SAUSA300_RS07730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/9);  total (9/9)

GAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTGGTCAGGATAAAATTCTTTGAATATTATCTGTGGTGTGCTAACAACTGCGTCTTGAGCAAATTCTTTATCTTTAAAATTAAACATGTTGTGAAATTCACTATTTTTAAAACTTGATTCAATCGCTTCTTT  >  NC_007793/1579371‑1579559                                                                                                 |                                                                                             gAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTa                                                                                           <  1:3701064/100‑1 (MQ=255)           aTAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCgc                                                                                 <  1:190237/100‑2 (MQ=255)            tAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCgcc                                                                                <  1:1149751/100‑3 (MQ=255)              gtatatAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCtg                                                                              <  1:3847484/98‑2 (MQ=255)               tatatAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTg                                                                             <  1:1733830/99‑1 (MQ=255)                atatAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGa                                                                            <  1:2293079/100‑1 (MQ=255)                  ataAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCttagat                                                                          <  1:2678909/100‑5 (MQ=255)                        tgtaaGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATtc                                                                    >  1:1861665/4‑100 (MQ=255)                                aTGTCTATCGCTATCACGATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATtctttatctt                                                            <  1:5723018/100‑1 (MQ=255)                                                  aTTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACAttttgt                                          >  1:345257/1‑100 (MQ=255)                                                    tGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGa                                        >  1:2331530/1‑100 (MQ=255)                                                    tGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGa                                        >  1:5145958/1‑100 (MQ=255)                                                     tccAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGaa                                       >  1:1428089/2‑100 (MQ=255)                                                     gCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGaa                                       >  1:3859205/1‑100 (MQ=255)                                                     gCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGaa                                       >  1:5329101/1‑100 (MQ=255)                                                                                     tatCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGAAATTCTGTATCTTTAAAACTTGATTCAATCGct       <  1:2706254/100‑1 (MQ=255)                                                                                         tgtgGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGAAATTCTGTATCTTTAAAACTTGATTCAATCGcttctt   >  1:2658032/1‑100 (MQ=255)                                                                                          gtggtgTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGAAATTCTGTATCTTTAAAACTTGATTCAATCGCTTCttt  >  1:919418/1‑100 (MQ=255)                                                                                                 |                                                                                             GAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTGGTCAGGATAAAATTCTTTGAATATTATCTGTGGTGTGCTAACAACTGCGTCTTGAGCAAATTCTTTATCTTTAAAATTAAACATGTTGTGAAATTCACTATTTTTAAAACTTGATTCAATCGCTTCTTT  >  NC_007793/1579371‑1579559

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