Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	599,761	G→A	A350T (GCG→ACG)	cusA →	copper/silver export system RND permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	599,761	0	G	A	96.9%	103.3 / ‑6.4	32	A350T (GCG→ACG)	cusA	copper/silver export system RND permease
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (16/15);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGAAGAGTGGCAGCACCAGCATCCTGACGCCACGCTGGATA  >  NC_000913/599531‑600007                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGttt                                                                                                                                                                                                                                    <  1:209675/251‑1 (MQ=255)      aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTg                                                                                                                                                                                                                               <  1:122179/251‑1 (MQ=255)                  aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCGCTGGCATGTGCGCTc                                                                                                                                                                                                                   <  1:271248/251‑1 (MQ=255)                    cGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATTCTGCGATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAAAAACTGGAAAAGCTGAAAAGGAGGCTGCCGGGAGGCGTGGAGATAGTTACAAAATACGATCGCCGCCAGCTCCTTTACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGGAGAGGTTTTTGTGGTGGCGGGAGTCTGTGGGCTGTTTCTCTGGCATGTGCGCTCgg                                                                                                                                                                                                                 >  1:74355/1‑251 (MQ=255)                     ggcgaaggcgaaGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGgc                                                                                                                                                                                                                <  2:331284/251‑1 (MQ=255)                                   ggcgggcgGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGCGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGatt                                                                                                                                                                                                   >  2:296215/1‑250 (MQ=255)                                                               ggCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGc                                                                                                                                                                      >  1:142567/1‑251 (MQ=255)                                                                            gagaAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATg                                                                                                                                                         >  1:283988/1‑251 (MQ=255)                                                                            gagaAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTAGTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATg                                                                                                                                                         <  1:332071/251‑1 (MQ=255)                                                                                         cgccgTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGAc                                                                                                                                            <  2:5577/251‑1 (MQ=255)                                                                                                                              cTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGAt                                                                                                       >  2:148071/1‑251 (MQ=255)                                                                                                                                      aGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCg                                                                                               <  1:160415/251‑1 (MQ=255)                                                                                                                                      aGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCg                                                                                               >  2:216291/1‑251 (MQ=255)                                                                                                                                       ggCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCgg                                                                                              >  1:8521/1‑251 (MQ=255)                                                                                                                                          gTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGc                                                                                           >  1:342968/1‑251 (MQ=255)                                                                                                                                                   gTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGACGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCg                                                                                   >  2:177852/1‑250 (MQ=255)                                                                                                                                                    ttACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGAt                                                                                 >  1:69016/1‑251 (MQ=255)                                                                                                                                                           aTACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATgctgct                                                                           >  1:62411/1‑250 (MQ=255)                                                                                                                                                              cGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCg                                                                       <  1:272627/251‑1 (MQ=255)                                                                                                                                                                     agccagcTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGAt                                                                <  2:315166/251‑1 (MQ=255)                                                                                                                                                                      gccagcTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATc                                                               <  1:306249/251‑1 (MQ=255)                                                                                                                                                                          gcTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCTTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCATGTACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGTGCGCTGGTCGATGCTGCTATCGTCATGATCg                                                              >  1:5284/1‑248 (MQ=255)                                                                                                                                                                              agtgacggcgacatcagcacccTCTGCGGCAAGTTGCTGGAAGAGGTTATTGTGGTTACGGTAGTCCGTGCGCTGTTTCTCTGGCAGGTGCGCGCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATgc                                                       <  2:185574/231‑1 (MQ=255)                                                                                                                                                                                  aCCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCAt                                                   <  1:127246/251‑1 (MQ=255)                                                                                                                                                                                    cgcgCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATa                                                  >  2:20820/1‑250 (MQ=255)                                                                                                                                                                                         cATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGCGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGCTTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGc                                            >  1:241352/1‑251 (MQ=255)                                                                                                                                                                                              aCAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGa                                        >  2:152269/1‑250 (MQ=255)                                                                                                                                                                                                                           tttATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGAGCATAAACGGCTGGAAGAGTGGCAGCACCAGCATCCTGACGCCAc          >  2:364828/1‑251 (MQ=255)                                                                                                                                                                                                                            ttATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGCTTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGAAGATTGGCAGCTCCAGCATCCTGACGCCAc          <  1:100142/250‑1 (MQ=255)                                                                                                                                                                                                                            ttATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGAAGAGTGGCAGCACCAGCATCCTGACGCCAc          <  1:162181/250‑1 (MQ=255)                                                                                                                                                                                                                            ttATTGTGGTGACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGAAGAGGGGCAGCACCAGCATCCTGACGCCAc          <  1:81231/250‑1 (MQ=255)                                                                                                                                                                                                                                 gtggtgACGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGAAGAGTGGCAGCACCAGCATCCTGACGCCACGCTgg     >  1:220205/1‑250 (MQ=255)                                                                                                                                                                                                                                   ggtgaCGGTAGTCTGTGCGCGGGTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGAAGAGTGGCAGCACCAGCATCCTGACGCCACGCTGGata  <  2:274796/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCAGCGGCAAGTTGCTGGAAGAGTTTATTGTGGTGGCGGTAGTCTGTGCGCTGTTTCTCTGGCATGTGCGCTCGGCGCTGGTGGCGATTATTTCGTTGCCGCTGGGGTTGTGTATTGCTTTTATTGTCATGCACTTCCAGGGACTGAATGCCAATATTATGTCGTTGGGTGGCATTGCGATTGCCGTCGGGGCGATGGTCGATGCTGCTATCGTCATGATCGAGAATGCGCATAAACGGCTGGAAGAGTGGCAGCACCAGCATCCTGACGCCACGCTGGATA  >  NC_000913/599531‑600007

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