Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,403,176	C→T	P380L (CCG→CTG)	hflK →	modulator for HflB protease specific for phage lambda cII repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,403,176	0	C	T	100.0%	184.6 / NA	50	P380L (CCG→CTG)	hflK	modulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/22);  total (28/22)

TGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTATCGC  >  NC_000913/4403029‑4403319                                                                                                                                                    |                                                                                                                                                tGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgctgcc                                                                                                                                              >  2:59121/1‑151 (MQ=255)    aaaaaGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGc                                                                                                                                           >  1:678874/1‑151 (MQ=255)        aGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCctct                                                                                                                                       >  1:836992/1‑151 (MQ=255)           gtTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTtc                                                                                                                                     <  2:301133/150‑1 (MQ=255)                 tAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCcac                                                                                                                               >  2:528820/1‑150 (MQ=255)                  aaCACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCcaca                                                                                                                              >  1:326467/1‑150 (MQ=255)                    cacCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAAcc                                                                                                                           >  2:435759/1‑151 (MQ=255)                      cccGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAAcc                                                                                                                           >  1:740632/1‑149 (MQ=255)                       ccGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAg                                                                                                                         >  2:141105/1‑150 (MQ=255)                       ccGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGt                                                                                                                        >  1:169674/1‑151 (MQ=255)                                gCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                                                                                               <  1:284474/151‑1 (MQ=255)                                 cTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                                                                                               >  2:2212/1‑150 (MQ=255)                                      tAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcg                                                                                                         >  2:589272/1‑151 (MQ=255)                                      tAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcg                                                                                                         >  1:273552/1‑151 (MQ=255)                                        aCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgtc                                                                                                       <  2:651656/151‑1 (MQ=255)                                            tAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGt                                                                                                                        <  2:762985/130‑1 (MQ=255)                                            tAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGt                                                                                                                        >  1:762928/1‑130 (MQ=255)                                             aaaGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcca                                                                                                  <  2:840361/151‑1 (MQ=255)                                                        ccTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAggg                                                                                        <  1:741096/150‑1 (MQ=255)                                                          tGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCg                                                                                      >  1:629537/1‑150 (MQ=255)                                                               gTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtatt                                                                                 >  1:443482/1‑150 (MQ=255)                                                                   tGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATgg                                                                             >  2:733017/1‑150 (MQ=255)                                                                   tGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATgg                                                                             >  2:187670/1‑150 (MQ=255)                                                                       gTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTTCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCTCCAGTCAGGGCGATATTATGGACCaa                                                                        <  1:477561/151‑1 (MQ=255)                                                                             ccAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAAcgccgc                                                                  <  2:796607/151‑1 (MQ=255)                                                                              cAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCcgcg                                                                 >  2:67281/1‑151 (MQ=255)                                                                                 aTGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCa                                                              <  2:118663/151‑1 (MQ=255)                                                                                          ggtggtAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGc                                                     <  2:723943/151‑1 (MQ=255)                                                                                                     cccTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTAc                                          >  1:230807/1‑151 (MQ=255)                                                                                                        tGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTAcc                                         >  2:64659/1‑149 (MQ=255)                                                                                                         gcggcgAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAg                                       <  1:827712/150‑1 (MQ=255)                                                                                                             cgAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCCACGCGCAGCGTAACGACTACCAGCGt                                    >  1:390738/1‑149 (MQ=255)                                                                                                              gAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTc                                   >  2:304548/1‑149 (MQ=255)                                                                                                              gAAGAGCCATAACGGTGCCAGCAGTCTGCTGCTTCTGCTGCCAGCCTCATCCTCCACAACCAGTGGAGCAAGCAACACATCGTACACCAGTCCGGGCGATATTAGGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAg                                 <  2:621465/151‑1 (MQ=255)                                                                                                                agagCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAggg                               <  1:485988/151‑1 (MQ=255)                                                                                                                   gCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAggggg                             >  1:407836/1‑150 (MQ=255)                                                                                                                    cGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGaa                           <  1:850917/151‑1 (MQ=255)                                                                                                                      aTAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATa                         <  2:237278/151‑1 (MQ=255)                                                                                                                           gttGCCAGCAACCTGCTGCGCCTTCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGCCGTCCTCCAGTCAGGGCGATATTATGGAGCAACGCCGCGCCAACGCGCAGCGGAACGACTTCCAGCGTCAGGGGGAAGAACGa                     <  1:725112/148‑1 (MQ=255)                                                                                                                                 aGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTaa              <  1:222452/151‑1 (MQ=255)                                                                                                                                 aGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTa               >  2:654189/1‑150 (MQ=255)                                                                                                                                   cAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGt            >  2:950/1‑151 (MQ=255)                                                                                                                                         gctgcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGGTGCGTAAGGCAGtt       >  1:35393/1‑150 (MQ=255)                                                                                                                                         gctgcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAg         <  1:210805/148‑1 (MQ=255)                                                                                                                                         gctgcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGtt       <  2:398046/150‑1 (MQ=255)                                                                                                                                         gctgcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTa      <  2:230831/151‑1 (MQ=255)                                                                                                                                         gctgcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAAAGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGtt       >  2:576495/1‑150 (MQ=255)                                                                                                                                            gcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTATc    <  2:726576/150‑1 (MQ=255)                                                                                                                                             cGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTATcgc  >  2:402278/1‑151 (MQ=255)                                                                                                                                               tCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCg                 <  2:711460/134‑1 (MQ=255)                                                                                                                                               tCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCg                 >  1:711417/1‑134 (MQ=255)                                                                                                                                                    |                                                                                                                                                TGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTATCGC  >  NC_000913/4403029‑4403319

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