Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,558,827	T→C	intergenic (‑148/‑31)	eutB ← / → intZ	ethanolamine ammonia‑lyase subunit alpha/CPZ‑55 prophage; putative phage integrase IntZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,558,827	0	T	C	100.0%	221.6 / NA	61	intergenic (‑148/‑31)	eutB/intZ	ethanolamine ammonia‑lyase subunit alpha/CPZ‑55 prophage; putative phage integrase IntZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/33);  total (28/33)

TCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTG  >  NC_000913/2558701‑2558959                                                                                                                               |                                                                                                                                     tCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCg                                                                                                                            >  1:326731/1‑137 (MQ=255) tCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCg                                                                                                                            <  2:326731/137‑1 (MQ=255)   gatgatTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTaaa                                                                                                                       >  2:86267/1‑140 (MQ=255)                  tCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGtt                                                                                                        <  1:44130/140‑1 (MQ=255)                  tCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGtt                                                                                                        <  1:813589/140‑1 (MQ=255)                    acaGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGtttt                                                                                                      <  1:1198086/140‑1 (MQ=255)                              aaCAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGaca                                                                                            <  2:258033/140‑1 (MQ=255)                                cAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGAcacc                                                                                          <  2:145566/140‑1 (MQ=255)                                  aTCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGAcaccac                                                                                        >  1:1213101/1‑140 (MQ=255)                                  aTCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGAcaccac                                                                                        >  2:134271/1‑140 (MQ=255)                                    cACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGaca                                                                                            <  2:469288/134‑1 (MQ=255)                                    cACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGaca                                                                                            >  1:469288/1‑134 (MQ=255)                                     aCGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTca                                                                                     >  2:1175604/1‑140 (MQ=255)                                     aCGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTca                                                                                     <  2:972234/140‑1 (MQ=255)                                       gTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTcaca                                                                                   >  2:840864/1‑140 (MQ=255)                                           cATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaa                                                                               <  1:512224/140‑1 (MQ=255)                                           cATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaa                                                                               >  2:335478/1‑140 (MQ=255)                                           cATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaa                                                                               <  1:335478/140‑1 (MQ=255)                                           cATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaa                                                                               <  1:941846/140‑1 (MQ=255)                                                  aGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCAt                                                                        <  1:820947/140‑1 (MQ=255)                                                     ccTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGa                                                                     <  2:892393/140‑1 (MQ=255)                                                            gtgtAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGca                                                              <  1:86267/140‑1 (MQ=255)                                                             tgtAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGcac                                                             <  1:134271/140‑1 (MQ=255)                                                               tAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAg                                                           >  2:123283/1‑140 (MQ=255)                                                                  ggtggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCg                                                        <  1:689876/140‑1 (MQ=255)                                                                  ggtggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCg                                                        <  1:840864/140‑1 (MQ=255)                                                                  ggtggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCg                                                        >  2:1142191/1‑140 (MQ=255)                                                                      gtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAgg                                                    >  2:977373/1‑140 (MQ=255)                                                                       ttGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCg                                                        >  1:316953/1‑135 (MQ=255)                                                                       ttGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCg                                                        <  2:316953/135‑1 (MQ=255)                                                                       ttGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGa                                                   >  1:962917/1‑140 (MQ=255)                                                                       ttGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGa                                                   <  1:1142191/140‑1 (MQ=255)                                                                       ttGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGa                                                   >  2:689933/1‑140 (MQ=255)                                                                        tGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGAc                                                  <  2:518487/140‑1 (MQ=255)                                                                        tGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGAc                                                  >  2:460512/1‑140 (MQ=255)                                                                          ccTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACaa                                                <  2:767534/140‑1 (MQ=255)                                                                          ccTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACaa                                                <  1:689933/140‑1 (MQ=255)                                                                                   gCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATaa                                       >  2:196900/1‑140 (MQ=255)                                                                                    cccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGAc                                                  >  2:408144/1‑128 (MQ=255)                                                                                    cccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGAc                                                  <  1:408144/128‑1 (MQ=255)                                                                                    cccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATaaa                                      >  2:762085/1‑140 (MQ=255)                                                                                     ccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGAc                                                  >  2:686695/1‑127 (MQ=255)                                                                                     ccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGAc                                                  <  1:686695/127‑1 (MQ=255)                                                                                        ttATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCa                                  >  2:780928/1‑140 (MQ=255)                                                                                         tATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATaaa                                      >  1:444524/1‑135 (MQ=255)                                                                                         tATTCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATaaa                                      <  2:444524/135‑1 (MQ=255)                                                                                           ttCATCACGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCg                               <  1:460512/140‑1 (MQ=255)                                                                                             catcaCGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGAt                             >  2:593308/1‑140 (MQ=255)                                                                                               tcaCGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTc                           >  2:181405/1‑140 (MQ=255)                                                                                                caCGCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCa                          <  1:196900/140‑1 (MQ=255)                                                                                                   gCCGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTg                       >  2:585711/1‑140 (MQ=255)                                                                                                     cGAACGCGCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGAc                     >  2:249709/1‑140 (MQ=255)                                                                                                         cgcgCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCg                 >  2:359580/1‑140 (MQ=255)                                                                                                         cgcgCCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCg                 >  2:843672/1‑140 (MQ=255)                                                                                                            gcCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCgggg              <  1:843672/140‑1 (MQ=255)                                                                                                            gcCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCgggg              <  1:359580/140‑1 (MQ=255)                                                                                                            gcCGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCgggg              <  2:962917/140‑1 (MQ=255)                                                                                                            gcCGGATACATGATCGGGGCTATCCAGTAGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCgggg              <  1:908920/140‑1 (MQ=255)                                                                                                              cGGATACATGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCt            <  1:388365/140‑1 (MQ=255)                                                                                                                      tGATCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGc    >  1:913959/1‑140 (MQ=255)                                                                                                                        aTCGGGGCTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTg  <  1:762085/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                     TCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTG  >  NC_000913/2558701‑2558959

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