Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	1,959,198	G→A	D23D (GAC→GAT)	cutC ←	copper homeostasis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	1,959,198	0	G	A	100.0%	93.9 / NA	30	D23D (GAC→GAT)	cutC	copper homeostasis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/19);  total (11/19)

TCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGCCGACTGCTCGCTGGCAACGATCTCTTCAATGGACCACGGATGAAACTTAATAGTTGT  >  NC_000913_3_vch_tpiA/1959055‑1959342                                                                                                                                                |                                                                                                                                                 tCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGGTGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCACCCCTCTTTTGGGGCTGCGCATAATTGATGtctatctga                                                                                                                                             >  1:384049/1‑148 (MQ=255)   aCCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTgcgc                                                                                                                                           >  1:488641/1‑149 (MQ=255)           tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTtgc                                                                                                                                   >  1:9299/1‑149 (MQ=255)                cAAAAATCACCACCGCGTCGGCGAATTATCGGACGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTgcgc                                                                                                                              <  1:57764/149‑1 (MQ=255)                       ccccaccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGt                                                                                                                       <  1:303538/147‑1 (MQ=255)                       caccacCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAAGTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGt                                                                                                                       <  2:235599/149‑1 (MQ=255)                        accaccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTg                                                                                                                      <  2:274552/149‑1 (MQ=255)                          caccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGca                                                                                                                    >  2:457425/1‑149 (MQ=255)                                   ggCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGAGTTAAGCCCCCCTCTTTCGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATg                                                                                                           >  1:444107/1‑149 (MQ=255)                                     cGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCt                                                                                                         <  1:362871/149‑1 (MQ=255)                                       aaTTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGt                                                                                                       <  1:457425/149‑1 (MQ=255)                                         ttATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTaa                                                                                                     <  2:488641/149‑1 (MQ=255)                                          tATCGGATGCACAGGGATCGTCACCCGCTGGCGGACTGATTTCAGTACACCCAACGACGGCGTTAAG‑CCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACCTTCCATGCTGTAACa                                                                                                   >  1:249598/1‑149 (MQ=255)                                          tATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGCTTCCTGTACGCTCTACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAAc                                                                                                    <  1:39138/149‑1 (MQ=255)                                           aTCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACa                                                                                                   <  2:190308/149‑1 (MQ=255)                                                    acaGGGATCGTCACCCCCTGGCGCACGGATTTCAGTACACCCAACGACGGCCTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCCACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTc                                                                                          <  2:263973/149‑1 (MQ=255)                                                       gggATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGt                                                                                                                           <  1:433682/113‑1 (MQ=255)                                                       gggATAGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGt                                                                                                                           >  2:433682/1‑113 (MQ=255)                                                                 cccGCTGGCGCACGGATTTCAGGACACCCAACAACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCAAAATTCAACTATATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATGCCATGCTGTGACAGCCAATTTCCAGTCATGCCAtc                                                                             >  2:9042/1‑148 (MQ=255)                                                                        gcgcACGGATTTCAGGACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATACAAttcgtgcgc                                                                      <  2:258475/149‑8 (MQ=255)                                                                         cgcACGGATTTCAGTACACCCACCGACGTCGTTAAG‑CCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCt                                                                    <  2:249598/149‑1 (MQ=255)                                                                               gATTTCAGTACACCCAACGACGGCGTTAAGCCGCCTTCTTTTGGGGCTGCGAATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAAtt                                                               <  2:159373/149‑1 (MQ=255)                                                                                           cccAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAGCTCTATCTGCGCCGATTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGCCGACTGCt                                                   >  1:439265/1‑149 (MQ=255)                                                                                                        ttAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAGCTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGCCGACTGCTCGCTGGCAACGAt                                      <  1:288165/149‑1 (MQ=255)                                                                                                            gCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGc                                                                                                 >  2:169489/1‑86 (MQ=255)                                                                                                            gCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGc                                                                                                 <  1:169489/86‑1 (MQ=255)                                                                                                             ccccccccTTTTGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGCCGACTGCTCGCTGGCAACGATCTCtt                                 <  1:234340/149‑1 (MQ=255)                                                                                                                     ttttGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGcc                                                          >  2:43753/1‑116 (MQ=255)                                                                                                                     ttttGGGGCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGcc                                                          <  1:43753/116‑1 (MQ=255)                                                                                                                         ggggCTGCGCATAATTCAACTCTATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGCCGACTGCTCGCTGGCAACGATCTCTTCAATGGACCACg                     <  2:9299/149‑1 (MQ=255)                                                                                                                                            ctctATCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGCCGACTGCTCGCTGCCAACGATCTCTTCATTGGACCACGGATGAAACTTAATAgttgt  >  1:355960/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                 TCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCATAATTCAACTCTGTCTGCGCCGTTTTGCTGCGCCGTTAGTGCACATTCCATGCTGTAACAGCAAATTTCCAGTAATGCCATTTTTACTCCTTAATTACGCCGACTGCTCGCTGGCAACGATCTCTTCAATGGACCACGGATGAAACTTAATAGTTGT  >  NC_000913_3_vch_tpiA/1959055‑1959342

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