Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,319,806	C→T	Q52* (CAA→TAA)	SAUSA300_RS06470 →	GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,319,806	0	C	T	100.0%	85.4 / NA	30	Q52* (CAA→TAA)	SAUSA300_RS06470	GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/3);  total (27/3)

TACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAA  >  NC_007793/1319714‑1319894                                                                                             |                                                                                         tACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTAc                                                                                   >  1:3709461/1‑100 (MQ=255)  aCATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACt                                                                                  >  1:5007825/1‑100 (MQ=255)  aCATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACt                                                                                  >  1:3590846/1‑100 (MQ=255)    aTGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCa                                                                                >  1:1266139/1‑100 (MQ=255)         cAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACa                                                                           >  1:2547993/1‑100 (MQ=255)          aaGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACag                                                                          <  1:5398160/100‑1 (MQ=255)                   aGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGACCgtgt                                                                 <  1:5719727/100‑1 (MQ=255)                                     atACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGt                                               >  1:5224902/2‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTCTGTTg                                             >  1:5434870/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:2365379/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:471019/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:5275750/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:1036093/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:2803124/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:2430025/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:2354714/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:2171694/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:122483/1‑100 (MQ=255)                                       aCAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTg                                             >  1:737856/1‑100 (MQ=255)                                                            ttATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGaaa                        >  1:1697839/1‑100 (MQ=255)                                                            ttATCAGAGACTTGCCAACTTGAAGGGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGaaa                        >  1:1168256/1‑100 (MQ=255)                                                             tATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAAt                       >  1:3119973/1‑100 (MQ=255)                                                             tATCAGAGACTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAAt                       >  1:1926983/1‑100 (MQ=255)                                                               tgtgagaCTTGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTc                     >  1:1299914/4‑100 (MQ=255)                                                                       ttGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCAttt             >  1:4615461/1‑100 (MQ=255)                                                                        tGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTa            >  1:484195/1‑100 (MQ=255)                                                                        tGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTa            >  1:4512237/1‑100 (MQ=255)                                                                        tGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTa            >  1:1396007/1‑100 (MQ=255)                                                                        tGCCAACTTGAAGTGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATCATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTa            >  1:1072841/1‑100 (MQ=255)                                                                                  aaGCGTTGGGTTAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCaa  <  1:5224448/100‑1 (MQ=255)                                                                                             |                                                                                         TACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCAA  >  NC_007793/1319714‑1319894

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