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%	133.1 / NA	42	*475* (TAG→TAA)	otsA	trehalose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/22);  total (20/22)

GCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACG  >  NC_000913/1980052‑1980314                                                                                                                                         |                                                                                                                               gCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTa                                                                                                                              >  2:886432/1‑139 (MQ=255)     catTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCa                                                                                                                          <  1:1003107/139‑1 (MQ=255)            aTTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCa                                                                                                                          <  1:656575/132‑1 (MQ=255)            aTTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCa                                                                                                                          >  2:656575/1‑132 (MQ=255)                                   aatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCgcgc                                                                                            <  2:649830/139‑1 (MQ=255)                                      cacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctgc                                                                                         <  1:32108/139‑1 (MQ=255)                                        caGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGctgctg                                                                                       >  2:850904/1‑139 (MQ=255)                                                     cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCtt                                                                          >  2:1038838/1‑139 (MQ=255)                                                              tGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACt                                                                 >  1:651942/1‑139 (MQ=255)                                                               gCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTa                                                                <  2:127318/139‑1 (MQ=255)                                                                   ggCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTg                                                            >  2:955148/1‑139 (MQ=255)                                                                           gagTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGt                                                    <  1:1101980/139‑1 (MQ=255)                                                                             gTTATTGCTTCAAGTGAGGTCGATGTGCTGGTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTGTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGGTTTAGGTCg                                                  <  1:667327/139‑1 (MQ=255)                                                                                      tCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAg                                         >  1:1245804/1‑139 (MQ=255)                                                                                      tCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAg                                         >  2:82419/1‑139 (MQ=255)                                                                                       cAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGc                                        <  1:695167/139‑1 (MQ=255)                                                                                        aaGTGATGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCa                                       <  1:987815/139‑1 (MQ=255)                                                                                          gTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACt                                     <  2:1245804/139‑1 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                <  1:82419/139‑1 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                >  1:5181/1‑139 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                >  1:1122580/1‑139 (MQ=255)                                                                                                 cGATGTGCTGTTATTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCa                              <  1:1020490/139‑1 (MQ=255)                                                                                                 cGATGTGCTGTGAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCa                              <  1:850904/139‑1 (MQ=255)                                                                                                    tgtgCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                           >  2:183631/1‑43 (MQ=255)                                                                                                    tgtgCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                           <  1:183631/43‑1 (MQ=255)                                                                                                       gCTGTTAGTTCCACTTACGGGAGGTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGtt                        >  1:374297/1‑139 (MQ=255)                                                                                                        cTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                >  1:62362/1‑130 (MQ=255)                                                                                                        cTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                <  2:62362/130‑1 (MQ=255)                                                                                                              gTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCg                 <  1:377213/139‑1 (MQ=255)                                                                                                                  cACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                             <  2:827662/75‑1 (MQ=255)                                                                                                                  cACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                             >  1:827662/1‑75 (MQ=255)                                                                                                                  cACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGtttt             >  1:334267/1‑139 (MQ=255)                                                                                                                         gggAGGTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCTGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGAt      <  2:374297/139‑1 (MQ=255)                                                                                                                             gaTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACg  >  1:168774/1‑139 (MQ=255)                                                                                                                              aTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCg                 >  1:1123239/1‑123 (MQ=255)                                                                                                                              aTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCg                 <  2:1123239/123‑1 (MQ=255)                                                                                                                                tAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATcg                                                                                               >  1:1210508/1‑43 (MQ=255)                                                                                                                                tAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATcg                                                                                               <  2:1210508/43‑1 (MQ=255)                                                                                                                                     gCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                                                                                          >  2:1219002/1‑43 (MQ=255)                                                                                                                                     gCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                                                                                          <  1:1219002/43‑1 (MQ=255)                                                                                                                                     gCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                             >  1:1163356/1‑56 (MQ=255)                                                                                                                                     gCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                             <  2:1163356/56‑1 (MQ=255)                                                                                                                                         |                                                                                                                               GCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACG  >  NC_000913/1980052‑1980314

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