Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	135,028	A→C	N487H (AAT→CAT)	cueO →	multicopper oxidase CueO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	135,028	0	A	C	100.0%	87.0 / NA	28	N487H (AAT→CAT)	cueO	multicopper oxidase CueO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/19);  total (9/19)

AACCGCCAGCGGCTCATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACG  >  NZ_CP009273/134948‑135114                                                                                 |                                                                                       aaCCGCCAGCGGCTCATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATg                                                                               <  1:114927/90‑1 (MQ=255)        aGCGGCTCATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGaa                                                                        <  1:18383/90‑1 (MQ=255)             cTCATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAAc                                                                   <  1:186286/90‑1 (MQ=255)             cTCATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAAc                                                                   <  1:264568/90‑1 (MQ=255)             cTCATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAAc                                                                   <  1:455500/90‑1 (MQ=255)               cATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACAt                                                                 <  1:119845/90‑1 (MQ=255)               cATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACAt                                                                 <  1:346271/90‑1 (MQ=255)                            ggAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGgcgc                                                    >  2:175438/1‑90 (MQ=255)                              aaaGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCAc                                                  <  1:175438/90‑1 (MQ=255)                                  aTACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGcc                                              >  1:213229/1‑90 (MQ=255)                                      cGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATct                                          >  1:292448/1‑90 (MQ=255)                                      cGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATct                                          >  1:36706/1‑90 (MQ=255)                                               aGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCa                                 >  2:379608/1‑90 (MQ=255)                                               aGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCa                                 >  1:78440/1‑90 (MQ=255)                                                gAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCAt                                <  1:36695/90‑1 (MQ=255)                                                  aGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGa                              <  1:297842/90‑1 (MQ=255)                                                  aGGTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGa                              <  1:182810/90‑1 (MQ=255)                                                   ggTAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGaa                             <  2:326352/90‑1 (MQ=255)                                                     tAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCAt                                >  1:373365/1‑85 (MQ=255)                                                     tAATGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCAt                                <  2:373365/85‑1 (MQ=255)                                                       aTGTCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATa                         >  1:215736/1‑90 (MQ=255)                                                          tCAGCGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACgg                      <  2:100277/90‑1 (MQ=255)                                                              cGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGgat                  <  1:216296/90‑1 (MQ=255)                                                              cGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGgat                  <  2:354220/90‑1 (MQ=255)                                                              cGAAGTGCTGGTGAAGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGgat                  <  2:95968/90‑1 (MQ=255)                                                                            aGTTTCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTa    >  2:35936/1‑90 (MQ=255)                                                                              tttCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACg  <  1:112061/90‑1 (MQ=255)                                                                              tttCATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACg  <  2:215736/90‑1 (MQ=255)                                                                                 |                                                                                       AACCGCCAGCGGCTCATCGCGCGGGCTGGAAAGATACCGTTAAGGTAGAAGGTAATGTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACG  >  NZ_CP009273/134948‑135114

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