Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,165,614	C→A	intergenic (+93/‑54)	yegP → / → yegQ	DUF1508 domain‑containing protein YegP/putative peptidase YegQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,165,614	0	C	A	100.0%	184.1 / NA	52	intergenic (+93/‑54)	yegP/yegQ	DUF1508 domain‑containing protein YegP/putative peptidase YegQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/30);  total (22/30)

CAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGC  >  NC_000913/2165493‑2165744                                                                                                                          |                                                                                                                                   caAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTg                                                                                                                  >  2:1031825/1‑140 (MQ=255) caAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTg                                                                                                                  >  2:1271844/1‑140 (MQ=255)    gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAgg                                                                                                                    >  1:154750/1‑135 (MQ=255)    gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAgg                                                                                                                    <  2:154750/135‑1 (MQ=255)    gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAg                                                                                                               <  2:607982/140‑1 (MQ=255)    gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAg                                                                                                               <  1:297394/140‑1 (MQ=255)      cAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCg                                                                                                             <  1:1271844/140‑1 (MQ=255)      cAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCg                                                                                                             <  2:619388/140‑1 (MQ=255)      cAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCg                                                                                                             <  1:262445/140‑1 (MQ=255)      cAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCg                                                                                                             <  1:346950/140‑1 (MQ=255)       aGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCgg                                                                                                            >  2:73129/1‑140 (MQ=255)        gACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGa                                                                                                           >  1:49008/1‑140 (MQ=255)         aCGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGAt                                                                                                          >  2:1081715/1‑140 (MQ=255)             tGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGa                                                                                                      <  1:1141271/140‑1 (MQ=255)              gAAAGACAATACGTAATTACGGTGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGAc                                                                                                     <  2:1493514/140‑1 (MQ=255)               aaaGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGAcc                                                                                                    >  2:887665/1‑140 (MQ=255)                 aGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTg                                                                                                  <  2:49008/140‑1 (MQ=255)                       tACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtcatca                                                                                            >  2:206829/1‑140 (MQ=255)                            aaTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACg                                                                                       <  2:106507/140‑1 (MQ=255)                                 cGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagaga                                                                                  >  1:1060823/1‑140 (MQ=255)                                 cGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagaga                                                                                  >  2:346/1‑140 (MQ=255)                                 cGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagaga                                                                                  >  2:729462/1‑140 (MQ=255)                                   aTGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAAt                                                                                >  2:501633/1‑140 (MQ=255)                                   aTGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAAt                                                                                >  1:139820/1‑140 (MQ=255)                                           cGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTa                                                                        <  2:1190740/140‑1 (MQ=255)                                            gTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                                                                       <  2:172140/140‑1 (MQ=255)                                              gggTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc                                                                     <  2:628765/140‑1 (MQ=255)                                              gggTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc                                                                     <  2:929904/140‑1 (MQ=255)                                                gttgttGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCg                                                                   <  1:729462/140‑1 (MQ=255)                                                   gttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGaa                                                                <  1:346/140‑1 (MQ=255)                                                    ttgGCGTACGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAAc                                                               >  2:1103005/1‑140 (MQ=255)                                                                   ttACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCcggcgg                                                <  1:501633/140‑1 (MQ=255)                                                                     aCCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGa                                              <  1:287748/140‑1 (MQ=255)                                                                     aCCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGa                                              <  2:1060823/140‑1 (MQ=255)                                                                            ccccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGa                                       >  2:493420/1‑140 (MQ=255)                                                                              ccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGaaa                                     <  2:656640/140‑1 (MQ=255)                                                                               cccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGaaaa                                    >  2:1122490/1‑140 (MQ=255)                                                                                    aGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATg                               <  1:919180/140‑1 (MQ=255)                                                                                      cgcTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCg                             >  1:1358684/1‑140 (MQ=255)                                                                                        cTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGtt                           >  1:56143/1‑140 (MQ=255)                                                                                           cAATGCCCGCCCTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACg                        <  1:1103005/140‑1 (MQ=255)                                                                                                cccgcccTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTc                   >  1:1388930/1‑140 (MQ=255)                                                                                                   gcccTTAAAGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCt                <  1:493420/140‑1 (MQ=255)                                                                                                          aaGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg         <  1:686629/140‑1 (MQ=255)                                                                                                          aaGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg         <  1:887665/140‑1 (MQ=255)                                                                                                          aaGTGGGGGCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg         <  1:976564/140‑1 (MQ=255)                                                                                                              gggggCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGa     <  2:1388930/140‑1 (MQ=255)                                                                                                              gggggCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGa     <  2:842924/140‑1 (MQ=255)                                                                                                               ggggCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGAt    >  1:146543/1‑140 (MQ=255)                                                                                                               ggggCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGAt    >  1:1115322/1‑140 (MQ=255)                                                                                                                 ggCACTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGc  >  1:598521/1‑140 (MQ=255)                                                                                                                     cTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg         >  1:72664/1‑129 (MQ=255)                                                                                                                     cTCCCATAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg         <  2:72664/129‑1 (MQ=255)                                                                                                                       cccaTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGa     <  1:1335497/131‑1 (MQ=255)                                                                                                                          |                                                                                                                                   CAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGC  >  NC_000913/2165493‑2165744

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