Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,192,212	A→G	D140D (GAT→GAC)	hflD ←	putative lysogenization regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,192,212	0	A	G	100.0%	197.3 / NA	52	D140D (GAT→GAC)	hflD	putative lysogenization regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/30);  total (22/30)

GACGTCCGCCGCCGACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCAC  >  NC_000913/1192067‑1192352                                                                                                                                                  |                                                                                                                                             gACGTCCGCCGCCGACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAAc                                                                                                                                         >  1:408041/1‑151 (MQ=255)   cGTCCGCCGCCGACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACat                                                                                                                                       <  1:311839/151‑1 (MQ=255)       cgccgccgACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACat                                                                                                                                       >  1:536445/1‑147 (MQ=255)        gccgccgACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACatagata                                                                                                                                  <  1:827243/151‑1 (MQ=255)           gccgACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATagca                                                                                                                               <  1:344392/151‑1 (MQ=255)            ccgACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATagcag                                                                                                                              <  1:14416/151‑1 (MQ=255)             cgACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATagcagc                                                                                                                             >  1:381412/1‑151 (MQ=255)                ccTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCAt                                                                                                                          >  2:84618/1‑151 (MQ=255)                   ggTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTTACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATcg                                                                                                                        <  2:509397/150‑1 (MQ=255)                       ccAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCtc                                                                                                                   <  2:852545/151‑1 (MQ=255)                       ccAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCtc                                                                                                                   >  1:963432/1‑151 (MQ=255)                         agagCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCtcat                                                                                                                 <  1:144178/151‑1 (MQ=255)                         agagCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCtcat                                                                                                                 >  1:424672/1‑151 (MQ=255)                            gCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAg                                                                                                              <  1:848653/151‑1 (MQ=255)                               cggcggCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAGCATAGATAGCAGCCATCGCGCTCATCAGc                                                                                                             <  2:259222/149‑1 (MQ=255)                               cggcggCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGt                                                                                                           >  2:492953/1‑151 (MQ=255)                                    gcgcGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCgg                                                                                                      <  1:579/151‑1 (MQ=255)                                       cgAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACt                                                                                                   <  2:968282/151‑1 (MQ=255)                                        gAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTg                                                                                                  >  1:450784/1‑151 (MQ=255)                                          aTGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTa                                                                                                <  1:346921/151‑1 (MQ=255)                                            gcctgccAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTaa                                                                                               <  2:275591/150‑1 (MQ=255)                                             cctgccAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAAt                                                                                             >  1:398952/1‑151 (MQ=255)                                             cctgccAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAAt                                                                                             <  2:807087/151‑1 (MQ=255)                                              ctgccAGCAGGGATGCGGGAACTTTCTCCTGCACGTGCGGGCTTTGCAGTACAGCAAGGGAACCTGTGCCCTGAATGCGCGGGCCAAACGGGCTAATAACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGAATGTAAATc                                                                                            >  2:520663/1‑151 (MQ=255)                                                    gcagGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAgtg                                                                                      <  1:269634/151‑1 (MQ=255)                                                            gcgcGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCt                                                                              <  1:16940/151‑1 (MQ=255)                                                             cgcgAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTg                                                                             >  1:753310/1‑151 (MQ=255)                                                             cgcgAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTg                                                                             <  1:422921/151‑1 (MQ=255)                                                                  aCTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGtt                                                                        <  2:152773/151‑1 (MQ=255)                                                                      tCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCCCGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCCAAGTGTTCGAGCTGGCGTTGCa                                                                     <  2:194021/150‑1 (MQ=255)                                                                             cACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTg                                                             >  2:654534/1‑151 (MQ=255)                                                                              aCTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTg                                                             >  2:693625/1‑150 (MQ=255)                                                                              aCTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGa                                                            <  1:692233/151‑1 (MQ=255)                                                                              aCTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGa                                                            <  1:791825/151‑1 (MQ=255)                                                                                           gCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCa                                                <  1:831364/150‑1 (MQ=255)                                                                                            cAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCaga                                              <  1:352912/151‑1 (MQ=255)                                                                                             aGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCag                                               >  2:54872/1‑149 (MQ=255)                                                                                                aCAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAgtgt                                          >  1:810420/1‑151 (MQ=255)                                                                                                       gggAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAgcgc                                   <  2:155988/151‑1 (MQ=255)                                                                                                             cGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCttt                             <  2:363731/151‑1 (MQ=255)                                                                                                                 gACCTGAATGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCt                         >  2:418992/1‑151 (MQ=255)                                                                                                                       aaTGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTgaggag                   <  1:419003/151‑1 (MQ=255)                                                                                                                        aTGCGCGGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTgagg                     >  2:152405/1‑148 (MQ=255)                                                                                                                             cgGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTg             <  1:605562/151‑1 (MQ=255)                                                                                                                             cgGGCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTg             <  2:923506/151‑1 (MQ=255)                                                                                                                                gCCAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTgcgc          <  2:587170/151‑1 (MQ=255)                                                                                                                                 ccAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCt         >  2:676181/1‑151 (MQ=255)                                                                                                                                 ccAAGCGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCt         <  2:1037260/151‑1 (MQ=255)                                                                                                                                      cGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGc    >  2:395351/1‑151 (MQ=255)                                                                                                                                      cGGGCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGc    >  2:313987/1‑151 (MQ=255)                                                                                                                                       gggCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGc    >  2:995817/1‑150 (MQ=255)                                                                                                                                        ggCTAATCACGTCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCAc  >  1:915419/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                             GACGTCCGCCGCCGACCTGGTGCCAGAGCACGGCGGCGCGAATGCCTGCCAGCAGGGTTGCGCGAACTTTCGCCTGCACTTGTGGGCTTTGCAGTACAGCAGGGGAACCGGTGACCTGAATGCGCGGGCCAAGCGGGCTAATCACATCAACATAGATAGCAGCCATCGCGCTCATCAGCGTTTCGGACTGTAAATCGAAGTGTTCGAGCTGGCGTTGCAGGCCGTTGATCCGGTTGCCCAGAGTGTCGAGCGCGCCTTTCGCTGAGGAGAGTTTGCGCTCAAGCAC  >  NC_000913/1192067‑1192352

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