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

ATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTC  >  NC_000913/1980054‑1980316                                                                                                                                       |                                                                                                                                 atcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACg                                                                                                                              <  2:39765/139‑1 (MQ=255)             tGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGaa                                                                                                                  >  2:524290/1‑139 (MQ=255)                            atcacaatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                             >  2:662021/1‑113 (MQ=255)                            atcacaatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                             <  1:662021/113‑1 (MQ=255)                                   tcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAAATTTATCGCgctg                                                                                            <  1:1033269/139‑1 (MQ=255)                                      caGTTTCTTTTCGCGAGCATCCTGCTCGGTGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCgcgc                                                                                              >  1:1348320/1‑134 (MQ=255)                                      caGTTTCTTTTCGCGAGCATCCTGCTCGGTGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCgcgc                                                                                              <  2:1348320/134‑1 (MQ=255)                                        gTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                                                                                            >  2:603686/1‑134 (MQ=255)                                        gTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                                                                                            <  1:603686/134‑1 (MQ=255)                                           tCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCttt                                                                                    >  2:404054/1‑139 (MQ=255)                                                tCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                               <  1:964500/139‑1 (MQ=255)                                                         tCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGa                                                                      >  2:427461/1‑139 (MQ=255)                                                              cTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                 <  1:427461/139‑1 (MQ=255)                                                                           gTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCg                                                    <  2:1057227/139‑1 (MQ=255)                                                                                  ccTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGa                                             >  1:1252370/1‑139 (MQ=255)                                                                                    tCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAg                                           >  2:977506/1‑139 (MQ=255)                                                                                          gAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTcc                                     >  1:68345/1‑139 (MQ=255)                                                                                             gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                  <  1:524290/139‑1 (MQ=255)                                                                                                  tgtgCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAg                                                                                                            <  1:964917/60‑1 (MQ=255)                                                                                                  tgtgCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAg                                                                                                            >  2:964917/1‑60 (MQ=255)                                                                                                  tgtgCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTg                             <  2:68345/139‑1 (MQ=255)                                                                                                      cTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTa                         <  2:1252370/139‑1 (MQ=255)                                                                                                       tGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTa                                                 >  1:115014/1‑114 (MQ=255)                                                                                                       tGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTa                                                 <  2:115014/114‑1 (MQ=255)                                                                                                              tCCACTTACGGGAGTTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGATGGCTTTCCGCCCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTGATGCCAGTGGTTAATATCGtt                 >  1:1300132/1‑139 (MQ=255)                                                                                                              tCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGa                                                                                                                   <  1:697676/41‑1 (MQ=255)                                                                                                              tCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGa                                                                                                                   >  2:697676/1‑41 (MQ=255)                                                                                                              tCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGc                                                                                                           <  1:371220/49‑1 (MQ=255)                                                                                                              tCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGc                                                                                                           >  2:371220/1‑49 (MQ=255)                                                                                                              tCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGtt                 >  1:60805/1‑139 (MQ=255)                                                                                                                   ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAg                                           >  2:958948/1‑108 (MQ=255)                                                                                                                   ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAg                                           <  1:958948/108‑1 (MQ=255)                                                                                                                     aCGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                               <  1:116869/70‑1 (MQ=255)                                                                                                                     aCGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgcg                                                                               >  2:116869/1‑70 (MQ=255)                                                                                                                           gaTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACg    <  1:404054/139‑1 (MQ=255)                                                                                                                             ttAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                                                                                            >  1:257302/1‑49 (MQ=255)                                                                                                                             ttAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                                                                                            <  2:257302/49‑1 (MQ=255)                                                                                                                             ttAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTc  >  2:422347/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                 ATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTC  >  NC_000913/1980054‑1980316

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