Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,138,166	G→A	intergenic (‑187/+167)	SAUSA300_RS10880 ← / ← SAUSA300_RS10885	membrane protein/phage terminase small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,138,166	0	G	A	100.0%	50.6 / NA	19	intergenic (‑187/+167)	SAUSA300_RS10880/SAUSA300_RS10885	membrane protein/phage terminase small subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (2/17);  total (2/17)

CCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑GAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTATTAGGATATTGCGTTTACAACTCTTTTA  >  NC_007793/2138103‑2138255                                                                 |                                                                                          cccTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGCAATGTTGATGCAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTccc                                                       >  1:3882722/1‑100 (MQ=255) cccTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTcat                                                       <  1:6011046/100‑3 (MQ=255)  ccTTGTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCt                                                      <  1:1617395/100‑1 (MQ=255)   cTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTa                                                     <  1:1426381/100‑1 (MQ=255)       ttCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAatc                                                 <  1:6000914/100‑1 (MQ=255)         cTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAatcat                                               <  1:462812/100‑1 (MQ=255)           tAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAatcagcc                                             <  1:1420680/100‑4 (MQ=255)           tAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAAtcatca                                             <  1:876428/100‑1 (MQ=255)           tAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTA‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAAtcatca                                             <  1:2599433/100‑1 (MQ=255)            aaaTCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCaa                                            <  1:4233879/100‑1 (MQ=255)            aaaTCATTTTATTATATATTACAATAATCAATCTGAAATGTCGATGTAATTTAAAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAAtcatca                                             <  1:737379/100‑1 (MQ=255)             aaTCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAat                                           <  1:3459109/100‑1 (MQ=255)                         atatatTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCGCCCTAAATCATCAATATGAAAATCAGt                               <  1:1871691/100‑1 (MQ=255)                            tatTACAATAATCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTAtt                            <  1:4607775/100‑1 (MQ=255)                                   ataatCAATCTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTATTAGGatat                     <  1:1199983/100‑1 (MQ=255)                                           tcTGAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTATTAGGATATTGCGTTTa             <  1:4085323/100‑1 (MQ=255)                                              gAAATGTTGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTATTAGGATATTGCGTTTACaa          >  1:6198191/1‑100 (MQ=255)                                                     tGATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTATTAGGATATTGCGTTTACAACTCtttt   <  1:4571961/100‑1 (MQ=255)                                                      gATGTAATTT‑AAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTATTAGGATATTGCGTTTACAACTCTTTTa  <  1:3023527/100‑1 (MQ=255)                                                                 |                                                                                          CCCTTTTTCTTAAATCATTTTATTATATATTACAATAATCAATCTGAAATGTTGATGTAATTT‑GAAAAAAATATCATACTTTTACTCCTGAAAACCTCCCTAAATCATCAATATGAAAATCAGTATTAGGATATTGCGTTTACAACTCTTTTA  >  NC_007793/2138103‑2138255

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