Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,982,132	G→A	R86C (CGC→TGC)	otsB ←	trehalose‑6‑phosphate phosphatase, biosynthetic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,982,132	0	G	A	100.0%	124.0 / NA	32	R86C (CGC→TGC)	otsB	trehalose‑6‑phosphate phosphatase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/19);  total (13/19)

TGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCGCTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCTTGCAGAATATTGTC  >  NC_000913/1981984‑1982273                                                                                                                                                     |                                                                                                                                              tGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACt                                                                                                                                             >  1:553075/1‑151 (MQ=255)   aCGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTcc                                                                                                                                           <  1:471080/151‑1 (MQ=255)     gATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGc                                                                                                                                         >  1:562391/1‑151 (MQ=255)         aTGCGGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCa                                                                                                                                     <  2:484708/151‑1 (MQ=255)             aGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATg                                                                                                                                   <  1:393649/149‑1 (MQ=255)             aGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCa                                                                                                                                 <  1:89739/151‑1 (MQ=255)              gcgcAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCAc                                                                                                                                <  1:541392/151‑1 (MQ=255)                           cccTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACggg                                                                                                                   >  2:171157/1‑151 (MQ=255)                           cccTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACggg                                                                                                                   >  2:311934/1‑151 (MQ=255)                             cTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGa                                                                                                                  >  2:120928/1‑150 (MQ=255)                                   ccTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGAt                                                                                                            >  1:642295/1‑150 (MQ=255)                                                gcCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGt                                                                                              >  2:46103/1‑151 (MQ=255)                                                              cGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGACCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTccac                                                                                >  2:460633/1‑151 (MQ=255)                                                                    cTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTg                                                                           <  2:553059/150‑1 (MQ=255)                                                                     tGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTg                                                                           >  2:421858/1‑149 (MQ=255)                                                                                gCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGata                                                               <  1:81368/150‑1 (MQ=255)                                                                                  aCGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATca                                                            <  1:288296/151‑1 (MQ=255)                                                                                   cGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATcaa                                                           <  1:75578/151‑1 (MQ=255)                                                                                    gCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGTCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATcaat                                                          <  1:284182/151‑1 (MQ=255)                                                                                       aaTATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATcaatgcc                                                       >  2:220351/1‑151 (MQ=255)                                                                                       aaTATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATcaatgc                                                        <  2:589448/150‑1 (MQ=255)                                                                                            cACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATcaatgccaatgc                                                  <  1:604065/151‑1 (MQ=255)                                                                                                       gCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTg                                       >  2:354780/1‑151 (MQ=255)                                                                                                        cATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGc                                      >  2:520758/1‑151 (MQ=255)                                                                                                                aGATGAACGATATGTGTTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGcc                               >  1:458510/1‑150 (MQ=255)                                                                                                                             gtgtTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTcc                 <  1:171162/151‑1 (MQ=255)                                                                                                                               gtTTTACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCtt               <  1:44015/151‑1 (MQ=255)                                                                                                                                   tACCATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCTTGCa            <  2:109839/150‑1 (MQ=255)                                                                                                                                     ccATTGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCTTGCAGaa         <  2:339953/151‑1 (MQ=255)                                                                                                                                         tGATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCTTGCAGAatat      <  2:208103/150‑1 (MQ=255)                                                                                                                                          gATGTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCTTGCAGAATAtt     <  2:298311/150‑1 (MQ=255)                                                                                                                                             gTCACGGCACTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCTTGCAGAATATTGTc  <  1:658710/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              TGACGATAATGCAGCGCAAAAGCCATCCCTTTCGCCTCCAGCTCCGCGCCGGGATACTGAGCGATGACTGTATGCAGTTGCACGCTAATATCACGCGCAATCGCATCCGGCAGATGAACGATATGTGTTTTACCATTGATGTCACGGCGCTCCGCCCCATGCACGCCCGCTAACGGGAAGCGATAAGGTTTTGCCAGTGCGTCAAGCTCCACCATTGAGCGCCCTGATATCAATGCCAATGCACCATCACTTGCGGTTGCCAGTAGCTGTAGTCCTTGCAGAATATTGTC  >  NC_000913/1981984‑1982273

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