Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	592,602	C→T	Q6* (CAA→TAA)	cusF →	Cu(+)/Ag(+) efflux RND transporter periplasmic metallochaperone CusF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	592,602	0	C	T	100.0%	123.0 / NA	32	Q6* (CAA→TAA)	cusF	Cu(+)/Ag(+) efflux RND transporter periplasmic metallochaperone CusF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

TCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTA  >  NZ_CP009273/592470‑592718                                                                                                                                     |                                                                                                                     tctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCg                                                                                                                >  1:267811/1‑139 (MQ=255)    gcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAg                                                                                                                   <  2:665597/133‑1 (MQ=255)    gcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAg                                                                                                                   >  1:665597/1‑133 (MQ=255)     cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGc                                                                                                               <  2:247297/136‑1 (MQ=255)     cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGc                                                                                                               <  1:715186/136‑1 (MQ=255)     cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAAt                                                                                                            >  2:97433/1‑139 (MQ=255)     cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAAt                                                                                                            <  1:97674/139‑1 (MQ=255)              aCCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGt                                                                                                     <  2:310004/137‑1 (MQ=255)                   aGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTa                                                                                              <  2:631996/139‑1 (MQ=255)                    gTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTAc                                                                                             <  2:660423/139‑1 (MQ=255)                       acttacttATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCg                                                                                           <  2:26646/138‑1 (MQ=255)                                  tAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCttt                                                                                >  2:49159/1‑138 (MQ=255)                                     gCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATg                                                                            >  1:842387/1‑139 (MQ=255)                                              gATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTa                                                                   <  2:129571/139‑1 (MQ=255)                                                   aCTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGaa                                                              >  2:407391/1‑139 (MQ=255)                                                       gCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAAcat                                                           <  1:276608/138‑1 (MQ=255)                                                             cGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAAcatcatcat                                                     >  2:822368/1‑138 (MQ=255)                                                              gCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATg                                                    <  1:863753/138‑1 (MQ=255)                                                                              aaTTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGca                                   >  2:758315/1‑139 (MQ=255)                                                                                               cTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAg                  <  1:717874/139‑1 (MQ=255)                                                                                                  ttAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAgcgc               >  2:169418/1‑139 (MQ=255)                                                                                                  ttAACTCCAGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAgcgc               >  1:217746/1‑139 (MQ=255)                                                                                                          aGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAAcat                                                           <  1:730516/87‑1 (MQ=255)                                                                                                          aGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAAcat                                                           >  2:730516/1‑87 (MQ=255)                                                                                                          aGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATg                                                    <  2:626290/94‑1 (MQ=255)                                                                                                          aGGAGAGAATAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATg                                                    >  1:626290/1‑94 (MQ=255)                                                                                                                 aaTAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAg                                                                                                      >  2:632008/1‑37 (MQ=255)                                                                                                                 aaTAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAg                                                                                                      <  1:632008/37‑1 (MQ=255)                                                                                                                 aaTAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAgcgc               >  2:344929/1‑124 (MQ=255)                                                                                                                 aaTAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAgcgc               <  1:344929/124‑1 (MQ=255)                                                                                                                 aaTAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTa  >  2:707469/1‑137 (MQ=255)                                                                                                                 aaTAAATGAAAAAAGCACTGTAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTa  >  1:817466/1‑137 (MQ=255)                                                                                                                                     |                                                                                                                     TCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCTGCAACTCGAATTATTTAGAGTATGACTTTTAACTCCAGGAGAGAATAAATGAAAAAAGCACTGCAAGTCGCAATGTTCAGTCTGTTTACCGTTATTGGCTTTAATGCCCAGGCTAACGAACATCATCATGAAACCATGAGCGAAGCACAACCACAGGTTATTAGCGCCACTGGCGTGGTA  >  NZ_CP009273/592470‑592718

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