Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,608,682	T→C	intergenic (+33/+69)	zntA → / ← tusA	Zn(2(+))/Cd(2(+))/Pb(2(+)) exporting P‑type ATPase/sulfur transfer protein TusA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,608,682	0	T	C	100.0%	111.0 / NA	32	intergenic (+33/+69)	zntA/tusA	Zn(2(+))/Cd(2(+))/Pb(2(+)) exporting P‑type ATPase/sulfur transfer protein TusA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/18);  total (14/18)

TTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCA  >  NC_000913/3608554‑3608809                                                                                                                                 |                                                                                                                                ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGG‑GATAAGGCAAACAGATCGCAGCATAGAACGCGATCGGCCCCCt                                                                                                                            >  1:228751/1‑133 (MQ=255)    gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGTTCGCAACCTTGAGCCCGATCGGCCCCCTCg                                                                                                                          >  1:254114/1‑133 (MQ=255)    gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCg                                                                                                                          <  2:555212/133‑1 (MQ=255)      gATGACCGGGTTGTGGCTGGCACTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGcc                                                                                                                        <  1:269974/133‑1 (MQ=255)      gATGACCGGGTTGTGCCTTGCAATGCTGGCGGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGcc                                                                                                                        <  1:202861/133‑1 (MQ=255)          aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCctc                                                                                                                    >  2:65480/1‑133 (MQ=255)           ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCctct                                                                                                                   >  1:420450/1‑133 (MQ=255)                ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGga                                                                                                              <  2:601854/133‑1 (MQ=255)                ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGga                                                                                                              <  1:448513/133‑1 (MQ=255)                    ggctggcAGTGCTGGCAGCTACCGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAgg                                                                                                          <  1:207544/133‑1 (MQ=255)                                       tACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                       >  1:482215/1‑133 (MQ=255)                                       tACGGGGGCGACGGTCCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCCGGGGAGAGGGTTCGGGTGAGGGGAAAAg                                                                                       >  2:249373/1‑133 (MQ=255)                                            gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                  >  2:333404/1‑133 (MQ=255)                                            gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                  >  1:473825/1‑133 (MQ=255)                                              gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                <  2:67716/133‑1 (MQ=255)                                                        tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAAt                                                                      <  1:428499/133‑1 (MQ=255)                                                               aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                               >  2:7265/1‑133 (MQ=255)                                                               aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                               <  1:507669/133‑1 (MQ=255)                                                               aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                               >  1:326275/1‑133 (MQ=255)                                                                    aTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                          >  1:497840/1‑133 (MQ=255)                                                                                  ttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                            >  1:445207/1‑133 (MQ=255)                                                                                  ttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                            >  1:374335/1‑133 (MQ=255)                                                                                                  gCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                            <  1:65480/133‑1 (MQ=255)                                                                                                  gCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                            <  2:473825/133‑1 (MQ=255)                                                                                                  gCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                            <  1:333404/133‑1 (MQ=255)                                                                                                    aaaCCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                          <  1:119413/133‑1 (MQ=255)                                                                                                             gCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                 <  2:482215/133‑1 (MQ=255)                                                                                                              cAACATTGGGCGCGAACGTCCCGCTCGCCCTTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                <  1:188141/133‑1 (MQ=38)                                                                                                              cAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                <  2:48901/133‑1 (MQ=255)                                                                                                                        cgcgATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa      <  2:326275/133‑1 (MQ=255)                                                                                                                         gcgATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa     <  1:7265/133‑1 (MQ=255)                                                                                                                            aTCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa  >  1:379321/1‑133 (MQ=255)                                                                                                                                 |                                                                                                                                TTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCA  >  NC_000913/3608554‑3608809

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