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

GCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCAGCA  >  NC_000913/1980052‑1980321                                                                                                                                         |                                                                                                                                      gCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTa                                                                                                                                     >  2:665031/1‑139 (MQ=255)   atcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACg                                                                                                                                   >  2:1020465/1‑139 (MQ=255)    tcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                                  <  1:1002076/139‑1 (MQ=255)         tCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                                  >  1:983839/1‑134 (MQ=255)         tCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGc                                                                                                                                  <  2:983839/134‑1 (MQ=255)                 gTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAg                                                                                                                     <  2:808338/139‑1 (MQ=255)                                           tttCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGaa                                                                                                                       >  1:732275/1‑111 (MQ=255)                                           tttCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGaa                                                                                                                       <  2:732275/111‑1 (MQ=255)                                             tCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCa                                                                                                                                 <  1:852950/99‑1 (MQ=255)                                             tCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCa                                                                                                                                 >  2:852950/1‑99 (MQ=255)                                                tttCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcg                                                                                      <  1:1020465/139‑1 (MQ=255)                                                   cgcaAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCgcgc                                                                                   >  1:730684/1‑139 (MQ=255)                                                     cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCtt                                                                                 >  1:366362/1‑139 (MQ=255)                                                       aGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcg                                                                               >  2:545793/1‑139 (MQ=255)                                                                    gCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGc                                                                  >  1:197007/1‑139 (MQ=255)                                                                                   gCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATg                                                   >  2:937949/1‑139 (MQ=255)                                                                                      tCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      >  1:369585/1‑117 (MQ=255)                                                                                      tCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      <  2:369585/117‑1 (MQ=255)                                                                                      tCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAg                                                <  2:844370/139‑1 (MQ=255)                                                                                        aaGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCa                                              <  1:888547/139‑1 (MQ=255)                                                                                             aGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGACTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCt                                         >  1:35872/1‑139 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGaa                                                 <  1:489416/129‑1 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGaa                                                 >  2:489416/1‑129 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                       >  1:509520/1‑139 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                       >  2:436503/1‑139 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTACTGc                                       >  2:534395/1‑139 (MQ=255)                                                                                               gTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGACTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                       <  2:35872/139‑1 (MQ=255)                                                                                                 cGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCa                                     <  2:366362/139‑1 (MQ=255)                                                                                                      tgCTGTTAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGt                                <  2:509520/139‑1 (MQ=255)                                                                                                           ttAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGc                                                                                                                <  2:506531/54‑1 (MQ=255)                                                                                                           ttAGTTCCACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGc                                                                                                                >  1:506531/1‑54 (MQ=255)                                                                                                                  cACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      <  2:229491/89‑1 (MQ=255)                                                                                                                  cACTTACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      >  1:229491/1‑89 (MQ=255)                                                                                                                     ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      >  2:312511/1‑86 (MQ=255)                                                                                                                     ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      <  1:312511/86‑1 (MQ=255)                                                                                                                     ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      <  2:451561/86‑1 (MQ=255)                                                                                                                     ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt                                                                      >  1:451561/1‑86 (MQ=255)                                                                                                                     ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                       >  2:339438/1‑117 (MQ=255)                                                                                                                     ttACGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGc                                       <  1:339438/117‑1 (MQ=255)                                                                                                                       aCGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGt                                                           >  1:880444/1‑95 (MQ=255)                                                                                                                       aCGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGt                                                           <  2:880444/95‑1 (MQ=255)                                                                                                                       aCGGGAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACg               >  1:10996/1‑139 (MQ=255)                                                                                                                         gggAGATTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGAt             >  2:302987/1‑139 (MQ=255)                                                                                                                           gagaTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCtt                                                                <  1:305832/86‑1 (MQ=255)                                                                                                                           gagaTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCtt                                                                >  2:305832/1‑86 (MQ=255)                                                                                                                             gaTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACg         <  1:56074/139‑1 (MQ=255)                                                                                                                             gaTTAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACATTATCGCGCTGCTGGCTTTACGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACg         >  2:709028/1‑139 (MQ=255)                                                                                                                               ttAACCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTGTCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTc       <  1:534395/139‑1 (MQ=255)                                                                                                                                  aCCGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCAg    <  1:436503/139‑1 (MQ=255)                                                                                                                                   ccGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGa                                                                           <  1:217426/67‑1 (MQ=255)                                                                                                                                   ccGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGa                                                                           >  2:217426/1‑67 (MQ=255)                                                                                                                                   ccGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCa     >  2:786282/1‑137 (MQ=255)                                                                                                                                   ccGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCa     <  1:786282/137‑1 (MQ=255)                                                                                                                                   ccGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCAGc   >  1:554942/1‑139 (MQ=255)                                                                                                                                    cGCTCTTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCAGCa  <  1:545793/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      GCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTAGGTCGCTAATGAAGCACTCCTGCCAGTGGTTAATATCGTTTTTCACGATAACGTCCAGCA  >  NC_000913/1980052‑1980321

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