Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	597,758	(G)5→6	coding (280/1224 nt)	cusB →	copper/silver efflux system, membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	597,753	1	.	G	97.6%	151.9 / ‑2.9	42	coding (275/1224 nt)	cusB	copper/silver efflux system, membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (31/10);  total (31/11)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCT‑GGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTACCGT  >  NC_000913/597617‑597897                                                                                                                                          |                                                                                                                                                 gACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACg                                                                                                                                     <  1:4665/151‑1 (MQ=255)  aCCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACgg                                                                                                                                    <  1:309205/151‑1 (MQ=255)  aCCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACgg                                                                                                                                    >  2:547925/1‑151 (MQ=255)   cccAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGc                                                                                                                                   >  1:130333/1‑151 (MQ=255)      aaTGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTAc                                                                                                                                >  1:533162/1‑151 (MQ=255)       aTGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTAcc                                                                                                                               >  2:514829/1‑151 (MQ=255)       aTGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTAcc                                                                                                                               >  2:674178/1‑151 (MQ=255)                aaTACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCAcgcgc                                                                                                                      <  2:765049/151‑1 (MQ=255)                    cgcgGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACAGCTACCGTCACGCGCGGAc                                                                                                                  <  1:498988/151‑1 (MQ=255)                        gtgcGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCt                                                                                                              <  2:131148/148‑1 (MQ=255)                         ttCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTg                                                                                                             >  1:4500/1‑151 (MQ=255)                         ttCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTg                                                                                                             >  1:800864/1‑151 (MQ=255)                              tAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACttt                                                                                                        >  2:382441/1‑151 (MQ=255)                              tAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACttt                                                                                                        >  1:17665/1‑151 (MQ=255)                                  ccAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGcc                                                                                                    >  2:178312/1‑151 (MQ=255)                                             gCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTccc                                                                                         >  1:806102/1‑151 (MQ=255)                                              ccGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCg                                                                                        >  1:317753/1‑151 (MQ=255)                                               cGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCgg                                                                                       <  1:849699/151‑1 (MQ=255)                                                    atggatatggatCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAAt                                                                                  >  2:656815/1‑151 (MQ=255)                                                          atggatCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGt                                                                            >  1:334215/1‑151 (MQ=255)                                                             gatCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTAc                                                                         >  2:294712/1‑151 (MQ=255)                                                                      ccGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTAt                                                                >  2:271144/1‑151 (MQ=255)                                                                              tGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGc                                                        >  2:248404/1‑151 (MQ=255)                                                                               gCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGcc                                                       >  2:764374/1‑151 (MQ=255)                                                                                    tGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATtg                                                   >  2:677418/1‑150 (MQ=255)                                                                                        gagagTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAg                                              <  2:806158/151‑1 (MQ=255)                                                                                           agTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGc                                            >  1:654811/1‑150 (MQ=255)                                                                                              tCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTcgc                                        >  2:325089/1‑151 (MQ=255)                                                                                                tGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCgcgc                                      >  2:791473/1‑151 (MQ=255)                                                                                                  cGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTg                                    <  2:418085/151‑1 (MQ=255)                                                                                                   gTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGc                                   >  1:444200/1‑151 (MQ=255)                                                                                                      tGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCgg                                >  2:424436/1‑151 (MQ=255)                                                                                                      tGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCgg                                >  2:620412/1‑151 (MQ=255)                                                                                                       ggTGCGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCggg                               >  2:212302/1‑151 (MQ=255)                                                                                                                         aCTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTAt             <  1:670072/151‑1 (MQ=255)                                                                                                                                cGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTa      <  1:303944/151‑1 (MQ=255)                                                                                                                                 gCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTa      >  1:460952/1‑150 (MQ=255)                                                                                                                                 gCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTAc     >  2:427829/1‑151 (MQ=255)                                                                                                                                  aaGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTAcc    >  2:538879/2‑151 (MQ=255)                                                                                                                                   aGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTACCg   >  2:255379/1‑151 (MQ=255)                                                                                                                                   aGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTACCg   >  2:587874/1‑151 (MQ=255)                                                                                                                                   aGAATCT‑GGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTACCGt  <  2:665602/151‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                 GACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCT‑GGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTACCGT  >  NC_000913/597617‑597897

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