Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,657,699	T→C	T413A (ACC→GCC)	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,699	0	T	C	100.0%	187.8 / NA	49	T413A (ACC→GCC)	hscA	DnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/24);  total (25/24)

AGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACG  >  NC_000913/2657551‑2657841                                                                                                                                                     |                                                                                                                                               aGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCAt                                                                                                                                              >  2:959913/1‑151 (MQ=255)   cgcCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTa                                                                                                                                            <  2:137141/151‑1 (MQ=255)     cccgAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTAcg                                                                                                                                          >  2:751502/4‑151 (MQ=255)     ccGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTAcg                                                                                                                                          <  2:90253/151‑1 (MQ=255)             ccACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATc                                                                                                                                   >  1:650396/1‑150 (MQ=255)                cGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACt                                                                                                                                <  2:846738/150‑1 (MQ=255)                         aaaCGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTccacca                                                                                                                      <  1:512546/151‑1 (MQ=255)                         aaaCGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTccacc                                                                                                                       >  2:442934/1‑150 (MQ=255)                          aaCGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAg                                                                                                                     <  1:136225/151‑1 (MQ=255)                           aCGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAgg                                                                                                                    >  1:195750/1‑151 (MQ=255)                                 cAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGccc                                                                                                              >  2:367841/1‑151 (MQ=255)                                  aGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGTAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCa                                                                                                             <  1:513408/151‑1 (MQ=255)                                     gAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATc                                                                                                           <  1:475582/150‑1 (MQ=255)                                           cAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCg                                                                                                     >  2:476962/1‑150 (MQ=255)                                                    aCCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCg                                                                                           <  2:521158/151‑1 (MQ=255)                                                       aGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGAc                                                                                         >  2:959977/1‑150 (MQ=255)                                                             cgcTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGa                                                                                  <  2:374027/151‑1 (MQ=255)                                                                 cACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATcac                                                                              >  2:417520/1‑151 (MQ=255)                                                                    ccTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATc                                                                           <  2:791331/151‑1 (MQ=255)                                                                    ccTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATc                                                                           <  1:556437/151‑1 (MQ=255)                                                                      tGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCa                                                                          >  1:196929/1‑150 (MQ=255)                                                                       gCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCaa                                                                         <  2:135305/150‑1 (MQ=255)                                                                       gCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCa                                                                          <  2:990470/149‑1 (MQ=255)                                                                         aTTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCa                                                                      <  2:283912/151‑1 (MQ=255)                                                                           tACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAAc                                                                    >  2:85035/1‑151 (MQ=255)                                                                                ggATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAt                                                               >  1:924477/1‑151 (MQ=255)                                                                                 gATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                                              <  2:541137/151‑1 (MQ=255)                                                                                   tAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTc                                                            >  2:384827/1‑151 (MQ=255)                                                                                    aGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTc                                                            <  1:587774/150‑1 (MQ=255)                                                                                     gACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGc                                                          <  2:3070/151‑1 (MQ=255)                                                                                              gTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCt                                                 <  2:99881/151‑1 (MQ=255)                                                                                               tCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCtt                                                >  2:73186/1‑151 (MQ=255)                                                                                                cTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTg                                               >  1:128079/1‑151 (MQ=255)                                                                                                 tGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGt                                              >  1:64924/1‑151 (MQ=255)                                                                                                     cATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAcc                                          >  1:545100/1‑151 (MQ=255)                                                                                                       tCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAccc                                         >  2:508710/1‑150 (MQ=255)                                                                                                          ttAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACccacc                                      >  2:742137/1‑150 (MQ=255)                                                                                                           tAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACccacca                                     <  1:520269/150‑1 (MQ=255)                                                                                                            aaaGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAg                                    <  1:705841/150‑1 (MQ=255)                                                                                                             aaGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGaa                                  <  2:861298/151‑1 (MQ=255)                                                                                                                  ggtgAAATCCTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATcc                             >  1:705325/1‑151 (MQ=255)                                                                                                                          ccTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGaa                      <  1:367884/150‑1 (MQ=255)                                                                                                                          ccTGAGCGCGGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAAt                     >  2:401647/1‑151 (MQ=255)                                                                                                                                 gcgGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCcgcg              >  2:735971/1‑151 (MQ=255)                                                                                                                                  cgGGCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCgcgc             >  2:202918/1‑151 (MQ=255)                                                                                                                                     gCCACCGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCaa          <  1:940276/151‑1 (MQ=255)                                                                                                                                          cGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGcg     <  2:786396/151‑1 (MQ=255)                                                                                                                                          cGGAATAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGcg     >  1:362806/1‑151 (MQ=255)                                                                                                                                             aaTAGTGGCATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGcgacg  >  2:105426/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                               AGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGCGCGGGCCACCGGAATAGTGGTATTACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACG  >  NC_000913/2657551‑2657841

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