Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,390,985	C→T	E175E (GAG→GAA)	rsgA ←	ribosome small subunit‑dependent GTPase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,390,985	0	C	T	100.0%	226.7 / NA	57	E175E (GAG→GAA)	rsgA	ribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/25);  total (32/25)

TCAGCAGGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAGAC  >  NC_000913/4390837‑4391125                                                                                                                                                     |                                                                                                                                             tCAGCAGGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTc                                                                                                                                            >  2:136181/1‑151 (MQ=255)     cagGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTg                                                                                                                                        <  2:185617/151‑1 (MQ=255)     cagGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTg                                                                                                                                        <  1:385704/151‑1 (MQ=255)      agGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGa                                                                                                                                       <  2:40285/151‑1 (MQ=255)        gCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACg                                                                                                                                     >  2:654283/1‑151 (MQ=255)           ggATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGaa                                                                                                                                   <  2:443587/150‑1 (MQ=255)                 gCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATg                                                                                                                            <  2:410407/151‑1 (MQ=255)                  ccGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGc                                                                                                                           >  1:491480/1‑151 (MQ=255)                         cAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCAACCCTGAGTTTGCCTGGAAACCACCATTAGCCGATAACCATTATTCCGGTAGATATCCATTGGTCCGTTGCGGAACGCCATGctt                                                                                                                         >  2:363543/1‑144 (MQ=255)                          aGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTtcgtcg                                                                                                                   <  2:659681/151‑1 (MQ=255)                            aCTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCgtcgtc                                                                                                                 >  2:675440/1‑151 (MQ=255)                             cTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTcc                                                                                                                >  2:635008/1‑151 (MQ=255)                                cccGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagc                                                                                                             <  1:306773/151‑1 (MQ=255)                                cccGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagc                                                                                                             >  2:61054/1‑151 (MQ=255)                                cccGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagc                                                                                                             >  1:475865/1‑151 (MQ=255)                                 ccGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCGGCGGTTTTAGCCCATCCTGAGCATGACTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGGCGCCcag                                                                                                              >  2:435772/1‑149 (MQ=255)                                          aTGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCt                                                                                                   <  2:647775/151‑1 (MQ=255)                                            gCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCttg                                                                                                 >  2:204196/1‑151 (MQ=255)                                               gATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCttgttg                                                                                              >  2:242287/1‑151 (MQ=255)                                                aTGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGa                                                                                             <  2:639675/151‑1 (MQ=255)                                                aTGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGa                                                                                             <  1:180738/151‑1 (MQ=255)                                                  gcgcCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAg                                                                                            <  1:437900/150‑1 (MQ=255)                                                    gcCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGcac                                                                                         >  1:446130/1‑151 (MQ=255)                                                     cccGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGcac                                                                                         >  2:730794/1‑150 (MQ=255)                                                       cGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACaat                                                                                      <  1:600765/151‑1 (MQ=255)                                                        ggTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACaata                                                                                     <  2:422983/151‑1 (MQ=255)                                                           cAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGCCTATCTTGTTGAGCACaataat                                                                                   >  1:353809/1‑150 (MQ=255)                                                                ccTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTc                                                                             <  1:749463/151‑1 (MQ=255)                                                                ccTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTc                                                                             >  1:359276/1‑151 (MQ=255)                                                                     tCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCt                                                                        >  1:669858/1‑151 (MQ=255)                                                                      ccAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCt                                                                        <  2:348965/150‑1 (MQ=255)                                                                          cGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATATCCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCa                                                                     >  2:162653/1‑149 (MQ=255)                                                                           ggTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAgg                                                                  >  1:547087/1‑151 (MQ=255)                                                                           ggTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAgg                                                                  <  2:715153/151‑1 (MQ=255)                                                                              tttAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGttt                                                               >  1:278876/1‑151 (MQ=255)                                                                                 aGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGc                                                            <  2:204826/151‑1 (MQ=255)                                                                                           gAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacc                                                    <  1:49728/149‑1 (MQ=255)                                                                                             gTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGt                                                <  1:755651/151‑1 (MQ=255)                                                                                              tATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGt                                                >  1:59520/1‑150 (MQ=255)                                                                                               aTGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTaa                                              <  2:310618/151‑1 (MQ=255)                                                                                               aTGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTa                                               >  2:371240/1‑150 (MQ=255)                                                                                                          ccATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGataat                                   <  2:266258/151‑1 (MQ=255)                                                                                                               aaTACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTg                               >  1:75403/1‑150 (MQ=255)                                                                                                                  aCGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGc                            <  2:398162/150‑1 (MQ=255)                                                                                                                   cgcgATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGa                          <  2:671690/151‑1 (MQ=255)                                                                                                                    gcgATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGAc                         >  1:123543/1‑151 (MQ=255)                                                                                                                    gcgATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGAc                         >  1:320270/1‑151 (MQ=255)                                                                                                                     cgATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGATCGACa                        <  2:740287/151‑1 (MQ=255)                                                                                                                     cgATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGAc                         >  2:679227/1‑150 (MQ=255)                                                                                                                         caCCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCACCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGCTCGATAATATTGAGCGAc                         >  2:403298/2‑146 (MQ=255)                                                                                                                         aaCCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTc                    >  2:216309/1‑151 (MQ=255)                                                                                                                         aaCCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTc                    >  2:464562/1‑151 (MQ=255)                                                                                                                         aaCCGATATTGCGGTAGATATACATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTc                    >  2:116364/1‑151 (MQ=255)                                                                                                                              atatTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGc                >  1:605078/1‑150 (MQ=255)                                                                                                                                        aGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGCCGATAATATTGAGCGACAGCTCCGGCAAAATGGCgg      >  2:440412/1‑150 (MQ=255)                                                                                                                                          atatCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGgaga   <  1:709954/151‑1 (MQ=255)                                                                                                                                           tatCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAgac  >  1:700572/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                             TCAGCAGGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAGAC  >  NC_000913/4390837‑4391125

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