Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,376,341	+G	coding (436/498 nt)	sspB ←	ClpXP protease specificity enhancing factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,376,339	1	.	G	100.0%	159.9 / NA	42	P146P (CCT→CCC)	sspB	ClpXP protease specificity enhancing factor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (25/17);  total (25/17)

GGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑GGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCC  >  NC_000913/3376192‑3376482                                                                                                                                                     |                                                                                                                                                ggTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑ggg                                                                                                                                              <  1:433916/151‑2 (MQ=255)       gtggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAgt                                                                                                                                         >  1:133663/1‑151 (MQ=255)               tatGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTcatcat                                                                                                                                 >  2:721315/1‑151 (MQ=255)                atGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCatcatc                                                                                                                                >  1:462999/1‑151 (MQ=255)                  gATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCg                                                                                                                               <  2:365970/150‑1 (MQ=255)                     tGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGa                                                                                                                            <  1:354761/150‑1 (MQ=255)                      gCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATc                                                                                                                          >  1:343529/1‑151 (MQ=255)                        cacaATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCt                                                                                                                         >  2:445194/1‑150 (MQ=255)                           aaTCACGCTATTTTATTAAGTCATCAAAAAGACAAAACAGGCCGCCTCGCCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCCCCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGCTGAGTGTCATCATCGTGATCTGGCt                                                                                                                     >  1:711338/1‑151 (MQ=255)                             tCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTg                                                                                                                   >  1:691233/1‑151 (MQ=255)                                cGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCg                                                                                                                >  2:476535/1‑151 (MQ=255)                                 gCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGc                                                                                                               >  2:592682/1‑151 (MQ=255)                                       ttattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                                                                                          >  1:110886/1‑150 (MQ=255)                                       ttattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                                                                                          >  2:402145/1‑150 (MQ=255)                                           taAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATaa                                                                                                      <  2:456220/150‑1 (MQ=255)                                            aaGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAAcc                                                                                                    <  1:530085/151‑1 (MQ=255)                                                cATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGAc                                                                                                 >  1:10978/1‑150 (MQ=255)                                                cATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGAc                                                                                                 <  2:630393/150‑1 (MQ=255)                                                cATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGAc                                                                                                 >  1:494280/1‑150 (MQ=255)                                                     aaaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACg                                                                                           >  1:190467/1‑151 (MQ=255)                                                          aCAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTc                                                                                      >  1:536414/1‑151 (MQ=255)                                                                cAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGt                                                                                 >  2:120915/1‑150 (MQ=255)                                                                 aGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCt                                                                               <  1:526589/151‑1 (MQ=255)                                                                  ggCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTg                                                                              <  1:368021/151‑1 (MQ=255)                                                                     cgccTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCg                                                                           >  1:530384/1‑151 (MQ=255)                                                                        cTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATg                                                                        >  1:168543/1‑151 (MQ=255)                                                                             ccTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCctct                                                                   <  2:457032/151‑1 (MQ=255)                                                                                 ttttGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTtcat                                                               >  2:164854/1‑151 (MQ=255)                                                                                   ttGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTtcatc                                                              >  1:632062/1‑150 (MQ=255)                                                                                   ttGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTtcatc                                                              <  2:575863/150‑1 (MQ=255)                                                                                      tATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCAtcatt                                                           <  2:438883/150‑1 (MQ=255)                                                                                             tCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGc                                                   >  2:182521/1‑151 (MQ=255)                                                                                                 aaCGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGt                                               <  2:498312/151‑1 (MQ=255)                                                                                                      gTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCtt                                          <  2:411045/151‑1 (MQ=255)                                                                                                        aaTGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTc                                         >  2:419840/1‑150 (MQ=255)                                                                                                        aaTGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTc                                         >  2:510050/1‑150 (MQ=255)                                                                                                                       ccGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTc                          <  1:646186/150‑1 (MQ=255)                                                                                                                          cgTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGc                      >  2:667724/1‑151 (MQ=255)                                                                                                                             ggtggCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTc                    <  2:169246/150‑1 (MQ=255)                                                                                                                             ggtggCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTc                    <  2:416029/150‑1 (MQ=255)                                                                                                                                     cGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGt           <  1:164821/151‑1 (MQ=255)                                                                                                                                      ggAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTg          <  2:709471/151‑1 (MQ=255)                                                                                                                                              aTCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGcc  >  1:412936/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                GGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑GGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCC  >  NC_000913/3376192‑3376482

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