Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	137,955	A→G	K291K (AAA→AAG)	cueO →	cuprous oxidase/multicopper oxidase with role in copper homeostasis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	137,955	0	A	G	100.0%	93.3 / NA	28	K291K (AAA→AAG)	cueO	cuprous oxidase/multicopper oxidase with role in copper homeostasis
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/13);  total (15/13)

CAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCTTTGCCAGACACATTAAGTA  >  NC_000913/137832‑138088                                                                                                                            |                                                                                                                                      cAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGa                                                                                                                       <  1:613845/140‑1 (MQ=255)        ccGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGcc                                                                                                                <  2:413903/140‑1 (MQ=255)           cTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGt                                                                                                             >  2:156180/1‑140 (MQ=255)                        ccAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGa                                                                                                     >  2:661992/1‑135 (MQ=255)                        ccAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGa                                                                                                     <  1:661992/135‑1 (MQ=255)                        ccAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATgggg                                                                                                >  2:847617/1‑140 (MQ=255)                        ccAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATgggg                                                                                                >  2:152576/1‑140 (MQ=255)                                            ccTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGTCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGa                                                                            >  1:571623/1‑140 (MQ=255)                                              tGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATa                                                                          >  2:140310/1‑140 (MQ=255)                                                 aCCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGc                                                                       <  1:797362/140‑1 (MQ=255)                                                    aGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTc                                                                    >  2:901807/1‑140 (MQ=255)                                                    aGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTc                                                                    >  2:208053/1‑140 (MQ=255)                                                               gCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGATGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCg                                                               >  2:44054/1‑134 (MQ=255)                                                               gCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGATGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCg                                                               <  1:44054/134‑1 (MQ=255)                                                               gCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATg                                                         >  2:371853/1‑140 (MQ=255)                                                                cGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGc                                                        >  2:65511/1‑140 (MQ=255)                                                                  aaCTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCgg                                                      <  1:901807/140‑1 (MQ=255)                                                                       ccGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCa                                                 >  2:955773/1‑140 (MQ=255)                                                                             cTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGAt                                           <  1:175309/140‑1 (MQ=255)                                                                                      gAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAg                                  <  1:152576/140‑1 (MQ=255)                                                                                       aGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGt                                 <  1:65511/140‑1 (MQ=255)                                                                                                    gCTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCtt                    >  2:667981/1‑140 (MQ=255)                                                                                                     cTGGTGGAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCttt                   >  2:263485/1‑140 (MQ=255)                                                                                                          ggAGGTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCTTTGCCAg              <  1:140310/140‑1 (MQ=255)                                                                                                              gTTAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCTTTGCCAGacac          <  1:667981/140‑1 (MQ=255)                                                                                                               ttAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCTTTGc                 <  1:783597/132‑1 (MQ=255)                                                                                                               ttAACGATAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCTTTGc                 >  2:783597/1‑132 (MQ=255)                                                                                                                      tAACAAGCCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCTTTGCCAGACACATTAAGTa  <  1:263485/140‑1 (MQ=255)                                                                                                                            |                                                                                                                                      CAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGATTGCTATTAGTGCCTCCGGTGCTTTGCCAGACACATTAAGTA  >  NC_000913/137832‑138088

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