Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,513,877	A→G	K16R (AAG→AGG)	ydcU →	putative ABC transporter membrane subunit YdcU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,513,877	0	A	G	100.0%	199.3 / NA	55	K16R (AAG→AGG)	ydcU	putative ABC transporter membrane subunit YdcU
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/26);  total (29/26)

GAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGA  >  NC_000913/1513744‑1514012                                                                                                                                      |                                                                                                                                        gaaCTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTc                                                                                                                                   <  2:283785/140‑1 (MQ=255)     tgcctgccACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGt                                                                                                                               <  2:1295365/140‑1 (MQ=255)            cACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTtcttct                                                                                                                        >  2:1289567/1‑140 (MQ=255)                cTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCa                                                                                                                    <  2:454407/140‑1 (MQ=255)                     gCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATc                                                                                                               >  2:1254494/1‑140 (MQ=255)                          gACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCgggg                                                                                                          >  1:940393/1‑140 (MQ=255)                           acaacaGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCggggc                                                                                                         <  1:1123062/140‑1 (MQ=255)                           acaacaGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCggggc                                                                                                         <  2:416295/140‑1 (MQ=255)                               caGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCggggctggg                                                                                                     >  1:378616/1‑140 (MQ=255)                               caGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCggggctggg                                                                                                     >  2:1387184/1‑140 (MQ=255)                                 ggTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCggggctggggc                                                                                                   >  2:588196/1‑140 (MQ=255)                                     aTGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCggggctgggg                                                                                                    <  1:134167/135‑1 (MQ=255)                                     aTGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCggggctgggg                                                                                                    >  2:134167/1‑135 (MQ=255)                                     aTGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGtt                                                                                               >  2:1390742/1‑140 (MQ=255)                                                 tCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTgctgct                                                                                   <  1:1254494/140‑1 (MQ=255)                                                    cgTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTgg                                                                                <  2:581329/140‑1 (MQ=255)                                                      tgatgtgatgGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTTCTTGGcc                                                                              <  2:940393/140‑1 (MQ=255)                                                             atgGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAAt                                                                       >  1:1074076/1‑140 (MQ=255)                                                                     gCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAAtg                                                                      <  1:178916/133‑1 (MQ=255)                                                                     gCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAAtg                                                                      >  2:178916/1‑133 (MQ=255)                                                                     gCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTg                                                               >  1:1134547/1‑140 (MQ=255)                                                                     gCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTg                                                               <  2:176255/140‑1 (MQ=255)                                                                        ggTTGAGGAGAGGTGCATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCa                                                            <  1:1343283/140‑1 (MQ=255)                                                                        ggTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCa                                                            <  1:588196/140‑1 (MQ=255)                                                                          ttGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCAtt                                                          >  2:200960/1‑140 (MQ=255)                                                                           tGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCa                                                            <  1:871337/137‑1 (MQ=255)                                                                           tGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCa                                                            >  2:871337/1‑137 (MQ=255)                                                                           tGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTg                                                         >  1:1294756/1‑140 (MQ=255)                                                                            gaggagAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGt                                                        <  2:1266564/140‑1 (MQ=255)                                                                                   ggTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTc                                                 >  2:1156020/1‑140 (MQ=255)                                                                                    gTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCg                                                >  1:673000/1‑140 (MQ=255)                                                                                       aaTGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGACTGGGCAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTAGGCt                                             >  2:1272728/1‑140 (MQ=255)                                                                                         tGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCg                                           >  1:1266213/1‑140 (MQ=255)                                                                                           gCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCgct                                         <  1:1289567/140‑1 (MQ=255)                                                                                            cGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCgctg                                        >  1:2133/1‑140 (MQ=255)                                                                                              aTGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCgc                                          <  2:952787/136‑1 (MQ=255)                                                                                              aTGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCgc                                          >  1:952787/1‑136 (MQ=255)                                                                                                 aaTGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGac                                   <  2:1074076/140‑1 (MQ=255)                                                                                                  aTGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGaca                                  <  1:1390742/140‑1 (MQ=255)                                                                                                        tGCAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTa                            >  1:409864/1‑140 (MQ=255)                                                                                                          cAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGa                                    <  2:831336/130‑1 (MQ=255)                                                                                                          cAATCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGa                                    >  1:831336/1‑130 (MQ=255)                                                                                                           aaTCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATgg                         >  2:413212/1‑140 (MQ=255)                                                                                                            aTCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGc                        <  2:1134547/140‑1 (MQ=255)                                                                                                            aTCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGc                        >  2:1129985/1‑140 (MQ=255)                                                                                                             tCACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCa                       >  2:772410/1‑140 (MQ=255)                                                                                                              cACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      <  2:378616/140‑1 (MQ=255)                                                                                                              cACCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      <  2:1294756/140‑1 (MQ=255)                                                                                                               aCCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      >  2:465512/1‑139 (MQ=255)                                                                                                               aCCTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      <  1:465512/139‑1 (MQ=255)                                                                                                                 cTTCACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGa                   >  1:250028/1‑140 (MQ=255)                                                                                                                    cACGTCCAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGttt                >  2:1324633/1‑140 (MQ=255)                                                                                                                         ccAGGTCTGGGTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATAc           >  1:173362/1‑140 (MQ=255)                                                                                                                                 gggTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTTACg   >  1:417950/1‑140 (MQ=255)                                                                                                                                 gggTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACg   <  2:1266213/140‑1 (MQ=255)                                                                                                                                  ggTAGGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTAGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGa  <  1:1272728/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        GAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGA  >  NC_000913/1513744‑1514012

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