Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,530,076	(A)7→8	coding (155/483 nt)	ydcD →	putative immunity protein for RhsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,530,069	1	.	A	100.0%	238.7 / NA	58	E50K (GAA→AAA)	ydcD	putative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (10/48);  total (10/48)

GCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAG  >  NC_000913/1529924‑1530214                                                                                                                                                   |                                                                                                                                                  gCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑aaaaa                                                                                                                                              >  2:695130/1‑151 (MQ=255)                   cAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaat                                                                                                                             <  2:391731/151‑1 (MQ=255)                                         aCCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                                                                       >  2:749616/1‑151 (MQ=255)                                          cccTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                                                                      >  1:295413/1‑151 (MQ=255)                                           ccTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                                                                     <  2:356889/151‑1 (MQ=255)                                           ccTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                                                                     <  2:1024161/151‑1 (MQ=255)                                                 tttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                                                                               <  2:1314458/151‑1 (MQ=255)                                                     tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                                                                           <  1:828995/151‑1 (MQ=255)                                                     tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                                                                           <  2:295499/151‑1 (MQ=255)                                                        gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                                                                        >  1:1328085/1‑151 (MQ=255)                                                           gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                                                                     >  2:981090/1‑151 (MQ=255)                                                                ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAACTTAGTAGAGGGAATGGTCATTTTATAAg                                                                                >  1:197206/1‑151 (MQ=255)                                                                ccTTAGTTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTTAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGATCATTTTATaac                                                                                <  1:569555/151‑2 (MQ=255)                                                                        ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                                                                        <  2:1135432/151‑1 (MQ=255)                                                                        ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                                                                        <  1:1403970/151‑1 (MQ=255)                                                                        ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                                                                        <  1:1295937/151‑1 (MQ=255)                                                                        ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                                                                        <  1:36337/151‑1 (MQ=255)                                                                               cAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa                                                                 <  1:1240077/151‑1 (MQ=255)                                                                                   tCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc                                                             <  2:1149129/151‑1 (MQ=255)                                                                                   tCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc                                                             <  1:805460/151‑1 (MQ=255)                                                                                     aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg                                                           <  2:1162515/151‑1 (MQ=255)                                                                                     aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg                                                           <  1:962752/151‑1 (MQ=255)                                                                                     aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg                                                           <  1:168255/151‑1 (MQ=255)                                                                                               aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACtgt                                                 <  2:14115/151‑1 (MQ=255)                                                                                                 cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTgg                                               <  2:538597/151‑1 (MQ=255)                                                                                                 cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTgg                                               <  1:1281372/151‑1 (MQ=255)                                                                                                  gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGa                                              >  2:1262226/1‑151 (MQ=255)                                                                                                   tATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAg                                             <  2:1186617/151‑1 (MQ=255)                                                                                                   tATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGAAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTGCAGGACGCAGGGTGGAg                                             <  1:354841/151‑1 (MQ=255)                                                                                                      tttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTgg                                          <  2:382267/151‑1 (MQ=255)                                                                                                      tttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTgg                                          <  1:1378103/151‑1 (MQ=255)                                                                                                      tttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTgg                                          <  1:1333877/151‑1 (MQ=255)                                                                                                      tttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTgg                                          <  1:245738/151‑1 (MQ=255)                                                                                                      tttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTgg                                          <  1:383737/151‑1 (MQ=255)                                                                                                      tttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTgg                                          <  1:57272/151‑1 (MQ=255)                                                                                                      tttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATGAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGtag                                          <  1:634583/151‑3 (MQ=255)                                                                                                             aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGatat                                   >  2:623023/1‑151 (MQ=255)                                                                                                             aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGatat                                   <  2:572012/151‑1 (MQ=255)                                                                                                             aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGatat                                   <  2:834800/151‑1 (MQ=255)                                                                                                             aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGatat                                   <  1:1394540/151‑1 (MQ=255)                                                                                                             aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGatat                                   <  1:1368488/151‑1 (MQ=255)                                                                                                                  aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATc                              <  2:489245/151‑1 (MQ=255)                                                                                                                    acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCtt                            >  1:271458/1‑151 (MQ=255)                                                                                                                                   gagaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGat             >  2:1165916/1‑151 (MQ=255)                                                                                                                                     gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGatat           >  1:677332/1‑151 (MQ=255)                                                                                                                                       tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTc         <  1:1168741/151‑1 (MQ=255)                                                                                                                                        atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCt        <  2:204038/151‑1 (MQ=255)                                                                                                                                            cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTtat    <  2:185241/151‑1 (MQ=255)                                                                                                                                            cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTtat    <  2:968621/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGTGAATGGTCATTTTATAAGTATAATGGTTAATAATAACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  1:233335/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACATGACGCACTGTGAAGTGGTTGATATAAATCTTAATTTATTCCATGATATTCTTTATAg  <  1:623438/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  2:294771/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  2:518662/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  2:551282/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  1:149093/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  2:633071/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  2:1047804/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  1:1233338/151‑1 (MQ=255)                                                                                                                                              ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAg  <  1:912704/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAG  >  NC_000913/1529924‑1530214

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