Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,183,006	G→A	G153S (GGT→AGT)	zupT →	zinc transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,183,006	0	G	A	100.0%	111.9 / NA	35	G153S (GGT→AGT)	zupT	zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/16);  total (19/16)

TTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAGGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGC  >  NC_000913/3182878‑3183131                                                                                                                                 |                                                                                                                              ttCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCtggcagt                                                                                                                     >  1:315676/1‑139 (MQ=255)                gCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTa                                                                                                      >  2:814105/1‑139 (MQ=255)                     aGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCgg                                                                                                 >  2:361303/1‑139 (MQ=255)                         tGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAAc                                                                                             <  1:361303/139‑1 (MQ=255)                                             aTTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCgg                                                                                                 <  1:154073/115‑1 (MQ=255)                                             aTTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCgg                                                                                                 >  2:154073/1‑115 (MQ=255)                                             aTTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTaaa                                                                                   <  1:208786/129‑1 (MQ=255)                                             aTTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTaaa                                                                                   >  2:208786/1‑129 (MQ=255)                                                  cACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCtg                                                                    >  2:809038/1‑139 (MQ=255)                                                          tCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCggg                                                            <  2:74741/139‑1 (MQ=255)                                                                  agcagcAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTg                                                    >  1:340285/1‑139 (MQ=255)                                                                               aGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGaa                                         >  1:518536/1‑137 (MQ=255)                                                                               aGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGaa                                         <  2:518536/137‑1 (MQ=255)                                                                                    ggATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtg                                  >  2:377067/1‑139 (MQ=255)                                                                                    ggATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtg                                  >  2:161638/1‑139 (MQ=255)                                                                                            cATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGtgtg                          <  1:232549/139‑1 (MQ=255)                                                                                                      gCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTaaa                                                                                   >  2:620238/1‑72 (MQ=255)                                                                                                      gCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTaaa                                                                                   <  1:620238/72‑1 (MQ=255)                                                                                                      gCCGTCGCGTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGt                >  1:713025/1‑139 (MQ=255)                                                                                                             cgTTGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCt         <  2:729342/139‑1 (MQ=255)                                                                                                                tGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGt                <  1:440667/129‑1 (MQ=255)                                                                                                                tGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGt                >  2:440667/1‑129 (MQ=255)                                                                                                                tGCACAATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCgg      >  1:541483/1‑139 (MQ=255)                                                                                                                  cacaATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtggtg                               >  1:565767/1‑112 (MQ=255)                                                                                                                  cacaATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTtggtg                               <  2:565767/112‑1 (MQ=255)                                                                                                                    caATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCaa                                                                                              <  1:498614/47‑1 (MQ=255)                                                                                                                    caATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCaa                                                                                              >  2:498614/1‑47 (MQ=255)                                                                                                                    caATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGc  <  1:809038/139‑1 (MQ=255)                                                                                                                    caATATCCCTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGc  <  2:713025/139‑1 (MQ=255)                                                                                                                          cccTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCgg      >  2:790603/1‑129 (MQ=255)                                                                                                                          cccTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCgg      <  1:790603/129‑1 (MQ=255)                                                                                                                            cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAAt            >  1:270308/1‑121 (MQ=255)                                                                                                                            cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAAt            <  2:270308/121‑1 (MQ=255)                                                                                                                            cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAg   >  1:829228/1‑130 (MQ=255)                                                                                                                            cTGAAAGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAg   <  2:829228/130‑1 (MQ=255)                                                                                                                                 |                                                                                                                              TTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCTGGAGCTGGGATTTGGCATCGCACTGGCCGTCGCGTTGCACAATATCCCTGAAGGTCTGGCAGTGGCAGGCCCGGTTTATGCGGCAACGGGGTCTAAACGCACCGCGATTCTGTGGGCGGGGATTTCTGGACTAGCAGAAATTCTTGGTGGTGTGCTGGCGTGGTTAATCCTCGGTAGC  >  NC_000913/3182878‑3183131

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