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

GTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑GGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCC  >  NC_000913/3376189‑3376482                                                                                                                                                        |                                                                                                                                                gTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCa                                                                                                                                                  <  2:596231/151‑1 (MQ=255)   cGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑gg                                                                                                                                               >  1:740008/1‑151 (MQ=255)    ggTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑ggg                                                                                                                                              >  1:806575/1‑150 (MQ=255)            ggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAgtgt                                                                                                                                       >  2:804944/1‑151 (MQ=255)                   atGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCatcatc                                                                                                                                >  2:812894/1‑151 (MQ=255)                     gATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCg                                                                                                                               >  2:495923/1‑150 (MQ=255)                        tGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGa                                                                                                                            >  2:661865/1‑150 (MQ=255)                         gCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGa                                                                                                                            <  1:721911/149‑1 (MQ=255)                         gCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGAt                                                                                                                           <  2:636530/150‑1 (MQ=255)                                 cACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGt                                                                                                                  <  2:610996/151‑1 (MQ=255)                                       aTTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCAt                                                                                                            >  2:374326/1‑151 (MQ=255)                                         tttATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                                                                                           >  1:454667/1‑150 (MQ=255)                                          ttattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                                                                                          >  2:650899/1‑150 (MQ=255)                                          ttattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                                                                                          <  1:495887/150‑1 (MQ=255)                                            attaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                                                                                        >  1:790180/1‑150 (MQ=255)                                            attaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATa                                                                                                       <  2:9076/151‑1 (MQ=255)                                                    aTTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACAt                                                                                               >  1:638911/1‑151 (MQ=255)                                                       aaaaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAAc                                                                                            >  1:687152/1‑151 (MQ=255)                                                       aaaaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAAc                                                                                            <  1:56905/151‑1 (MQ=255)                                                         aaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACgg                                                                                          >  2:100012/1‑151 (MQ=255)                                                         aaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACgg                                                                                          <  1:130429/151‑1 (MQ=255)                                                               aaaaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGt                                                                                    >  2:145765/1‑151 (MQ=255)                                                                  aCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGt                                                                                 >  2:362603/1‑151 (MQ=255)                                                                     ggCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTg                                                                              <  1:185817/151‑1 (MQ=255)                                                                        cgccTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGGGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGGGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCg                                                                           <  2:745814/151‑1 (MQ=255)                                                                             gggCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATg                                                                        <  1:804888/149‑1 (MQ=255)                                                                                  tGTTTTGTATTACTTCACAACGCTTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCctct                                                                   <  1:684006/149‑1 (MQ=255)                                                                                          aTTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCAtcattca                                                         <  2:122052/151‑1 (MQ=255)                                                                                                   caACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTgg                                                >  2:105040/1‑151 (MQ=255)                                                                                                   caACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTgg                                                >  2:819792/1‑151 (MQ=255)                                                                                                      cgcgTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTAt                                             <  2:11486/151‑1 (MQ=255)                                                                                                           aaTGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTc                                         <  2:370912/150‑1 (MQ=255)                                                                                                            aTGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCAt                                       >  2:207421/1‑151 (MQ=255)                                                                                                            aTGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCAt                                       <  1:339814/151‑1 (MQ=255)                                                                                                                 ggTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAg                                  <  2:843350/151‑1 (MQ=255)                                                                                                                    cGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGc                                >  1:751792/1‑150 (MQ=255)                                                                                                                    cGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGc                                <  1:455140/150‑1 (MQ=255)                                                                                                                     gACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGc                                >  1:547796/1‑149 (MQ=255)                                                                                                                          ccGCGTGGTGGCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCa                         <  1:730758/151‑1 (MQ=255)                                                                                                                                ggtggCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCa                   <  1:130209/151‑1 (MQ=255)                                                                                                                                    gCTGCGGAGGTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAAc                >  1:466912/1‑150 (MQ=255)                                                                                                                                             gTTCATCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTg      >  1:240734/1‑151 (MQ=255)                                                                                                                                                  tCGTCAGGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGcc  <  2:88058/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                GTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCA‑GGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCC  >  NC_000913/3376189‑3376482

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