Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	6,125,257	0	T	C	21.6%	89.6 / 11.5	60	T418A (ACC→GCC)	copA_4	Copper resistance protein A
Reads supporting (aligned to +/- strand):  ref base T (18/22);  new base C (0/11);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 4.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATC  >  ACB122/6125122‑6125392                                                                                                                                        |                                                                                                                                        cttcgcccATCCCAGCCTTCTTGCTGTGGTCCATTCCGGCCATGCCGCCTTGGTCCATCGCAGCCATTTTGCTGTGATCCATTGCGGCCATGCGGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTc                                                                                                                                        <  1:193591/131‑1 (MQ=14) cgtGGTCCATCCCATCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGCTCCATGCCGGTCATGCTGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTc                                                                                                                                        <  4:227240/135‑1 (MQ=14)    ggTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCGGCCGTGGTCCATCTCAGCCATTTTGAGGTGATCCATGCCGGCCATGATGCCATGGTCCATCCCGGCCATGATGCGGTGGTCCATACCGGTCATg                                                                                                                                     >  7:332617/1‑137 (MQ=14)      gccATCCCAGCCATCTTGCTGTGGTCCATTCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGtt                                                                                                                                   <  3:39081/136‑1 (MQ=255)       ccatcccatcCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTg                                                                                                                                  <  6:397760/137‑1 (MQ=255)       ccATCCCCTCCATCTTGCTTTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCGCAGCCATTTTGCTGTGCTCCATGCCGGCCATGCTGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTg                                                                                                                                  <  2:159564/137‑1 (MQ=17)           cccAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGt                                                                                                                              <  6:74617/137‑1 (MQ=255)             cAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTg                                                                                                                                  >  3:88366/1‑131 (MQ=255)             cAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTg                                                                                                                                  <  4:88366/131‑1 (MQ=255)             cAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCCGCCATGATCCCATGGTCCATCCCCGCCATGATTCTGTGGTGAATGCCGGACATGATGCGGTgg                                                                                                                            >  1:389721/1‑137 (MQ=17)             cAGCCATCTTGCTGTGGTCCATGCCGGCCATGCAGCCGTGGTCCATCACAGCCATTTTGCTGTGGTCCATTACGGCCATGTTGCCCTGGTCCATCAAAGCCATGTTGAGGTGGACGATGACGGACATGTTgaagtgt                                                                                                                            >  6:132843/1‑131 (MQ=2)                  aTCTTGCTTTGGTCCATCCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCAt                                                                                                                       <  8:249537/137‑1 (MQ=255)                      tGCTGTGGTCCATGCCGGCCATGGAGCCGTGGTCCATCACAGCCATTATGCTGTGATCCATACCGGCCATGCTGCCATGGTCCATACCAGCCAGTTTGCGTTGGTCCACGCCCGACATCATGCCTTCGTCAAACCCa                                                                                                                   >  6:405361/1‑137 (MQ=9)                       gCTGTGGTCCATCCACGCCATCCCGCCATGGTCCATCACAGCCATTTTGCTGTGGTCCATCCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGATGCCGTGGTCCATCCCAg                                                                                                                  <  2:384005/137‑1 (MQ=14)                         tgtgGTCCATGCCGGCCATCCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGcc                                                                                                                <  3:383566/137‑1 (MQ=255)                              cccATGCCGGCCATGCCGACGTGGTCCATCTGAGCCATTTTTCTGTGATACATGCCTGCCATGCTGCCATGGTCCATCGGAGCCATTTAGCTGTGGTCCATGCCGGCCATGTTGCCTTGGTCCATCCCAGCCAtttt                                                                                                           <  7:321905/136‑1 (MQ=14)                               ccATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGTTCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTg                                                                                                          <  3:335993/137‑1 (MQ=255)                               ccATCCCGGCCATGCCGCCGTTGTCCATCACAGCCATATTGCTGTGTTCCATCCCGGCCATGCGGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGGCATGTTGCCGTGGTCCATCCCAGCCATTTTg                                                                                                          <  8:331320/137‑1 (MQ=14)                                cATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTg                                                                                                                                  >  7:376320/1‑112 (MQ=255)                                cATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTg                                                                                                                                  <  8:376320/112‑1 (MQ=255)                                cATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATGTTGc                                                                                                         >  4:199594/1‑137 (MQ=255)                                          aTGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCAt                                                                                               >  1:55490/1‑137 (MQ=255)                                          aTGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATACCAGCCATTTTGCTGTGGTCCAt                                                                                               >  4:195140/1‑137 (MQ=255)                                          aTGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCAGGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCAt                                                                                               >  1:40300/1‑137 (MQ=255)                                               gccgTGTTCCATCTCAGCCATTTTGCTGTGCTCCATGCCGTCCATGCTGACCTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCgg                                                                                          <  1:111774/137‑1 (MQ=14)                                                 cgTGGTCCATCTCAGCCATTTTGCTGTGCTCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGcc                                                                                        <  2:55490/137‑1 (MQ=255)                                                 cgTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCCTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGcc                                                                                        <  8:34811/137‑1 (MQ=255)                                                   tGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                           <  5:67298/116‑1 (MQ=255)                                                   tGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                                                                           >  6:67298/1‑116 (MQ=255)                                                        cATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTATGCTGTGGTCCATCCAGGCCATGTTGc                                                                                 >  4:404043/1‑137 (MQ=255)                                                                  aTTTTGCTGTGCTCCATGCCGGCCATGCGGCCTTCGTCCATCCCAGCCATGATGCTTTTGTCCATGCCGGCCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCAt                                                                       <  3:286261/137‑1 (MQ=11)                                                                         tgtgGTCCATCCCGGACATGATGCCGTGGTCCATCCGAGCCATTTTGCTGTGGTCCATTACGGTCATGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGTCGGCCATGTTGCCGTGGTCCATCCCAGcc                                                                <  2:42642/137‑1 (MQ=11)                                                                              tcgaTGCCGGCCATTCGGCCATGCTCCATCCCAGCCATGTTGCTTTGGTCCATGCCGGCCATGGGGCCTTGGTCCATCCCAGCCATGTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                           <  7:436167/134‑1 (MQ=11)                                                                              tCCATGCCGGCCATGCGGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                           <  6:396780/137‑1 (MQ=255)                                                                              gcaatcccGCCCATGCTGCCTTGGTCCATCGCAGCCATGATGCTTTGGTCCATGCCGGACAGGCTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAtttt                                                           <  8:336546/131‑1 (MQ=11)                                                                              cccATGCCGGCCATGCTGCCTTGGTCCATCCCAGCCATCTTGCTCTGGTCCATTCCGGCCATGCAGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGATGCCGTGGTCCATCCCAGCCAtttt                                                           <  6:24379/136‑1 (MQ=11)                                                                                cATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCAGTTTgg                                                         >  3:271487/1‑136 (MQ=255)                                                                                cATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGACATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGc                                                         >  4:409695/1‑137 (MQ=255)                                                                                cATGCCGGCCATGCTGCCATGGTCCATCCCAGACATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCATTCCCAGCCATTTTGCTGTGGTCCATGCCGGAAATGTTGCCGTGGTCCATACCGGACATTTTgg                                                         >  1:163107/1‑136 (MQ=14)                                                                                    ccGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCGATGCCCGCCATTTTGCtgtg                                                     >  1:325088/1‑137 (MQ=255)                                                                                               gccctcgTCCATCCCCGCCATATTGCTGTGGTCCATGACTGACATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCgg                                          <  7:228677/131‑1 (MQ=14)                                                                                                 cttGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGGCATGGTGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGcc                                        <  5:435454/135‑1 (MQ=255)                                                                                                 cATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGcc                                        <  3:406458/137‑1 (MQ=255)                                                                                                      tCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGtt                                   <  3:409695/137‑1 (MQ=255)                                                                                                       ccagcacAGCCATAGTGCCGGGGTCCATGCCGGCCATGCCGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTAGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTg                                  <  5:272160/131‑1 (MQ=11)                                                                                                        cATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAg                                                                                                                  >  1:46665/1‑56 (MQ=35)                                                                                                        cATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAg                                                                                                                  <  2:46665/56‑1 (MQ=35)                                                                                                             cAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTTCTGTGGTCCATACCGGCCATGTTGCAGTgg                            >  3:159912/1‑137 (MQ=255)                                                                                                             cAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCAGGTTGCCGTgg                            >  1:32640/1‑137 (MQ=255)                                                                                                             cAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTACATACCGGCCATGTAGCAGTgg                            >  1:57720/1‑137 (MQ=255)                                                                                                                  aTTTTGCTGTGGTCCATGCCGGGCATGTTGCCGTGGTCCATCCCAGCCATTTTGATGTGGTCCATGCCGCGCATGTTGACTTGGTGCATCCCAGCCATATTGCATTCGTCCAGTCCGGCAATGATGACTTggagaaa                       >  6:238671/1‑132 (MQ=2)                                                                                                                    tttGCTGTTGTCCATGCCGCGCCTTTTGCCGTGGTCCATCCCAGCCCTTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGc                     <  5:132843/137‑1 (MQ=14)                                                                                                                              cccctGCCGGCCATGCCGCCGTGGCCCCTCCCCGCCCTCTTGATTTGGACCCTGCCGGCCATGAGGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATCCCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgct           <  8:332617/133‑1 (MQ=11)                                                                                                                               ccctcccGGTCATCTTGCCCTGCTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgctg          <  2:325088/132‑1 (MQ=17)                                                                                                                               ccATGCGGGTCATGGTGCCGTGGTCCATGCCCTCCTTTGTGCTTTGGTCCATGGCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgctg          <  5:282201/137‑1 (MQ=17)                                                                                                                               ccATCCCCGCCATGTGGCCGTGGTCCATCCCTGCCATGCTGCCGTGGTCCATCCCGGCCATGTTGCGGTGGTCCATCCCCGCCATGATGCCGTGGTCCATCCCTGCCATGTTGCCGTGGTCCATGCCGGCCAtgctg          <  8:125461/137‑1 (MQ=14)                                                                                                                                cATGCCGGTCATGTTGCCGTGGTCCATCCCACCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgctgc         >  4:209441/1‑137 (MQ=255)                                                                                                                                cATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATCCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgctgc         <  1:38367/137‑1 (MQ=255)                                                                                                                                aatcgcGGCCATGTTGCCGGGGTCCATCCCTTCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCAtgctgc         <  1:83736/132‑1 (MQ=14)                                                                                                                                       gTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCATGATc  <  3:404043/137‑1 (MQ=255)                                                                                                                                       ctCATGTTGCCGGGGGCCATCCCCGCCATTTTGCGGTGGGCCATGCCGGCCATGTTGCCTTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATc  <  7:308701/136‑1 (MQ=17)                                                                                                                                        |                                                                                                                                        CTTGGTCCATCCCAGCCATCTTGCTGTGGTCCATGCCGGCCATGCCGCCGTGGTCCATCTCAGCCATTTTGCTGTGATCCATGCCGGCCATGCTGCCATGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGTCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATGCCGGCCATGTTGCCGTGGTCCATCCCAGCCATTTTGCTGTGGTCCATACCGGCCATGTTGCCGTGGTCCATGCCGGCCATGCTGCCGTGATC  >  ACB122/6125122‑6125392

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