Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	510,057	C→T	7.4%	V441V (GTG→GTA)	copA ←	copper transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	510,057	0	C	T	7.4%	163.0 / 3.2	54	V441V (GTG→GTA)	copA	copper transporter
Reads supporting (aligned to +/- strand):  ref base C (23/27);  new base T (2/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TCGACATCGGCGTCGCCAGCCCCAGCGCACACGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACT  >  NC_000913/509909‑510204                                                                                                                                                     |                                                                                                                                                    tCGACATCGGCGTCGCCAGCCCCAGCGCACCCGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGCAATACCAGATTGCCGCCCTGACAAGCGCAATAACCACCCCTACCGGc                                                                                                                                                     >  2:177063/1‑149 (MQ=255)      aTCGGCGTCGCCAGCCCCAGCGCACACGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACaaa                                                                                                                                                <  6:166493/149‑1 (MQ=255)                  aGCCCCAGCGCACACGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGGAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGAt                                                                                                                                    <  8:270302/149‑1 (MQ=255)                    ccccAGCGCACACGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGTACAAATACGGCTGAGAttt                                                                                                                                  >  1:236084/1‑149 (MQ=255)                    ccccAGCGCACACGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGTACAAATACGGCTGAGAttt                                                                                                                                  >  7:98089/1‑149 (MQ=255)                        aGCGCACACGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATc                                                                                                                              <  7:3213/149‑1 (MQ=255)                             acaaGGACAGGAATTATTCAGTATCGTGGTGGAAATCACCAGTGTATAGACAATCTGGTGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGAGAATCGCAATAACCACCACTACCGGCACAAATTCGGCTGAGATTTTATCCGCCa                                                                                                                         <  4:5870/149‑1 (MQ=255)                                cGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCt                                                                                                                      <  4:85456/149‑1 (MQ=255)                                       gCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCTCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGAt                                                                                                               >  6:301348/1‑149 (MQ=255)                                         aataatCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCCCTGACAAGCGCAATACCCCCCACTACCGGCACAAATACGGCTGAGATTTTATCAGCCAGCTGACCGAttt                                                                                                             >  1:384210/1‑149 (MQ=255)                                             atCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAACTACCAGATTGCCGCACTGACAAGGGCAATAACCACCACTACCGGCAGAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTg                                                                                                          <  4:117854/148‑1 (MQ=255)                                             atCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTgg                                                                                                         <  8:181193/149‑1 (MQ=255)                                              tCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGc                                                                                                        >  7:174846/1‑149 (MQ=255)                                                 gTACCTTGTTGGCAATCACCATGGTATAGACAATCTGCGGTTCCGGACCACAGAAATACCAGATTGCCGCCCTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTtg                                                                                                     <  4:226290/149‑1 (MQ=255)                                                        ggtggCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGTACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGctct                                                                                              <  8:98089/149‑1 (MQ=255)                                                         gtggCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGTACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGctct                                                                                              <  2:236084/148‑1 (MQ=255)                                                          gggCAATCACCAGGGTATAGACAATCTGCCGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTgg                                                                                            <  6:206305/148‑1 (MQ=255)                                                            gCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGc                                                                                          >  2:38007/1‑149 (MQ=255)                                                                    cAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGgcg                                                                                    >  1:189595/1‑147 (MQ=255)                                                                    cAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGgcg                                                                                    <  2:189595/147‑1 (MQ=255)                                                                           taGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAACCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGc                                                                            >  7:231628/1‑148 (MQ=255)                                                                           taGACAATCTGCCGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCGGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCg                                                                           <  8:119939/149‑1 (MQ=255)                                                                             gACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGaa                                                                         >  1:362226/1‑149 (MQ=255)                                                                             gACAATCTGCGGTGCCGGACCAAAGAAATACCAGATAGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGaa                                                                         <  5:273411/149‑1 (MQ=255)                                                                                      cGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGt                                                                 <  7:302756/148‑1 (MQ=255)                                                                                           ccGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAg                                                            <  8:174846/148‑1 (MQ=255)                                                                                           ccGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGc                                                           >  4:97825/1‑149 (MQ=255)                                                                                             ggACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCg                                                          >  7:129123/1‑148 (MQ=255)                                                                                              gACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCgta                                                        <  1:38007/149‑1 (MQ=255)                                                                                               aCCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCgtag                                                       >  5:353818/1‑149 (MQ=255)                                                                                                ccAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCgtag                                                       >  6:254586/1‑148 (MQ=255)                                                                                                   aaGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATg                                                   <  1:272447/149‑1 (MQ=255)                                                                                                   aaGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATg                                                   >  1:204088/1‑149 (MQ=255)                                                                                                    aGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATgg                                                  >  8:14583/1‑149 (MQ=255)                                                                                                         tACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGgc                                                                                     <  8:211427/109‑1 (MQ=255)                                                                                                         tACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGgc                                                                                     >  7:211427/1‑109 (MQ=255)                                                                                                           ccAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACa                                                             >  8:180747/1‑131 (MQ=255)                                                                                                           ccAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACa                                                             <  7:180747/131‑1 (MQ=255)                                                                                                           ccAGAATGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATATCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCaa                                           <  4:275868/149‑1 (MQ=255)                                                                                                                    cGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTgg                                  >  6:139535/1‑149 (MQ=255)                                                                                                                     gCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTgg                                  >  1:218421/1‑148 (MQ=255)                                                                                                                     gCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGc                                 <  5:301348/149‑1 (MQ=255)                                                                                                                         tGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGCCAGCGTAGTATGGCTGCCCACCGCCCTGGCACg                              >  3:66673/1‑148 (MQ=255)                                                                                                                         tGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACg                              >  3:281962/1‑148 (MQ=255)                                                                                                                             aaGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACa                                                             <  7:125599/113‑1 (MQ=255)                                                                                                                             aaGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACa                                                             >  8:125599/1‑113 (MQ=255)                                                                                                                             aaGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGaaa                           >  7:9756/1‑147 (MQ=255)                                                                                                                             aaGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGaaa                           >  5:66918/1‑147 (MQ=255)                                                                                                                               gagcAATAACCACCACTACCGGCAGAAATTCGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGc                       <  8:231628/147‑1 (MQ=255)                                                                                                                                  cAATAACCACCACTAACGGCACAAATACGGCTGAGATTTTATCAGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGcaca                    <  2:384210/149‑1 (MQ=255)                                                                                                                                      aaCCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTtgctg                                                                                                  <  6:177809/67‑1 (MQ=255)                                                                                                                                             acTACCGGCACAAATACGGCTGCGATTTTATCCGCCCGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg         <  1:177063/149‑1 (MQ=255)                                                                                                                                                 ccGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTa       <  7:14583/147‑1 (MQ=255)                                                                                                                                                   ggCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCAc   >  3:123561/1‑149 (MQ=255)                                                                                                                                                    gCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt  <  6:66918/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCGACATCGGCGTCGCCAGCCCCAGCGCACACGGACAGGCAATAATCAGTACCGTGGTGGCAATCACCAGGGTATAGACAATCTGCGGTGCCGGACCAAAGAAATACCAGATTGCCGCACTGACAAGCGCAATAACCACCACTACCGGCACAAATACGGCTGAGATTTTATCCGCCAGCTGACCGATTTCTGGCTTGCTGCTCTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACT  >  NC_000913/509909‑510204

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