Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,042,112	T→C	intergenic (+87/‑256)	zinT → / → yodB	zinc and cadmium binding protein, periplasmic/cytochrome b561 homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,042,112	0	T	C	100.0%	200.6 / NA	52	intergenic (+87/‑256)	zinT/yodB	zinc and cadmium binding protein, periplasmic/cytochrome b561 homolog
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/33);  total (19/33)

AACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGT  >  NC_000913/2041965‑2042254                                                                                                                                                    |                                                                                                                                               aaCGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAagccagc                                                                                                                                             <  1:125877/151‑1 (MQ=255)   cGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAagccagct                                                                                                                                            >  2:426609/1‑150 (MQ=255)      attatCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTAc                                                                                                                                          >  2:325455/1‑149 (MQ=255)        tatCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGa                                                                                                                                       <  1:639527/150‑1 (MQ=255)           ccATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGacac                                                                                                                                    <  2:659595/150‑1 (MQ=255)                 cAGTTGAGTAGCGAAGAAGTGTTCGAGGAAATGATGTCTAATTGAGTATCCTCATGATAACGCCTCGATGCCCCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCtcat                                                                                                                              <  2:322733/150‑1 (MQ=255)                 cAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCtcatc                                                                                                                             <  1:426105/151‑1 (MQ=255)                          aGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACg                                                                                                                    >  1:877648/1‑151 (MQ=255)                            cGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACgtg                                                                                                                  >  1:283064/1‑151 (MQ=255)                              aagaagTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACgtgt                                                                                                                 <  1:669268/150‑1 (MQ=255)                                  agTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGACAGCTACAGACACGCTCATCTCTATGACGTGTCTg                                                                                                              <  1:256198/149‑1 (MQ=255)                                     ggTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTatat                                                                                                         <  1:671995/151‑1 (MQ=255)                                      gTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATAtt                                                                                                        <  1:167921/151‑1 (MQ=255)                                      gTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATAtt                                                                                                        <  1:792996/151‑1 (MQ=255)                                         gAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATtaat                                                                                                     >  1:499850/1‑151 (MQ=255)                                          aGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATtaat                                                                                                     >  1:424611/1‑150 (MQ=255)                                           ggAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATtaataa                                                                                                   <  1:118327/151‑1 (MQ=255)                                                tgatgTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTa                                                                                               >  2:580438/1‑150 (MQ=255)                                                 gatgTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAAc                                                                                             >  1:244035/1‑151 (MQ=255)                                                        cATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCa                                                                                        <  2:522881/149‑1 (MQ=255)                                                        cATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCAtt                                                                                      <  1:504954/151‑1 (MQ=255)                                                        cATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCAtt                                                                                      >  2:447437/1‑151 (MQ=255)                                                        cATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCAt                                                                                       <  1:787487/150‑1 (MQ=255)                                                          ttGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCCCAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTg                                                                                     <  1:817982/150‑1 (MQ=255)                                                            gAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAg                                                                                   <  1:580493/150‑1 (MQ=255)                                                             aGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGt                                                                                  <  2:634651/150‑1 (MQ=255)                                                               tATTCTCATGATAACGCCTCGCTGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTa                                                                                >  1:809829/1‑150 (MQ=255)                                                                aTTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAAc                                                                              <  2:907759/151‑1 (MQ=255)                                                                 ttCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAAcc                                                                             >  1:193052/1‑151 (MQ=255)                                                                     cATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAAt                                                                          <  1:1049710/150‑1 (MQ=255)                                                                      aTGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATaa                                                                        <  2:49892/151‑1 (MQ=255)                                                                           aaCGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGAt                                                                   <  2:998897/151‑1 (MQ=255)                                                                                 tCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCAt                                                              >  2:676882/1‑150 (MQ=255)                                                                                        cccTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACTCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAAGAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAAt                                                       <  2:698421/148‑1 (MQ=255)                                                                                         gCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATa                                                      >  2:895713/1‑150 (MQ=255)                                                                                            ttAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACg                                                    <  1:151573/149‑1 (MQ=255)                                                                                             tAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGc                                                  <  2:655698/150‑1 (MQ=255)                                                                                              aGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGcc                                                 >  1:161185/1‑150 (MQ=255)                                                                                                 aaGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGAAACGCTCATCACTATGACGTGTCTTTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCaa                                               <  2:611482/149‑1 (MQ=255)                                                                                                    ttATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACaga                                           <  2:526640/150‑1 (MQ=255)                                                                                                           aaCTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGAt                                    <  1:695589/150‑1 (MQ=255)                                                                                                            aCTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATAc                                  <  1:81105/151‑1 (MQ=255)                                                                                                               gCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCt                                <  2:845929/150‑1 (MQ=255)                                                                                                                ccACTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCtgt                              >  2:90104/1‑151 (MQ=255)                                                                                                                  aCTGGTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGc                            >  1:666565/1‑151 (MQ=255)                                                                                                                      gTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCAc                        <  1:362840/151‑1 (MQ=255)                                                                                                                      gTCCATCCACAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCAc                        >  2:723965/1‑151 (MQ=255)                                                                                                                       tccatccaCAAACGCCACTGAACGAAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATTTTAATAAGCAAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATAcgtcagcaca                         >  2:658846/1‑140 (MQ=255)                                                                                                                          atccaCAAACGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCa                    <  1:966615/151‑1 (MQ=255)                                                                                                                                aaaCGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCtt              >  1:1026751/1‑151 (MQ=255)                                                                                                                                   cGCCACTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTtata           <  2:506727/151‑1 (MQ=255)                                                                                                                                       aCTGAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTgg       <  2:41068/151‑1 (MQ=255)                                                                                                                                          gAACGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCt     <  1:494573/150‑1 (MQ=255)                                                                                                                                            aCGCAAGCCAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGt  >  2:694866/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                               AACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGT  >  NC_000913/2041965‑2042254

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