Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	509,517	C→T	E621E (GAG→GAA)	copA ←	Cu(+) exporting P‑type ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	509,517	0	C	T	100.0%	52.6 / NA	19	E621E (GAG→GAA)	copA	Cu(+) exporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/5);  total (14/5)

CTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTAC‑GCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAAT  >  NC_000913/509276‑509733                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                         cTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTCGCCGATCGCATTGGCCGTGGCTGGGTTATCCCCGGTAAACATCAACAGACGATTTTCCGATTTCTTCCCGCGTTGCCGCGCCGCCCCGCTACCCCTCC‑GCAACGGATCGCGTACTGCCAGCCGGGCTACCGCTTTCCCGCCAACCGCCAGCACCACAGGCGTTGCCCCCTGCGATGCCTGAGCAGTAATTTCCGCTTct                                                                                                                                                                                                                 >  1:290445/1‑250 (MQ=255)                           accACCTCATCAACCCCTGCTTCTTTGACGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTAC‑GCAACGGATTGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGc                                                                                                                                                                                                                     <  1:116832/219‑1 (MQ=255)                             naCCTCATCAATCCCTGCTTCTTTGACGATCGCATTGGCGGTGGGTGGGGTATCCCCGGGCAACATCACCAGGCGATATCCCGCTTTATGCCGGCGTTGCAGCGCCGCCACGCTATCACTACCGCA‑CGGATCGCGTACTGCCAGCAGGGATACCGCTTTCCCGTCAACCGCCCGCAGAACCGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTTCGATTCGGTAGCTTTGGTAACCACCTGTTGCCCAt                                                                                                                                                                                     >  1:43668/2‑251 (MQ=255)                                                                                                                                    gccgccACGCTATCACTAC‑GCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTACCCGTCAAACCCCAGCAGCACAGGCGTTGCCACTTGCGATGCCTGCACAGTAATTTCCGCTTCGATAGCTTTGGTACCAAACTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCAc                                                                                                                                     >  1:88265/1‑196 (MQ=255)                                                                                                                                             cTATCACTAC‑GCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGc                                                                     >  1:219483/1‑251 (MQ=255)                                                                                                                                                      c‑GCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTg                                                                                                                                                                                          <  1:142663/125‑1 (MQ=255)                                                                                                                                                         cAACGGATCGCATACTGCCAGCCGGGCTACCGCTTTCCCGTCAACCGCAAGCAGCACAGGCTTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGCCCTTCAGCTTCCCCGCTCACGCCCAGCCCGCGCACTGTGCGGAAAcc                                                                                                   >  1:184574/1‑210 (MQ=255)                                                                                                                                                             ggATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGc                                                                                                                                                                                         <  1:439373/120‑1 (MQ=255)                                                                                                                                                                 cncgTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCa                                                                                                                                                                                      >  1:396333/1‑119 (MQ=255)                                                                                                                                                                       cTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCCCAGGCGTTGCACCTTGCGATGACTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCa                                                                               >  1:267795/1‑216 (MQ=255)                                                                                                                                                                         gCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTan                                                                                                                                                                                                    <  1:313657/97‑2 (MQ=255)                                                                                                                                                                         gCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTAc                                                                                                                                                                                                    >  2:313657/1‑97 (MQ=255)                                                                                                                                                                         gCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAAATTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTAATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGa                                                                                                                                                  >  1:328340/1‑147 (MQ=255)                                                                                                                                                                                ggtcTACCGCTTTCCCGTCAACCGCCAGCCGCACAGGCCTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTa                                                                                                                                                                                                     >  1:455786/4‑89 (MQ=255)                                                                                                                                                                                        gCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTCGCGATGCCTGAGCAGTCATTTCCGCTTCGATACCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATCCCGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCAGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCCGCGGATGGCTGGAAcc                              >  1:288603/1‑248 (MQ=255)                                                                                                                                                                                                       nccAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCCCCGCTCACGCCCAGCCcgcg                                                                                                                  >  1:123346/2‑149 (MQ=255)                                                                                                                                                                                                                           cccTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGTCCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAAt  >  1:496627/1‑241 (MQ=255)                                                                                                                                                                                                                            ccTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCa                                                                                                                                                                >  2:34939/1‑82 (MQ=255)                                                                                                                                                                                                                            ccTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGACTGATTGCCCa                                                                                                                                                                <  1:34939/82‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                         CTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTAC‑GCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAAT  >  NC_000913/509276‑509733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑