Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,521,593	C→T	*295* (TAG→TAA)	xapR ←	transcriptional activator of xapAB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,521,593	0	C	T	100.0%	121.9 / NA	38	*295* (TAG→TAA)	xapR	transcriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/24);  total (14/24)

AAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGCCAGTTCAT  >  NC_000913/2521458‑2521718                                                                                                                                        |                                                                                                                              aagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTAc                                                                                                                            <  2:435475/139‑1 (MQ=255)                 aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGTGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCa                                                                                                            <  1:775835/139‑1 (MQ=255)                               tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATc                                                                                              >  2:612945/1‑139 (MQ=255)                                       cgcCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGt                                                                                      <  2:780204/139‑1 (MQ=255)                                       cgcCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGt                                                                                      >  1:521189/1‑139 (MQ=255)                                       cgcCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGt                                                                                      >  1:54658/1‑139 (MQ=255)                                             aaaGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACctg                                                                                <  2:54658/139‑1 (MQ=255)                                                    cTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGtt                                                                         >  2:173543/1‑139 (MQ=255)                                                        tttGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCAt                                                                     <  2:965559/139‑1 (MQ=255)                                                         ttGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATa                                                                    <  1:1023078/139‑1 (MQ=255)                                                         ttGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATa                                                                    <  1:28471/139‑1 (MQ=255)                                                           gCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                  <  2:521189/139‑1 (MQ=255)                                                             tATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACa                                                                <  2:341658/139‑1 (MQ=255)                                                               tCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCtgt                                                                            >  2:902035/1‑125 (MQ=255)                                                               tCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCtgt                                                                            <  1:902035/125‑1 (MQ=255)                                                                  gCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAatata                                                           <  2:46087/139‑1 (MQ=255)                                                                     aTTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATaa                                                        <  1:846237/139‑1 (MQ=255)                                                                                    tGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATgcg                                         >  1:276245/1‑139 (MQ=255)                                                                                            ttCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                  <  1:722507/106‑1 (MQ=255)                                                                                            ttCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                  >  2:722507/1‑106 (MQ=255)                                                                                                cTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAggg                                              <  1:940201/122‑1 (MQ=255)                                                                                                cTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAggg                                              >  2:940201/1‑122 (MQ=255)                                                                                                     cAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAatatat                                                          >  1:870236/1‑105 (MQ=255)                                                                                                     cAGACTACTATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAatatat                                                          <  2:870236/105‑1 (MQ=255)                                                                                                        actactATTTCCGGCAATTCCTGTCTCCTCACTTa                                                                                                                             <  2:993032/35‑1 (MQ=255)                                                                                                        actactATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                  >  1:742950/1‑94 (MQ=255)                                                                                                        actactATTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                  <  2:742950/94‑1 (MQ=255)                                                                                                           actaTTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAAt                                                      >  2:997048/1‑103 (MQ=255)                                                                                                           actaTTTCCGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAAt                                                      <  1:997048/103‑1 (MQ=255)                                                                                                                   cGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGc          <  2:553000/139‑1 (MQ=255)                                                                                                                   cGGCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGc          <  1:173543/139‑1 (MQ=255)                                                                                                                    ggCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCa                                                                                                            <  2:133981/40‑1 (MQ=255)                                                                                                                    ggCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCa                                                                                                            >  1:133981/1‑40 (MQ=255)                                                                                                                    ggCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCtgt                                                                            <  2:431929/72‑1 (MQ=255)                                                                                                                    ggCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCtgt                                                                            >  1:431929/1‑72 (MQ=255)                                                                                                                    ggCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                  <  2:34671/82‑1 (MQ=255)                                                                                                                    ggCAATTCCTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATaaa                                                                  >  1:34671/1‑82 (MQ=255)                                                                                                                           ccTGTCTCCTCACTTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGCCAGTTCAt  <  1:612945/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                              AAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTCATAAACAATATATAAATCTGCAGGGATGCGCTGTTTGAGCGGACGGAAAATGACACCTGGCCAGTTCAT  >  NC_000913/2521458‑2521718

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