Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	383,259	C→T	R226Q (CGG→CAG)	yaiP ←	putative family 2 glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	383,259	0	C	T	100.0%	218.5 / NA	68	R226Q (CGG→CAG)	yaiP	putative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/37);  total (31/37)

ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGT  >  NC_000913/383123‑383390                                                                                                                                         |                                                                                                                                    aCCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGa                                                                                                                                   >  2:1630166/1‑139 (MQ=255)          gCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAg                                                                                                                          >  2:209157/1‑139 (MQ=255)           cTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGc                                                                                                                         >  1:867139/1‑139 (MQ=255)                  aaaGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                                  <  1:209157/139‑1 (MQ=255)                     gTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccac                                                                                                               <  1:525531/139‑1 (MQ=255)                       cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  2:573160/139‑1 (MQ=255)                       cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  1:74582/139‑1 (MQ=255)                        ttttATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCaa                                                                                                            >  1:366510/1‑139 (MQ=255)                           tATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGt                                                                                                         <  2:867139/139‑1 (MQ=255)                            aTGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                        <  2:1318347/139‑1 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             >  2:1241022/1‑118 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  1:1241022/118‑1 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                        >  1:192077/1‑139 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                        >  1:490869/1‑139 (MQ=255)                                              cgcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCt                                                                                      >  2:193947/1‑139 (MQ=255)                                              cgcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCt                                                                                      >  1:3504/1‑139 (MQ=255)                                               gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                     <  1:1836213/139‑1 (MQ=255)                                               gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                     <  1:186301/139‑1 (MQ=255)                                               gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                     <  2:3504/139‑1 (MQ=255)                                                cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  2:1379756/114‑1 (MQ=255)                                                cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             >  1:1379756/1‑114 (MQ=255)                                                cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                                    <  2:1006535/139‑1 (MQ=255)                                                cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                                    <  2:490869/139‑1 (MQ=255)                                                  tATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGa                                                                                  <  2:366510/139‑1 (MQ=255)                                                               aaCGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                          >  1:1647315/1‑134 (MQ=255)                                                               aaCGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                          <  2:1647315/134‑1 (MQ=255)                                                                      ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGc                                                                                                                         <  1:1774466/80‑1 (MQ=255)                                                                      ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGc                                                                                                                         >  2:1774466/1‑80 (MQ=255)                                                                      ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGaaa                                                              >  1:1024437/1‑139 (MQ=255)                                                                               ccACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCGCCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                     <  2:1289111/139‑1 (MQ=255)                                                                                 aCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                   <  2:350112/139‑1 (MQ=255)                                                                                       cACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATAc                                             <  1:1630166/139‑1 (MQ=255)                                                                                            gggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             >  1:1786091/1‑70 (MQ=255)                                                                                            gggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  2:1786091/70‑1 (MQ=255)                                                                                               cTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACg                                                                        <  2:35877/104‑1 (MQ=255)                                                                                               cTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACg                                                                        >  1:35877/1‑104 (MQ=255)                                                                                                  ttGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGa                                  >  1:1608151/1‑139 (MQ=255)                                                                                                    gCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAAc                                <  1:193947/139‑1 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                     <  2:423723/78‑1 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                     >  1:423723/1‑78 (MQ=255)                                                                                                               ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             >  2:1155002/1‑51 (MQ=255)                                                                                                               ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  1:1155002/51‑1 (MQ=255)                                                                                                                 tataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCAc                   <  2:1024437/139‑1 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                           <  2:402791/112‑1 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                           >  1:402791/1‑112 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCa                    >  1:1160994/1‑135 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCa                    <  2:1160994/135‑1 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTg                                                                                                                   <  2:1331747/37‑1 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTg                                                                                                                   >  1:1331747/1‑37 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                     <  2:202314/67‑1 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                     >  1:202314/1‑67 (MQ=255)                                                                                                                         gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             >  2:1411202/1‑41 (MQ=255)                                                                                                                         gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  1:1411202/41‑1 (MQ=255)                                                                                                                         gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                        >  2:710720/1‑46 (MQ=255)                                                                                                                         gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                        <  1:710720/46‑1 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                                <  2:542579/69‑1 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                                >  1:542579/1‑69 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                     >  2:254319/1‑96 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                     <  1:254319/96‑1 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGataa          >  2:1180093/1‑139 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             <  2:1560020/39‑1 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                             >  1:1560020/1‑39 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCa                                         >  2:531581/1‑107 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCa                                         <  1:531581/107‑1 (MQ=255)                                                                                                                            gCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGataata        <  2:1608151/139‑1 (MQ=255)                                                                                                                                  cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATAc                                             >  2:305167/1‑96 (MQ=255)                                                                                                                                  cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATAc                                             <  1:305167/96‑1 (MQ=255)                                                                                                                                  cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGt  >  1:1167330/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                    ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGT  >  NC_000913/383123‑383390

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