Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,606,456	G→T	S2S (TCG→TCT)	zntA →	zinc, cobalt and lead efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,606,456	0	G	T	100.0%	155.3 / NA	48	S2S (TCG→TCT)	zntA	zinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (27/21);  total (27/21)

CGGCCTGCTACTTTGCCGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCGACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCACGCCAACTCTCTCTGAAAACGTCTCCGGC  >  NC_000913/3606325‑3606591                                                                                                                                    |                                                                                                                                        cGGCCTGCTACTTTGCCGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTg                                                                                                                                  >  2:746712/1‑139 (MQ=255)      tGCTACTTTGCCGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             <  2:995556/139‑1 (MQ=255)      tGCTACTTTGCCGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             >  2:1484114/1‑139 (MQ=255)      tGCTACTTTGCCGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             <  1:1484114/139‑1 (MQ=255)      tGCTACTTTGCCGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             <  1:746712/139‑1 (MQ=255)                ccGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGGCTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             <  1:951362/129‑1 (MQ=255)                ccGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             >  2:951362/1‑129 (MQ=255)                  ggTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGaaa                                                                                                                 >  2:41411/1‑139 (MQ=255)                     cACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcc                                                                                                              <  1:1268451/139‑1 (MQ=255)                                  tCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTtgctg                                                                                                 >  2:1209219/1‑139 (MQ=255)                                  tCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTtgctg                                                                                                 >  1:807398/1‑139 (MQ=255)                                  tCCGCTCACTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcccc                                                                                                            >  1:1073715/1‑128 (MQ=255)                                  tCCGCTCACTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcccc                                                                                                            <  2:1073715/128‑1 (MQ=255)                                     gctcgctGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGt                                                                                              >  1:1387743/1‑139 (MQ=255)                                                  ctGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAAc                                                                                 <  1:1209219/139‑1 (MQ=255)                                                   tGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAAcc                                                                                <  2:807398/139‑1 (MQ=255)                                                     aTAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             >  1:5407/1‑92 (MQ=255)                                                     aTAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             <  2:5407/92‑1 (MQ=255)                                                          aCTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTAc                                                                         >  2:1504494/1‑139 (MQ=255)                                                             acaCTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             >  2:1502562/3‑84 (MQ=255)                                                             acaCTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAAt                                                                                                                             <  1:1502562/82‑1 (MQ=255)                                                                   tGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGaa                                                                     <  1:1439666/134‑1 (MQ=255)                                                                   tGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGaa                                                                     >  2:1439666/1‑134 (MQ=255)                                                                       gTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGa                                                            <  2:1387743/139‑1 (MQ=255)                                                                           aCTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGt                                                        <  1:682834/139‑1 (MQ=255)                                                                                agagTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGa                                                                      <  1:314644/120‑1 (MQ=255)                                                                                agagTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGa                                                                      >  2:314644/1‑120 (MQ=255)                                                                                         ccTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGgcgc                                          >  2:520149/1‑139 (MQ=255)                                                                                            tCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcc                                                                                                              >  1:224542/1‑68 (MQ=255)                                                                                            tCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcc                                                                                                              <  2:224542/68‑1 (MQ=255)                                                                                             cGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGt                                      >  2:1151991/1‑139 (MQ=255)                                                                                              ggTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAg                                                                                                                >  2:456294/1‑64 (MQ=255)                                                                                              ggTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAg                                                                                                                <  1:456294/64‑1 (MQ=255)                                                                                              ggTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGtt                                     >  1:239919/1‑139 (MQ=255)                                                                                              ggTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGtt                                     >  2:620988/1‑139 (MQ=255)                                                                                                ttAATGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTcc                                   >  1:1283034/1‑139 (MQ=255)                                                                                                  aaTGAGAAAAAACTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAg                                 >  1:859873/1‑139 (MQ=255)                                                                                                         aaaaaCTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcccc                                                                                                            >  1:792644/1‑57 (MQ=255)                                                                                                         aaaaaCTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcccc                                                                                                            <  2:792644/57‑1 (MQ=255)                                                                                                             aCTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcc                                                                                                              >  1:1391686/1‑51 (MQ=255)                                                                                                             aCTTAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGcc                                                                                                              <  2:1391686/51‑1 (MQ=255)                                                                                                                tAACCGGAGGATGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCATGCCAActctct                   >  2:463544/1‑139 (MQ=255)                                                                                                                       aggaTGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCACGCCAACTCTCTCTGaaaa            >  2:711905/1‑139 (MQ=255)                                                                                                                         gaTGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCt                                                                                    >  2:1440720/1‑65 (MQ=255)                                                                                                                         gaTGCCATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCt                                                                                    <  1:1440720/65‑1 (MQ=255)                                                                                                                             ccATGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCACGCCAACTCTCTCTGAAAACGtctc      <  1:711905/139‑1 (MQ=255)                                                                                                                               aTGTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCACGCCAACTCTCTCTGAAAACGTCTCCg    <  2:859873/139‑1 (MQ=255)                                                                                                                                 gTCTACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCACGCCAACTCTCTCTGAAAACGTCTCCGGc  >  2:743533/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                        CGGCCTGCTACTTTGCCGGTCACTTCCTGATCGTCCGCTCGCTGTATCTCTGATAAAACTTGACTCTGGAGTCGACTCCAGAGTGTATCCTTCGGTTAATGAGAAAAAACTTAACCGGAGGATGCCATGTCGACTCCTGACAATCACGGCAAGAAAGCCCCTCAATTTGCTGCGTTCAAACCGCTAACCACGGTACAGAACGCCAACGACTGTTGCTGCGACGGCGCATGTTCCAGCACGCCAACTCTCTCTGAAAACGTCTCCGGC  >  NC_000913/3606325‑3606591

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