Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	184,257	C→T	*271* (TAG→TAA)	yaeI ←	phosphodiesterase with model substrate bis‑pNPP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	184,257	0	C	T	100.0%	242.5 / NA	74	*271* (TAG→TAA)	yaeI	phosphodiesterase with model substrate bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/34);  total (40/34)

GGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTT  >  NC_000913/184143‑184381                                                                                                                   |                                                                                                                             ggcctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACt                                                                                                      >  2:646837/1‑139 (MQ=255)      ggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCgg                                                                                                 >  1:1600437/1‑139 (MQ=255)      ggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCgg                                                                                                 >  1:1398633/1‑139 (MQ=255)          cgcccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCTTGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGGTAAAATTCACCATTTACACCAGTTCCAGCATCGTTACTTCCGggcgg                                                                                             <  2:1307416/139‑1 (MQ=255)            cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTaca                                                                                                                          <  2:162460/108‑1 (MQ=255)            cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTaca                                                                                                                          >  1:162460/1‑108 (MQ=255)            cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGggcgg                                                                                             >  1:541154/1‑137 (MQ=255)            cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGggcgg                                                                                             <  2:541154/137‑1 (MQ=255)            cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCa                                                                                           >  1:1203373/1‑139 (MQ=255)            cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGCTCCAGCATCGTTACTTCCGGGCGGCa                                                                                           >  1:1307416/1‑139 (MQ=255)                 cTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCa                                                                                      >  1:833063/1‑139 (MQ=255)                      cTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGc                                                                                 <  1:2008242/139‑1 (MQ=255)                         tgATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACg                                                                                  >  2:1112468/1‑135 (MQ=255)                         tgATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACg                                                                                  <  1:1112468/135‑1 (MQ=255)                           aTCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                               <  1:1472095/136‑1 (MQ=255)                           aTCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                               >  2:1472095/1‑136 (MQ=255)                              gccgtaaGCCTGTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATCACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATa                                                                         >  1:1193903/6‑139 (MQ=255)                              gCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                               >  2:963312/1‑101 (MQ=255)                              gCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                               <  1:963312/101‑1 (MQ=255)                                aGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCa                                                                                >  2:821719/1‑130 (MQ=255)                                aGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCa                                                                                <  1:821719/130‑1 (MQ=255)                                       cTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGcc                                                                <  2:1193903/139‑1 (MQ=255)                                               tcCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccac                                                        <  1:646837/139‑1 (MQ=255)                                                ccAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccacg                                                       <  2:1600437/139‑1 (MQ=255)                                                  aGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                     <  2:1203373/139‑1 (MQ=255)                                                  aGACTTCTATCATGGCCTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATTCAAACTGCCCACGCCAAggg                                                     <  2:1398633/139‑1 (MQ=255)                                                         tATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgcc                                                          >  2:1829944/1‑127 (MQ=255)                                                         tATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgcc                                                          <  1:1829944/127‑1 (MQ=255)                                                              tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                               >  2:1040841/1‑69 (MQ=255)                                                              tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                               <  1:1040841/69‑1 (MQ=255)                                                              tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAtgt                                         >  2:889624/1‑139 (MQ=255)                                                              tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAtgt                                         >  1:925371/1‑139 (MQ=255)                                                                gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCa                                                                                                              >  1:859351/1‑68 (MQ=255)                                                                gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCa                                                                                                              <  2:859351/68‑1 (MQ=255)                                                                gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                               <  1:946580/99‑1 (MQ=255)                                                                gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                               >  2:946580/1‑99 (MQ=255)                                                                gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATgtgt                                       >  1:1429461/1‑139 (MQ=255)                                                                  gTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGtt                                                                                        <  1:882257/88‑1 (MQ=255)                                                                  gTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGtt                                                                                        >  2:882257/1‑88 (MQ=255)                                                                   tAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCtt                                    >  2:1722531/1‑139 (MQ=255)                                                                   tAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCtt                                    >  2:1023359/1‑139 (MQ=255)                                                                   tAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCtt                                    >  1:190351/1‑139 (MQ=255)                                                                        aTCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGcc                               <  2:190351/139‑1 (MQ=255)                                                                         tCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCa                              <  1:1722531/139‑1 (MQ=255)                                                                           aaaaaaTTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCaaa                            <  2:833063/139‑1 (MQ=255)                                                                                        cTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCg               <  2:925371/139‑1 (MQ=255)                                                                                              tCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACg         >  1:1341604/1‑139 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACg                                                                                  <  2:1509454/62‑1 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACg                                                                                  >  1:1509454/1‑62 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                               >  1:535760/1‑65 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                               <  2:535760/65‑1 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCa                                                                                >  2:555832/1‑64 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCa                                                                                <  1:555832/64‑1 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaaa                                                                     >  1:921538/1‑75 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaaa                                                                     <  2:921538/75‑1 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                     >  1:748475/1‑91 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                     <  2:748475/91‑1 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTaa                   <  2:1859740/125‑1 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTaa                   >  1:1859740/1‑125 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCgg              >  1:1648718/1‑130 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCgg              <  2:1648718/130‑1 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAAc     >  2:854463/1‑139 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAAc     >  1:1055340/1‑139 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAAc     >  1:615297/1‑139 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCA‑TTCCAGCATCGTTACTTCCGGGCGGCAGtt                                                                                        >  1:13892/1‑55 (MQ=255)                                                                                                  gTGATAAAATTCACGATTTACACCA‑TTCCAGCATCGTTACTTCCGGGCGGCAGtt                                                                                        <  2:13892/55‑1 (MQ=255)                                                                                                    gATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGt   <  2:1429461/139‑1 (MQ=255)                                                                                                     aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGtt  >  2:429356/1‑139 (MQ=255)                                                                                                     aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACAGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCaa                             <  2:1979857/112‑1 (MQ=255)                                                                                                     aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACAGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCaa                             >  1:1979857/1‑112 (MQ=255)                                                                                                          aTTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACgg                                                      >  1:1867941/1‑82 (MQ=255)                                                                                                          aTTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACgg                                                      <  2:1867941/82‑1 (MQ=255)                                                                                                            tCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGc                                                                                 <  2:1873357/53‑1 (MQ=255)                                                                                                            tCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGc                                                                                 >  1:1873357/1‑53 (MQ=255)                                                                                                                   |                                                                                                                             GGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTT  >  NC_000913/184143‑184381

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