Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,075,345	A→G	N226S (AAC→AGC)	yihW →	putative DNA‑binding transcriptional regulator YihW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,075,345	0	A	G	100.0%	200.0 / NA	58	N226S (AAC→AGC)	yihW	putative DNA‑binding transcriptional regulator YihW
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/33);  total (25/33)

TGAGAGCGATGGCGCGTTGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCA  >  NC_000913/4075208‑4075481                                                                                                                                          |                                                                                                                                         tGAGAGCGATGGCGCGTTGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCg                                                                                                                                        >  1:297876/1‑140 (MQ=255)   agagCGATGGCGCGTTGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTg                                                                                                                                      <  1:249170/140‑1 (MQ=255)          tGGCGCGTTGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAgg                                                                                                                               >  2:240325/1‑140 (MQ=255)          tGGCGCGTTGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAgg                                                                                                                               >  1:600875/1‑140 (MQ=255)                  tGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcg                                                                                                                       <  2:600875/140‑1 (MQ=255)                  tGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcg                                                                                                                       <  2:297876/140‑1 (MQ=255)                  tGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcg                                                                                                                       >  2:255434/1‑140 (MQ=255)                  tGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcg                                                                                                                       <  1:321124/140‑1 (MQ=255)                   gATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcgc                                                                                                                      <  2:181250/140‑1 (MQ=255)                   gATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcgc                                                                                                                      <  2:176534/140‑1 (MQ=255)                         gTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTa                                                                                                                >  1:344691/1‑140 (MQ=255)                            tGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACcg                                                                                                             >  2:277665/1‑140 (MQ=255)                              aTGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcg                                                                                                                       <  2:826969/128‑1 (MQ=255)                              aTGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTgcg                                                                                                                       >  1:826969/1‑128 (MQ=255)                               tGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACcgacg                                                                                                          <  2:51274/140‑1 (MQ=255)                                    aCGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTg                                                                                                     >  2:741245/1‑140 (MQ=255)                                        aGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTTGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTgccgc                                                                                                 >  1:801959/1‑140 (MQ=255)                                        aGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTgccgc                                                                                                 >  2:73811/1‑140 (MQ=255)                                                  gtgAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCt                                                                                       >  2:79647/1‑140 (MQ=255)                                                   tgAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTa                                                                                      >  2:656006/1‑140 (MQ=255)                                                        aCGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATc                                                                                 <  1:255434/140‑1 (MQ=255)                                                          gatgatGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCAc                                                                               <  2:66050/140‑1 (MQ=255)                                                          gatgatGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCAc                                                                               <  1:711960/140‑1 (MQ=255)                                                          gatgatGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCAc                                                                               <  1:584732/140‑1 (MQ=255)                                                          gatgatGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCAc                                                                               <  2:574839/140‑1 (MQ=255)                                                           atgatgGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACg                                                                              <  2:525205/140‑1 (MQ=255)                                                            tgatgGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGc                                                                             >  2:516384/1‑140 (MQ=255)                                                                      cgcgAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACa                                                                   >  1:194175/1‑140 (MQ=255)                                                                         gagaAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGcc                                                                <  2:194175/140‑1 (MQ=255)                                                                                atTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGa                                                         <  2:798229/140‑1 (MQ=255)                                                                                 ttCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGaa                                                        <  1:516384/140‑1 (MQ=255)                                                                                       tgGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCAtt                                                  <  1:89338/140‑1 (MQ=255)                                                                                       tgGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCAtt                                                  <  1:738645/140‑1 (MQ=255)                                                                                        ggTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCAttc                                                 <  1:741245/140‑1 (MQ=255)                                                                                        ggTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCAttc                                                 <  2:344691/140‑1 (MQ=255)                                                                                            gCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCGCGCTTACAAGACAGCCAAATTGAAATCATTCTTcc                                             <  1:240325/140‑1 (MQ=255)                                                                                            gCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTcc                                             <  2:671762/140‑1 (MQ=255)                                                                                                          gTATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGt                               >  1:27875/1‑140 (MQ=255)                                                                                                           tATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAgg                                        >  1:773769/1‑130 (MQ=255)                                                                                                           tATCATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAgg                                        <  2:773769/130‑1 (MQ=255)                                                                                                              cATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGAt                           >  1:233630/1‑140 (MQ=255)                                                                                                              cATGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGAt                           >  1:499815/1‑140 (MQ=255)                                                                                                               aTGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGAtt                          >  2:807491/1‑140 (MQ=255)                                                                                                                tGCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGAttt                         >  1:749976/1‑140 (MQ=255)                                                                                                                 gCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGAtttt                        >  2:428538/1‑140 (MQ=255)                                                                                                                 gCTTCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGAtttt                        >  1:627880/1‑140 (MQ=255)                                                                                                                   ttCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTtg                      >  1:412749/1‑140 (MQ=255)                                                                                                                    tCGGCGGCGGTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTtgt                     <  2:749976/140‑1 (MQ=255)                                                                                                                     cggcggcggTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTtgtg                    <  2:627880/140‑1 (MQ=255)                                                                                                                     cggcggcggTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTtgtg                    >  2:650841/1‑140 (MQ=255)                                                                                                                           cggTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGAcc                 >  2:632999/1‑137 (MQ=255)                                                                                                                           cggTTGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGAcc                 <  1:632999/137‑1 (MQ=255)                                                                                                                               tGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAAcc          <  1:16917/140‑1 (MQ=255)                                                                                                                               tGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAAcc          <  1:73811/140‑1 (MQ=255)                                                                                                                               tGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAAcc          <  1:289501/140‑1 (MQ=255)                                                                                                                               tGAAATTGGTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAAcc          <  2:233630/140‑1 (MQ=255)                                                                                                                                       gTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCa  <  1:656006/140‑1 (MQ=255)                                                                                                                                       gTAGCGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCa  <  2:27875/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         TGAGAGCGATGGCGCGTTGATGGAGTTTGATGTAAACGAAGCTAACGTGGTGAAAACGATGATGGCGCACGCGAGAAATATTCTGCTGGTCGCCGATCACACTAAGTATCATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCA  >  NC_000913/4075208‑4075481

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