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%	140.7 / NA	37	P146P (CCT→CCC)	sspB	ClpXP protease specificity enhancing factor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (18/19);  total (18/19)

CGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑GGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTG  >  NC_000913/3376191‑3376478                                                                                                                                                      |                                                                                                                                            cGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑gg                                                                                                                                           <  1:693131/151‑1 (MQ=255) cGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑g                                                                                                                                            <  1:483874/150‑1 (MQ=255)        gtggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAgt                                                                                                                                     >  2:612299/1‑151 (MQ=255)        gtggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAg                                                                                                                                      >  1:563665/1‑150 (MQ=255)          ggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAgtg                                                                                                                                    <  1:684746/150‑1 (MQ=255)           gtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGtc                                                                                                                                  <  2:813302/151‑1 (MQ=255)                 atGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTcatcat                                                                                                                             <  1:891494/150‑1 (MQ=255)                  tGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCg                                                                                                                           >  1:494499/1‑151 (MQ=255)                   gATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCg                                                                                                                           <  1:420293/150‑1 (MQ=255)                    aTTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTg                                                                                                                         >  1:386751/1‑151 (MQ=255)                            aaTCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCt                                                                                                                 <  2:58588/151‑1 (MQ=255)                                     aTTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGcc                                                                                                          <  2:529035/149‑1 (MQ=255)                                        ttattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                                                                                      <  1:202478/150‑1 (MQ=255)                                              aGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCg                                                                                               <  2:381267/151‑1 (MQ=255)                                                tCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGAc                                                                                             >  2:897258/1‑151 (MQ=255)                                                 cATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACa                                                                                            <  2:76641/151‑1 (MQ=255)                                                       aaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACgg                                                                                      >  2:190577/1‑151 (MQ=255)                                                       aaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACgg                                                                                      >  1:382285/1‑151 (MQ=255)                                                          gACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGttt                                                                                   >  1:534218/1‑151 (MQ=255)                                                           aCAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTc                                                                                  <  2:77152/151‑1 (MQ=255)                                                               aaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTg                                                                              <  1:881754/151‑1 (MQ=255)                                                                    gccgccTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGc                                                                         <  1:556595/151‑1 (MQ=255)                                                                            ggCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCct                                                                 <  1:471541/151‑1 (MQ=255)                                                                               cTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCtt                                                              >  1:842685/1‑151 (MQ=255)                                                                                    ttGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTtcatc                                                          >  2:38639/1‑150 (MQ=255)                                                                                          tACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCAtcattcat                                                    <  1:531314/150‑1 (MQ=255)                                                                                            cTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCatgat                                                 >  2:579389/1‑151 (MQ=255)                                                                                                acaACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTg                                             >  2:711309/1‑151 (MQ=255)                                                                                                      cgTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCt                                       <  1:645585/151‑1 (MQ=255)                                                                                                            gCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATc                                  <  2:503450/150‑1 (MQ=255)                                                                                                             ccGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGt                                <  1:268754/151‑1 (MQ=255)                                                                                                              cGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTa                               <  1:379079/151‑1 (MQ=255)                                                                                                                   gACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGc                            >  2:367194/1‑149 (MQ=255)                                                                                                                       accGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTc                      >  2:114242/1‑151 (MQ=255)                                                                                                                           cgTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGc                  >  1:227288/1‑151 (MQ=255)                                                                                                                            gtggtgGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCt                 >  2:169571/1‑151 (MQ=255)                                                                                                                                       ggAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTg      >  2:487486/1‑151 (MQ=255)                                                                                                                                          ggTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCt   >  2:516603/1‑151 (MQ=255)                                                                                                                                           gTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTg  <  2:649944/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                            CGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑GGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTG  >  NC_000913/3376191‑3376478

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