Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	453,610	T→G	N81K (AAT→AAG)	SAUSA300_RS02150 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	453,610	0	T	G	100.0%	75.0 / NA	28	N81K (AAT→AAG)	SAUSA300_RS02150	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (26/2);  total (26/2)

CTACTCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAATCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAATATAAAGAAGGCATTGAAGGCAAAGATGTCTTTGTGGTAAAAGAATTAA  >  NC_007793/453514‑453695                                                                                                 |                                                                                      atactCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCAc                                                                                    >  1:2740212/2‑100 (MQ=255)  tactCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCAcc                                                                                   >  1:1687417/1‑100 (MQ=255)  tactCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGGAAGGATGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCAcc                                                                                   >  1:4561335/1‑100 (MQ=255)   actCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTCAGGTTGATAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCa                                                                                  >  1:1047603/1‑100 (MQ=255)   actCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCATGAAAATCAAAAGCACCa                                                                                  >  1:1870887/1‑100 (MQ=255)   actCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCa                                                                                  >  1:4087507/1‑100 (MQ=255)    ctCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCaa                                                                                 >  1:5506121/1‑100 (MQ=255)      catcaGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAatt                                                                               >  1:3031741/1‑100 (MQ=255)      catcaGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAatt                                                                               >  1:3813533/1‑100 (MQ=255)               gTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCttatt                                                                      <  1:5069948/100‑1 (MQ=255)                 cctGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCttattgg                                                                    >  1:2593556/3‑98 (MQ=255)                       ttCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGtaaag                                                              >  1:4088483/1‑100 (MQ=255)                         cAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGtaaagat                                                            >  1:225687/1‑100 (MQ=255)                         cAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGtaaagat                                                            >  1:3492371/1‑100 (MQ=255)                         cAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGtaaagat                                                            >  1:4109663/1‑100 (MQ=255)                           gTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATAAGGtaagatcgg                                                          >  1:1735989/1‑94 (MQ=255)                            tAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTGTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGtacgagatcg                                                         >  1:3920663/1‑92 (MQ=255)                             gatatTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGtaaagatcgga                                                        >  1:4848958/2‑96 (MQ=255)                                 ctAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGtaaagatcggaagag                                                    >  1:3765248/2‑92 (MQ=255)                                  tAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTTGGtaaagatcggaagagc                                                   >  1:4035252/1‑91 (MQ=255)                                   aGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAat                                                  >  1:1398129/1‑100 (MQ=255)                                   aGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAat                                                  >  1:1175340/1‑100 (MQ=255)                                   aGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAAGtacagatcggaagagca                                                  >  1:70903/1‑90 (MQ=255)                                     tGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGTAAAGATAAAgagatcgga                                                >  1:1128867/1‑95 (MQ=255)                                                   aTTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAATATAAAGAAGGCATTGa                                  <  1:2400885/100‑1 (MQ=255)                                                     tataACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAATATAAAGAAGGCATTGAAg                                >  1:3600277/1‑100 (MQ=255)                                                                 aaCGTTGTACGCTTTAACCAAGAAAATCAAAAGCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAATATAAAGAAGGCATTGAGGGCAAAGATGTCt                    >  1:5220013/1‑100 (MQ=255)                                                                                   cAAGAAAATCAAAAGCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAATATAAAGAAGGCATTGAAGGCAAAGATGTCTTTGTGGTAAAAGaattaa  >  1:5589174/1‑100 (MQ=255)                                                                                                 |                                                                                      CTACTCATCAGAAAGTTTTGAATTCAGTAATATTAGTGGTAAGGTTGAAAATTATAACGGTTCTAACGTTGTACGCTTTAACCAAGAAAATCAAAATCACCAATTATTCTTATTAGGTAAAGATAAAGAGAAATATAAAGAAGGCATTGAAGGCAAAGATGTCTTTGTGGTAAAAGAATTAA  >  NC_007793/453514‑453695

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