Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,980,188	C→T	*475* (TAG→TAA)	otsA ←	trehalose‑6‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,980,188	0	C	T	100.0%	62.8 / NA	21	*475* (TAG→TAA)	otsA	trehalose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/10);  total (11/10)

GCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCA  >  NC_000913/1980052‑1980318                                                                                                                                         |                                                                                                                                   gCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTa                                                                                                                                  >  1:167222/1‑139 (MQ=255)   atcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACg                                                                                                                                <  1:458698/139‑1 (MQ=255)    tcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                               <  2:311297/139‑1 (MQ=255)       tttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAg                                                                                                                            <  2:34040/139‑1 (MQ=255)                                   aatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGCGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCgcgc                                                                                                >  2:501892/1‑139 (MQ=255)                                                  tCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                                 <  2:167222/139‑1 (MQ=255)                                                                              ttATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTATCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGc                                                     <  1:716251/139‑1 (MQ=255)                                                                                        aaGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGct                                           <  1:501892/139‑2 (MQ=255)                                                                                        aaGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCa                                           >  2:86566/1‑139 (MQ=255)                                                                                         aGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCAc                                          >  1:100938/1‑139 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcg                                                                            >  1:524797/1‑99 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcg                                                                            <  2:524797/99‑1 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGc                                            <  1:725272/131‑1 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGc                                            >  2:725272/1‑131 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                    <  2:100938/139‑1 (MQ=255)                                                                                                                tCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAatat                       >  2:403786/1‑135 (MQ=255)                                                                                                                tCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAatat                       <  1:403786/135‑1 (MQ=255)                                                                                                                   aCTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGttttt                >  1:783930/1‑139 (MQ=255)                                                                                                                     ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCa              >  2:205476/1‑139 (MQ=255)                                                                                                                                 aaCCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCa  >  1:737444/1‑139 (MQ=255)                                                                                                                                 aaCCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCa  >  2:487866/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                   GCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCA  >  NC_000913/1980052‑1980318

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