Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	784,869	A→G	intergenic (‑46/+68)	zitB ← / ← ybgS	Zn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter/PF13985 family protein YbgS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	784,869	0	A	G	100.0%	179.0 / NA	52	intergenic (‑46/+68)	zitB/ybgS	Zn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter/PF13985 family protein YbgS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/25);  total (27/25)

AACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCACCGGTTTGTACCTTTTT  >  NC_000913/784732‑785005                                                                                                                                          |                                                                                                                                         aaCCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCt                                                                                                                                        <  1:158648/140‑1 (MQ=255)       gCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCtttt                                                                                                                                  >  1:414895/1‑140 (MQ=255)       gCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCtttt                                                                                                                                  >  1:439752/1‑140 (MQ=255)               gCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGa                                                                                                                          <  1:359757/140‑1 (MQ=255)                         tACAACAAGCGACGAGCATTATTATCTTCTAGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTaaaa                                                                                                                >  1:787956/1‑140 (MQ=255)                                 gCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGgagaga                                                                                                        >  2:662873/1‑140 (MQ=255)                                      gaGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCg                                                                                                   <  2:248384/140‑1 (MQ=255)                                       aGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGc                                                                                                  <  1:371027/140‑1 (MQ=255)                                       aGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGc                                                                                                  <  1:76781/140‑1 (MQ=255)                                        gCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTCTTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGct                                                                                                 <  2:787956/140‑1 (MQ=255)                                        gCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGct                                                                                                 <  2:414895/140‑1 (MQ=255)                                                   ttctGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCattatt                                                                                      >  1:12317/1‑140 (MQ=255)                                                   ttctGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCattatt                                                                                      >  2:64880/1‑140 (MQ=255)                                                      tGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCa                                                                                           >  2:154525/1‑132 (MQ=255)                                                      tGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCa                                                                                           <  1:154525/132‑1 (MQ=255)                                                      tGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGc                                                                                   >  1:499827/1‑140 (MQ=255)                                                      tGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGc                                                                                   >  1:669860/1‑140 (MQ=255)                                                      tGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGc                                                                                   >  2:399110/1‑140 (MQ=255)                                                        gCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAATTTAAAAGGGAGAGCCTCCGCTCTCCCATTATTGGCtt                                                                                 >  1:59533/1‑139 (MQ=255)                                                             tgtgAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAGAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTg                                                                            <  2:762004/140‑1 (MQ=255)                                                                   gACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGt                                                                      <  1:399110/140‑1 (MQ=255)                                                                   gACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGt                                                                      <  1:238000/140‑1 (MQ=255)                                                                      gtgtgtGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAggg                                                                       <  1:388528/136‑1 (MQ=255)                                                                      gtgtgtGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAggg                                                                       >  2:388528/1‑136 (MQ=255)                                                                       tgtgtgAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTa                                                                    >  2:696976/1‑138 (MQ=255)                                                                       tgtgtgAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTa                                                                    <  1:696976/138‑1 (MQ=255)                                                                       tgtgtgAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACt                                                                  <  2:190969/140‑1 (MQ=255)                                                                       tgtgtgAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACt                                                                  >  2:605130/1‑140 (MQ=255)                                                                          gtgagtgtgagtGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCg                                                               <  2:542586/140‑1 (MQ=255)                                                                            gagtgtgagtgCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTg                                                             <  2:669860/140‑1 (MQ=255)                                                                                 tgagtgCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTAcc                                                        >  2:4616/1‑140 (MQ=255)                                                                                       cgcCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGt                                                  >  2:717769/1‑140 (MQ=255)                                                                                        gcCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGtt                                                 >  2:770213/1‑140 (MQ=255)                                                                                         ccATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGttt                                                >  2:570158/1‑140 (MQ=255)                                                                                            tGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTgg                                             >  1:105963/1‑140 (MQ=255)                                                                                               cGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTAt                                          >  2:532975/1‑140 (MQ=255)                                                                                                  gTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGa                                       <  1:570158/140‑1 (MQ=255)                                                                                                       cccTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGAcatca                                  <  2:53692/140‑1 (MQ=255)                                                                                                       cccTTAATTATTCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGAcatca                                  <  1:717769/140‑1 (MQ=255)                                                                                                               tattCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTAt                                          <  1:512304/124‑1 (MQ=255)                                                                                                               tattCTCACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGGAt                                          >  2:512304/1‑124 (MQ=255)                                                                                                                   ctcACTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCa                      >  2:674789/1‑140 (MQ=255)                                                                                                                      aCTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCa                   >  1:354909/1‑140 (MQ=255)                                                                                                                       cTTAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCAc                  <  1:605130/140‑1 (MQ=255)                                                                                                                        ttAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCa                   >  2:857847/1‑138 (MQ=255)                                                                                                                        ttAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCa                   <  1:857847/138‑1 (MQ=255)                                                                                                                         tAACTTTATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCACCg                <  2:12317/140‑1 (MQ=255)                                                                                                                              ttATGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCACCGGTTTg           <  1:487557/140‑1 (MQ=255)                                                                                                                                aTGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGTTTTTGGTATCGACATCATTATTGATGCCATCACCGGTTTGTa         >  2:800735/1‑140 (MQ=255)                                                                                                                                aTGACATCATGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCACCGGTTTGTa         <  1:662873/140‑1 (MQ=255)                                                                                                                                    catcatGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCACCGGTTTGTACCtt     >  1:792678/1‑140 (MQ=255)                                                                                                                                       catGCTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCACCACCGGTTTGTACCttttt  >  2:632128/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                         AACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGGCTATTTTGCAGGGTTACTGCGTGGTACCGTCGGTTTTGGTATCGACATCATTATTGATGCCATCACCGGTTTGTACCTTTTT  >  NC_000913/784732‑785005

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