Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	599,576	0	T	G	42.0%	80.4 / 17.7	50	V288G (GTG→GGG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (0/29);  new base G (20/1);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 4.46e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGC  >  NC_000913/599433‑599712                                                                                                                                                |                                                                                                                                          cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATc                                                                                                                                    <  1:176037/151‑1 (MQ=255)  gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATcc                                                                                                                                   <  2:180566/151‑1 (MQ=255)         ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑c                                                                                                                           >  2:803049/1‑151 (MQ=255)                   cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCaaaaa                                                                                                                  <  1:659586/150‑1 (MQ=255)                     ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCCCGGG‑CCGGGaaaaa                                                                                                                  >  1:694768/1‑148 (MQ=255)                      tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGc                                                                                                               <  1:413145/150‑1 (MQ=255)                         aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGGT‑CCGGCaaaaaaccc                                                                                                              >  2:391484/1‑144 (MQ=255)                           aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCAAAAACGccc                                                                                                             >  2:965713/1‑147 (MQ=255)                            aGCCAGTGAAAATGGCGTGCCCGTTTATTTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGCCCTGCGAC‑CCGGCcaaaaataccg                                                                                                            >  2:170550/1‑141 (MQ=255)                                aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGa                                                                                                    <  2:71475/151‑1 (MQ=255)                                   gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATc                                                                                                  <  1:926370/150‑1 (MQ=255)                                   gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATc                                                                                                  <  2:1006946/150‑1 (MQ=255)                                                 cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                                                                   <  2:809362/151‑1 (MQ=255)                                                 cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGT‑CCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGGAgg                                                                                     >  2:122384/1‑149 (MQ=255)                                                    ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                 <  1:773360/150‑1 (MQ=255)                                                         tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTGCGGGCAAAAAC‑CCCGCGAAGGGGTCGCCGCCGCGGAGGGCAAAATgg                                                                            >  1:116881/1‑150 (MQ=255)                                                          gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                            <  1:794253/149‑1 (MQ=255)                                                          gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAT‑CCGGCAAAAAAGCCCCGGAAGTGATCGCCGCCGTGAAGGACAAACAgg                                                                            >  1:672351/1‑149 (MQ=255)                                                                tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                     >  2:146348/1‑150 (MQ=255)                                                                 gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                    <  2:109882/150‑1 (MQ=255)                                                                    gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                 <  2:978665/150‑1 (MQ=255)                                                                     cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGCGAAGGACAAACTGGAAACGCTGaaa                                                                 >  1:9048/1‑149 (MQ=255)                                                                      gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                              <  1:199315/151‑1 (MQ=255)                                                                            ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                         <  2:117605/150‑1 (MQ=255)                                                                            ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                        <  1:888594/151‑1 (MQ=255)                                                                              aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGGT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGGAGGGCAAACAGGAAACGCTGAAAAGTAGTCTGcc                                                      >  1:550761/1‑151 (MQ=255)                                                                                   ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CCGGCAAAAAAGCCCGAGAAGGGAACGCCCCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                    >  1:901293/1‑148 (MQ=255)                                                                                   ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGGCTGCCgg                                                    >  2:246049/1‑148 (MQ=255)                                                                                         gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                           <  2:559071/151‑1 (MQ=255)                                                                                           gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                         <  2:23924/151‑1 (MQ=255)                                                                                           gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                         <  1:504992/151‑1 (MQ=255)                                                                                              gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                       <  2:526247/150‑1 (MQ=255)                                                                                               cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                     <  2:543411/151‑1 (MQ=255)                                                                                                gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGAT‑CCGGCAAAAAAACCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGATa                                      >  2:214711/1‑149 (MQ=255)                                                                                                 ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGtt                                   <  1:717750/151‑1 (MQ=255)                                                                                                          tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGAT‑CCGGCAAAAACACCCGGGAAGTGATCGCCGCCGTGGAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTACAACATa                           >  1:487907/1‑150 (MQ=255)                                                                                                          tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGGCAAACTGGAAAAGCTGGAAAGTAGTCTGCCGGAAGGCGGGGGGATAGTTACAACAAAc                          >  2:1033492/1‑151 (MQ=255)                                                                                                             cGAACTAAACGGCGAAGGGGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                       <  1:86606/151‑1 (MQ=255)                                                                                                              gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                      <  1:55228/151‑1 (MQ=255)                                                                                                              gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                      <  1:83033/151‑1 (MQ=255)                                                                                                              gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGGCAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGGGGGGAGATAGTTACAACATACGATc                      >  2:348426/1‑151 (MQ=255)                                                                                                              gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGAAAAGTAGGCTGCCGGAAGGCGGGGGGGTAGTTACAACATACg                         >  1:174299/1‑148 (MQ=255)                                                                                                                   aaccGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcag                  <  1:442129/147‑1 (MQ=255)                                                                                                                          gaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGAT‑CCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTACAACATACGATCGCAGCCAGCTc           >  2:139025/1‑150 (MQ=255)                                                                                                                           aaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGCAGCCAGCTCAt         >  2:100933/1‑151 (MQ=255)                                                                                                                              gcgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGa      <  2:670445/151‑1 (MQ=255)                                                                                                                                gaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc    <  2:694983/151‑1 (MQ=255)                                                                                                                                gaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc     <  1:232477/150‑1 (MQ=255)                                                                                                                                 aagTGGCGGGCGGGGGGGTGATCCGGCGAT‑CCGG‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc    <  1:391342/149‑1 (MQ=255)                                                                                                                                  agTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc  <  1:46913/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                          CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGC  >  NC_000913/599433‑599712

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