Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,659,365	T→G	K35Q (AAA→CAA)	hscB ←	co‑chaperone for [Fe‑S] cluster biosynthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,659,365	0	T	G	100.0%	204.0 / NA	53	K35Q (AAA→CAA)	hscB	co‑chaperone for [Fe‑S] cluster biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (24/29);  total (24/29)

GGCGAGATCAAAGCCGTGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGAT  >  NC_000913/2659228‑2659497                                                                                                                                          |                                                                                                                                     ggCGAGATCAAAGACGTGCAAAGAAAGCAAATATTCCGCGAGCATTAACGGATCACGCCGCCTTTGCCAGGCCTGGCTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTCCCGCCTGGATTCCGCTGGCGAATTGAt                                                                                                                                    <  2:30002/140‑1 (MQ=255)     agaTCAAAGCCGTGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAgg                                                                                                                                <  1:575838/140‑1 (MQ=255)      gaTCAAAGCCGTGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGga                                                                                                                               <  1:828210/140‑1 (MQ=255)      gaTCAAAGCCGTGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGga                                                                                                                               <  1:1242683/140‑1 (MQ=255)      gaTCAAAGCCGTGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGga                                                                                                                               <  1:1270481/140‑1 (MQ=255)              ccGTGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTg                                                                                                                       <  1:52965/140‑1 (MQ=255)                 tGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACg                                                                                                                    <  2:432006/140‑1 (MQ=255)                 tGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACg                                                                                                                    >  1:827351/1‑140 (MQ=255)                 tGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACg                                                                                                                    >  1:1050444/1‑140 (MQ=255)                          agCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATc                                                                                                           >  2:938215/1‑140 (MQ=255)                           gCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATcc                                                                                                          >  1:942419/1‑140 (MQ=255)                                    ccGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGc                                                                                                 >  1:746635/1‑140 (MQ=255)                                         cgcATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCt                                                                                            <  1:1141054/140‑1 (MQ=255)                                         cgcATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCt                                                                                            <  1:827925/140‑1 (MQ=255)                                         cgcATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCt                                                                                            <  1:582181/140‑1 (MQ=255)                                            aTTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAg                                                                                         >  2:1008487/1‑140 (MQ=255)                                            aTTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAg                                                                                         >  2:22668/1‑140 (MQ=255)                                            aTTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAg                                                                                         <  2:746635/140‑1 (MQ=255)                                            aTTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAg                                                                                         >  1:277178/1‑140 (MQ=255)                                                aCGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGcc                                                                                     <  2:923554/140‑1 (MQ=255)                                                 cGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCt                                                                                    <  2:277178/140‑1 (MQ=255)                                                   gATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTgg                                                                                  <  2:628578/140‑1 (MQ=255)                                                   gATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTgg                                                                                  <  1:22668/140‑1 (MQ=255)                                                   gATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTgg                                                                                  <  1:1008487/140‑1 (MQ=255)                                                   gATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTgg                                                                                  <  1:938215/140‑1 (MQ=255)                                                      gACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTa                                                                               >  1:249311/1‑140 (MQ=255)                                                        cGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATc                                                                             >  2:207927/1‑140 (MQ=255)                                                        cGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATc                                                                             >  2:189914/1‑140 (MQ=255)                                                                 tGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATa                                                                    <  2:249311/140‑1 (MQ=255)                                                                       gCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggc                                                              >  1:369896/1‑140 (MQ=255)                                                                         cTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGGGTTGATAGCGggcag                                                            <  2:942419/140‑1 (MQ=255)                                                                           ggTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggcaggc                                                          >  1:1113681/1‑140 (MQ=255)                                                                                 tGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCa                                                    >  1:1316792/1‑140 (MQ=255)                                                                                  ggTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCaa                                                   <  2:905877/140‑1 (MQ=255)                                                                                   gTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCaaa                                                  <  2:1113681/140‑1 (MQ=255)                                                                                      gCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggcagg                                                           >  2:1044585/1‑128 (MQ=255)                                                                                      gCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggcagg                                                           <  1:1044585/128‑1 (MQ=255)                                                                                      gCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAagag                                               >  2:707661/1‑140 (MQ=255)                                                                                        aGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAggg                                             <  2:1050444/140‑1 (MQ=255)                                                                                             gCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggcaggca                                                         <  2:265015/123‑1 (MQ=255)                                                                                             gCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggcaggca                                                         >  1:265015/1‑123 (MQ=255)                                                                                                 tACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAAt                                    <  2:1316792/140‑1 (MQ=255)                                                                                                 tACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAAt                                    <  1:189914/140‑1 (MQ=255)                                                                                                      cggcgAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATa                               >  1:567826/1‑140 (MQ=255)                                                                                                      cggcgAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATa                               >  2:544532/1‑140 (MQ=255)                                                                                                      cggcgAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCAAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATa                               >  1:143263/1‑140 (MQ=255)                                                                                                        gcgAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAAc                             >  1:831914/1‑140 (MQ=255)                                                                                                           aGTTGTTCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACAtt                          >  2:297419/1‑140 (MQ=255)                                                                                                             ttgttCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATtct                        <  2:827351/140‑1 (MQ=255)                                                                                                               gttCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCa                      <  2:831914/140‑1 (MQ=255)                                                                                                               gttCCGCCTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCa                      <  1:207927/140‑1 (MQ=255)                                                                                                                     ccTGGCTTCCGCTGGCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTg                >  1:1018357/1‑140 (MQ=255)                                                                                                                                   gCGAATTGATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGAt  >  1:1316849/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                     GGCGAGATCAAAGCCGTGCAAAGAAAGCAAATATTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGAT  >  NC_000913/2659228‑2659497

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