Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	518,721	C→A	M143I (ATG→ATT)	ybbO ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	518,721	0	C	A	100.0%	227.1 / NA	69	M143I (ATG→ATT)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (39/30);  total (39/30)

TCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATA  >  NC_000913/518587‑518852                                                                                                                                       |                                                                                                                                    tCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtaatca                                                                                                                                 <  1:1719199/139‑1 (MQ=255)   gATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATcaca                                                                                                                               <  2:456617/139‑1 (MQ=255)         gCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             <  1:632414/135‑1 (MQ=255)         gCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             >  2:632414/1‑135 (MQ=255)         gCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGc                                                                                                                         >  1:352476/1‑139 (MQ=255)             aCTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtaat                                                                                                                                   >  1:813601/1‑125 (MQ=255)               tttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACTACGCCTCCAGCGCATATTTACTGGCCGCGTACGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                   <  1:1337649/139‑1 (MQ=255)                    ttCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAAT‑CCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGc                                                                                                             >  2:806283/1‑139 (MQ=255)                            tgGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             >  1:820637/1‑116 (MQ=255)                            tgGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             <  2:820637/116‑1 (MQ=255)                                            gcAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgcg                                                                                      <  1:554600/139‑1 (MQ=255)                                              aGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             <  1:111191/98‑1 (MQ=255)                                              aGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             >  2:111191/1‑98 (MQ=255)                                                     cTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAAACGCCCTTCACCGTg                                                                                                              >  1:570476/1‑106 (MQ=255)                                                     cTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAAACGCCCTTCACCGTg                                                                                                              <  2:570476/106‑1 (MQ=255)                                                        aCCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTg                                                                          >  2:1352035/1‑139 (MQ=255)                                                               cTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTgcg                                                                   >  2:1086936/1‑139 (MQ=255)                                                                  cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             <  1:1539423/78‑1 (MQ=255)                                                                  cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             >  2:1539423/1‑78 (MQ=255)                                                                  cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCg                                                                >  2:1552998/1‑139 (MQ=255)                                                                  cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCg                                                                >  2:750940/1‑139 (MQ=255)                                                                       cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                <  1:1378517/86‑1 (MQ=255)                                                                       cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                >  2:1378517/1‑86 (MQ=255)                                                                                  gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTa                                                                                                         >  2:122170/1‑82 (MQ=255)                                                                                  gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTa                                                                                                         <  1:122170/82‑1 (MQ=255)                                                                                    cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                      >  1:41220/1‑83 (MQ=255)                                                                                    cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                      <  2:1385804/83‑1 (MQ=255)                                                                                    cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                      <  2:41220/83‑1 (MQ=255)                                                                                    cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                      >  1:1385804/1‑83 (MQ=255)                                                                                    cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                              >  2:1403002/1‑139 (MQ=255)                                                                                    cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                              >  2:1209165/1‑139 (MQ=255)                                                                                    cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                              >  1:1922869/1‑139 (MQ=255)                                                                                     gcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                             >  2:1428858/1‑139 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCAc                                                                                                                  >  1:500533/1‑69 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCAc                                                                                                                  <  2:500533/69‑1 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCg                                                    <  1:1224612/131‑1 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCg                                                    >  2:1224612/1‑131 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActg                                            >  2:384238/1‑139 (MQ=255)                                                                                       gTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgc                                           >  1:64683/1‑139 (MQ=255)                                                                                        tAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgct                                          >  2:239896/1‑139 (MQ=255)                                                                                               cACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             >  1:1950495/1‑49 (MQ=255)                                                                                               cACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             <  2:1950495/49‑1 (MQ=255)                                                                                                 cGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCt                                 >  1:1724982/1‑139 (MQ=255)                                                                                                 cGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCt                                 >  2:84485/1‑139 (MQ=255)                                                                                                      cGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCa                            <  1:239896/139‑1 (MQ=255)                                                                                                        gCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGt                                                                                              >  1:9086/1‑71 (MQ=255)                                                                                                        gCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGt                                                                                              <  2:9086/71‑1 (MQ=255)                                                                                                        gCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACg                          <  1:84485/139‑1 (MQ=255)                                                                                                         cGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACgg                         >  2:1610694/1‑139 (MQ=255)                                                                                                          gTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGc                        <  1:1097692/139‑1 (MQ=255)                                                                                                             gagaTTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGa                     <  2:352476/139‑1 (MQ=255)                                                                                                                  tAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtggtg                >  1:1812356/1‑139 (MQ=255)                                                                                                                   aaTCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                   <  1:1017529/39‑1 (MQ=255)                                                                                                                   aaTCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                   >  2:1017529/1‑39 (MQ=255)                                                                                                                    aatccATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGGGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGa              <  1:806283/136‑1 (MQ=255)                                                                                                                      cccATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTACGCACGGCTGATGGTGGaaa            <  1:1209165/139‑1 (MQ=255)                                                                                                                       ccATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAg           <  2:64683/139‑1 (MQ=255)                                                                                                                        cATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAgg                                                                                        <  1:598462/61‑1 (MQ=255)                                                                                                                        cATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAgg                                                                                        >  2:598462/1‑61 (MQ=255)                                                                                                                           cACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAggggg       <  2:1922869/139‑1 (MQ=255)                                                                                                                              cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGa              >  2:1339531/1‑129 (MQ=255)                                                                                                                              cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGa              <  1:1339531/129‑1 (MQ=255)                                                                                                                              cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCa    <  1:384238/139‑1 (MQ=255)                                                                                                                              cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCa    >  2:147412/1‑139 (MQ=255)                                                                                                                               gatgatGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTg                      <  1:2028054/120‑1 (MQ=255)                                                                                                                               gatgatGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTg                      >  2:2028054/1‑120 (MQ=255)                                                                                                                               gatgatGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCAt   >  1:1002831/1‑139 (MQ=255)                                                                                                                               gatgatGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCAt   >  1:1033250/1‑139 (MQ=255)                                                                                                                                atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATa  >  1:551560/1‑139 (MQ=255)                                                                                                                                   atgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTgcgc                             >  1:1357967/1‑109 (MQ=255)                                                                                                                                   atgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTgcgc                             <  2:1357967/109‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    TCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATA  >  NC_000913/518587‑518852

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