Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	34,876	A→G	C166R (TGT→CGT)	caiE ←	stimulator of CaiD and CaiB enzyme activities

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	34,876	0	A	G	100.0%	108.4 / NA	31	C166R (TGT→CGT)	caiE	stimulator of CaiD and CaiB enzyme activities
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/14);  total (17/14)

CATCTTATCCAGCAACCAGGTCGCATCCGGCAAGATCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACAGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGCCCGCTTTGACAAAGCTC  >  NC_000913/34746‑35018                                                                                                                                   |                                                                                                                                               cATCTTATCCAGCAACCAGGTCGCATCCGGCAAGATCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGAt                                                                                                                            <  2:51943/151‑1 (MQ=255)         ccAGCAACCAGGTCGCATCCGGCAAGATCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTg                                                                                                                    <  2:133725/151‑1 (MQ=255)              aaCCAGGTCGCATCCGGCAAGATCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTAtt                                                                                                               >  1:61026/1‑151 (MQ=255)                   ggTCGCATCCGGCAAGATCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAAcc                                                                                                          <  1:232523/151‑1 (MQ=255)                        cATCCGGCAAGATCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTc                                                                                                     >  2:214319/1‑151 (MQ=255)                                   aTCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACtc                                                                                          >  2:612771/1‑151 (MQ=255)                                     cACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACtcg                                                                                         <  1:605441/150‑1 (MQ=255)                                          tttAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCgtcgtc                                                                                    >  1:185580/1‑150 (MQ=255)                                              ggCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTaa                                                                               <  1:526951/151‑1 (MQ=255)                                                cGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAAc                                                                              >  1:621359/1‑150 (MQ=255)                                                  tCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACAtt                                                                           >  1:457510/1‑151 (MQ=255)                                                   cacaTCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACAtt                                                                           <  1:283203/150‑1 (MQ=255)                                                   cacaTCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTg                                                                          <  1:61172/151‑1 (MQ=255)                                                      aTCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGt                                                                       >  2:183290/1‑151 (MQ=255)                                                                cccTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCggg                                                             >  1:247956/1‑151 (MQ=255)                                                                   tGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTg                                                          >  1:215435/1‑151 (MQ=255)                                                                      aaaCGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCg                                                       >  1:210915/1‑151 (MQ=255)                                                                           ggggCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTAccc                                                  <  1:52661/151‑1 (MQ=255)                                                                            gggCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTAccc                                                  >  1:563403/1‑150 (MQ=255)                                                                                 aTTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATcaac                                            <  2:289281/151‑1 (MQ=255)                                                                                    ttCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCaacaac                                         >  2:317750/1‑151 (MQ=255)                                                                                               gCCTCAGTGGCTGCGTTTCATGTAACGATACATGGCGGCGCCCGACAAGATCCTGATGCTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTtctc                              <  2:422172/151‑1 (MQ=255)                                                                                                     gTGGCTGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCcgcg                        <  2:621087/151‑1 (MQ=255)                                                                                                          tGCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGc                   >  1:386784/1‑151 (MQ=255)                                                                                                           gCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGc                   <  2:112316/150‑1 (MQ=255)                                                                                                           gCGTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGc                   <  2:481270/150‑1 (MQ=255)                                                                                                             gTTTCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGCCCg                >  1:129255/1‑151 (MQ=255)                                                                                                                tCATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGCCCGCtt             >  2:620578/1‑151 (MQ=255)                                                                                                                 cATGTAACGATACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGCCCGCtt             <  1:612788/150‑1 (MQ=255)                                                                                                                          aTACATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGCCCGCTTTGACAAAGCt   >  2:165420/1‑151 (MQ=255)                                                                                                                            aCATGACGGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGCCCGCTTTGACAAAGCTc  >  1:294960/1‑150 (MQ=255)                                                                                                                                   |                                                                                                                                               CATCTTATCCAGCAACCAGGTCGCATCCGGCAAGATCACCGTTTAGGCGTCACATCCGTCGTCCCCTGCAAACGGGGGCGATTTTCCTCCATTTGCCTCAGTGGCTGCGTTTCATGTAACGATACATGACAGCGCCCGACAAGATCCTGATACTCTTTGGTATTCAACCGTTTCCAGTGTAACTCGTCGTCACTAACATTGCGTACAGCGCGGGCTGGCGTACCCATCAACAACTGGCGTTTCTCGCCGCGAAAGCCCGCTTTGACAAAGCTC  >  NC_000913/34746‑35018

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