Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,657,657	T→G	K427Q (AAA→CAA)	hscA ←	DnaK‑like molecular chaperone specific for IscU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,657,657	0	T	G	100.0%	219.6 / NA	57	K427Q (AAA→CAA)	hscA	DnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/27);  total (30/27)

GACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCC  >  NC_000913/2657511‑2657802                                                                                                                                                   |                                                                                                                                                  gACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAg                                                                                                                                               <  2:386926/151‑1 (MQ=255)      ggAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAggtggt                                                                                                                                          >  2:817604/1‑151 (MQ=255)      ggAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAggtgg                                                                                                                                           <  2:598904/150‑1 (MQ=255)       gAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGgtggtg                                                                                                                                         >  1:646790/1‑151 (MQ=255)         aCGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCGTTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGaa                                                                                                                                       <  2:800710/151‑1 (MQ=255)         aCGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGaa                                                                                                                                       <  1:228302/151‑1 (MQ=255)             cACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATc                                                                                                                                    <  2:577000/150‑1 (MQ=255)              aCGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATc                                                                                                                                    >  1:485589/1‑149 (MQ=255)                gcgAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGa                                                                                                                                <  2:245881/151‑1 (MQ=255)                 cgAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTg                                                                                                                                 >  2:816725/1‑149 (MQ=255)                   aaTATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAgcg                                                                                                                             >  2:562639/1‑151 (MQ=255)                         cgcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGcc                                                                                                                       >  2:502798/1‑151 (MQ=255)                          gcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCa                                                                                                                      <  1:662366/151‑1 (MQ=255)                          gcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCa                                                                                                                      <  2:216393/151‑1 (MQ=255)                              cGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCgg                                                                                                                  >  1:690549/1‑151 (MQ=255)                               gCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCgg                                                                                                                  >  1:277341/1‑150 (MQ=255)                                     cGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTg                                                                                                           >  1:567911/1‑151 (MQ=255)                                         aGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTAt                                                                                                       <  1:674763/151‑1 (MQ=255)                                           cgcCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTa                                                                                                     <  2:263281/151‑1 (MQ=255)                                            gcCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTAc                                                                                                    >  1:325398/1‑151 (MQ=255)                                            gcCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTAc                                                                                                    >  1:394213/1‑151 (MQ=255)                                                   tACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGCGGTATTACGCGGaa                                                                                              >  2:690881/1‑150 (MQ=255)                                                    aCCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAAt                                                                                             <  2:260797/150‑1 (MQ=255)                                                     ccACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATc                                                                                            >  2:85253/1‑150 (MQ=255)                                                     ccACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCa                                                                                           <  1:703595/151‑1 (MQ=255)                                                        cGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCGCCGGAATAGTGGTATTACGCGGAATCACtt                                                                                        <  1:677849/151‑1 (MQ=255)                                                        cGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACt                                                                                         <  1:102735/150‑1 (MQ=255)                                                               cAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTccac                                                                                 <  1:296432/151‑1 (MQ=255)                                                                  aaCGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTccacc                                                                                >  1:280296/1‑149 (MQ=255)                                                                    cgcgCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGgc                                                                            >  1:603294/1‑151 (MQ=255)                                                                               gCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGtt                                                                 <  2:364782/151‑1 (MQ=255)                                                                                  gCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCg                                                              <  2:238479/151‑1 (MQ=255)                                                                                    aGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCg                                                              >  2:174722/1‑149 (MQ=255)                                                                                     gTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAgg                                                           <  1:144154/151‑1 (MQ=255)                                                                                     gTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAgg                                                           >  1:899710/1‑151 (MQ=255)                                                                                         tGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGccc                                                        >  2:810003/1‑150 (MQ=255)                                                                                          gCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCa                                                       <  1:559642/150‑1 (MQ=255)                                                                                           cACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGc                                                     >  1:234569/1‑151 (MQ=255)                                                                                                  tCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTg                                              >  2:667841/1‑151 (MQ=255)                                                                                                    gcgcTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTggg                                            <  1:455451/151‑1 (MQ=255)                                                                                                         cACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATc                                         <  2:439942/149‑1 (MQ=255)                                                                                                          aCCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATcac                                       >  1:353116/1‑150 (MQ=255)                                                                                                           cccTGCATTACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACAt                                     >  1:814986/1‑151 (MQ=255)                                                                                                                   tACATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCaa                              <  1:406085/150‑1 (MQ=255)                                                                                                                    aCATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACa                            <  2:443005/151‑1 (MQ=255)                                                                                                                     cATGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAAc                             >  2:526336/1‑149 (MQ=255)                                                                                                                      aTGGATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGc                          >  2:890329/1‑151 (MQ=255)                                                                                                                         gATAGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                       <  1:820119/151‑1 (MQ=255)                                                                                                                            aGACATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCg                    <  2:426972/151‑1 (MQ=255)                                                                                                                               cATCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                 <  2:603313/151‑1 (MQ=255)                                                                                                                                aTCGCCGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGt                >  2:888476/1‑151 (MQ=255)                                                                                                                                     cGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGc           >  2:74876/1‑151 (MQ=255)                                                                                                                                     cGTCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTAGc           <  1:289186/151‑1 (MQ=255)                                                                                                                                       tCTGACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCtt         >  2:635962/1‑151 (MQ=255)                                                                                                                                          gACCATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGtt      >  2:178378/1‑151 (MQ=255)                                                                                                                                            ccATCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAc    >  2:780309/1‑151 (MQ=255)                                                                                                                                               tCTTGAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAccc  >  2:314700/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCC  >  NC_000913/2657511‑2657802

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