Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,376,403	T→C	D125G (GAC→GGC)	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,403	0	T	C	100.0%	105.4 / NA	30	D125G (GAC→GGC)	sspB	ClpXP protease specificity enhancing factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/16);  total (14/16)

CAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGACGCGGAATGCCACCA  >  NC_000913/3376255‑3376549                                                                                                                                                     |                                                                                                                                                   cAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCt                                                                                                                                                  >  1:92611/1‑151 (MQ=255)   ggCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGc                                                                                                                                                >  1:242810/1‑151 (MQ=255)                 gTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTtcat                                                                                                                                  >  2:428855/1‑151 (MQ=255)                          tACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCatg                                                                                                                          <  1:95182/150‑1 (MQ=255)                           aCTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCatga                                                                                                                         <  1:29643/150‑1 (MQ=255)                              tCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCt                                                                                                                     >  1:249728/1‑151 (MQ=255)                                     gcgTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCt                                                                                                              <  1:277442/151‑1 (MQ=255)                                         aaTGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCa                                                                                                           <  2:515259/150‑1 (MQ=255)                                             ccGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTa                                                                                                      <  2:92612/151‑1 (MQ=255)                                                     ccaccGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCtt                                                                                              >  2:411330/1‑151 (MQ=255)                                                                 ggCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAAc                                                                                   <  2:581700/150‑1 (MQ=255)                                                                          ggTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTg                                                                         <  2:295905/151‑1 (MQ=255)                                                                                             tgtCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCTCggg                                                          <  1:284200/147‑1 (MQ=255)                                                                                              gtCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTa                                                      <  1:668949/150‑1 (MQ=255)                                                                                                catcatCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGAt                                                   <  2:575670/151‑1 (MQ=255)                                                                                                      cGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCa                                              <  1:61956/150‑1 (MQ=255)                                                                                                       gTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCa                                              >  2:290977/1‑149 (MQ=255)                                                                                                           tCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACgg                                        >  2:614606/1‑151 (MQ=255)                                                                                                              ggCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGcag                                     <  1:669959/151‑1 (MQ=255)                                                                                                              ggCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGcag                                     <  2:208491/151‑1 (MQ=255)                                                                                                                  tGTCGCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGcagcc                                   <  2:25089/149‑1 (MQ=255)                                                                                                                      gCCATCAATAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGc                             >  1:633526/1‑151 (MQ=255)                                                                                                                              tAACCGACATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAAc                     >  1:16822/1‑151 (MQ=255)                                                                                                                                     cATAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGAcgc              <  2:68324/151‑1 (MQ=255)                                                                                                                                      aTAACGGTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGAcgcg             >  2:73536/1‑151 (MQ=255)                                                                                                                                           ggTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGACGCGGAATg        >  2:107063/1‑151 (MQ=255)                                                                                                                                            gTTTCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGACGCGGAATGc       >  2:623504/1‑151 (MQ=255)                                                                                                                                             tttCGTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGACGCGGAATGcc      >  1:335531/1‑151 (MQ=255)                                                                                                                                                 gTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGACGCGGAATGccacca  <  2:615371/151‑1 (MQ=255)                                                                                                                                                 gTTGCCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGACGCGGAATGccacc   >  2:477389/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATGCCTCTTCATCATTCATGATGCTGGTATCTTCATCGTAGGCAGCTTCAGGCTCAAACATCGTGCCTGCGCCATTTTCACGGGCGTAGATAGCCAGCACGGCAGCCAGCGGCACAGAAACCTGACGCGGAATGCCACCA  >  NC_000913/3376255‑3376549

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