Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	716,517	G→A	intergenic (+23/+142)	SAUSA300_RS03435 → / ← SAUSA300_RS03440	lipase LipA/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	716,517	0	G	A	100.0%	65.0 / NA	26	intergenic (+23/+142)	SAUSA300_RS03435/SAUSA300_RS03440	lipase LipA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/1);  total (25/1)

CAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTGCTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAAT  >  NC_007793/716429‑716600                                                                                         |                                                                                    cAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACa                                                                          >  1:30297/1‑100 (MQ=255)  aaaCTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACAt                                                                         >  1:391823/1‑100 (MQ=255)    ccTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACAttt                                                                       >  1:2982726/2‑100 (MQ=255)            gaAACCACAAATAGAAAATCCGAGTAATGAATTACCTTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGAtttt                                                               >  1:1277412/1‑100 (MQ=255)            gaAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATCTAAGAtttt                                                               >  1:2273580/1‑100 (MQ=255)             taaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTc                                                              >  1:1290899/2‑100 (MQ=255)             gaaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTCTTACTTCAGCAACATTTAAGATTTTc                                                              >  1:5248035/2‑100 (MQ=255)             caaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTc                                                              >  1:1848531/2‑100 (MQ=255)             aaaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTGATTTATTACTTAAGCAACATTTAAGCTTTTc                                                              >  1:4155860/1‑100 (MQ=255)             aaaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTc                                                              >  1:3309975/1‑100 (MQ=255)              taCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAGCATTTAAGATTTTCa                                                             >  1:2103713/2‑100 (MQ=255)              caCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCa                                                             >  1:6614638/2‑100 (MQ=255)              aaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCa                                                             >  1:2629083/1‑100 (MQ=255)              aaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCa                                                             >  1:5883382/1‑100 (MQ=255)               cccACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaa                                                            >  1:3747954/2‑100 (MQ=255)               aCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaa                                                            >  1:3085948/1‑100 (MQ=255)                gcccaAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaaa                                                           >  1:5852109/4‑100 (MQ=255)                ccACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaaa                                                           >  1:6334452/1‑100 (MQ=255)                ccACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaaa                                                           >  1:1795925/1‑100 (MQ=255)                  ccaAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAAtt                                                         >  1:3756171/2‑100 (MQ=255)                              gccGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACgaagat                                             >  1:974401/2‑99 (MQ=255)                                          aTTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACAtgg                                 <  1:1180623/100‑2 (MQ=255)                                                   caaTCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTaga                        >  1:1879981/2‑98 (MQ=255)                                                   aaaTCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGAGGTGCTAAt                        >  1:6643256/1‑100 (MQ=255)                                                                 aaTGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTATAagagatcgg          >  1:5908595/1‑95 (MQ=255)                                                                         aCAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATataat  >  1:3721520/1‑100 (MQ=255)                                                                                         |                                                                                    CAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTGCTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAAT  >  NC_007793/716429‑716600

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