Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,086	0	A	C	48.1%	83.8 / 11.5	52	Y580D (TAC→GAC)	ggt	gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (26/1);  new base C (0/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.09e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTA  >  NC_000913/3584937‑3585228                                                                                                                                                       |                                                                                                                                               ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAc                                                                                                                                               >  2:751475/1‑151 (MQ=255)       aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                         >  2:194721/1‑151 (MQ=255)           aGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt                                                                                                                                      >  2:913224/1‑150 (MQ=255)           aGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt                                                                                                                                      >  1:754105/1‑150 (MQ=255)                      aaaTTCTTTCTGGGGGAACGAAAAACCCGCGGAAAGAATATCTACAGGATTACA‑AAAAAAAAAAGGCTCCCTTCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTACCCCGCCGTTAAATCATCCAcc                                                                                                                          <  2:83778/151‑1 (MQ=255)                         ttcattcTCGGGGAACGGAAAACCCGCTGGAAAAATATCTACAGAATTACA‑AAAAAAAAAGGGCTCCTTTCGGCTTGCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGGCCGCTACCCTTATTCCCCCGCCGTTAAATCATCCACCGa                                                                                                                        <  1:53749/150‑1 (MQ=255)                             ttcTCGGGGAACGAAAAACCCGCGGAAAGAATATCTACAGAATTACA‑AAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAACGCCCCTTTCCCCACGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAgc                                                                                                                      <  1:250977/148‑1 (MQ=255)                                tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                                                >  1:567152/1‑151 (MQ=255)                                   gggACCGGAAAAACCGCTGAAAGAATTTCTAAAGAATTAAG‑AAAAAAAAAGGGCCCCCTCCGTCTTCCCCGGAAAAAATAGCCCTCTCCCCCCGAGGAGGCCCCCTAACCTTGGTCCCCCGCCGTTAAATCATCCACCGAGcgcg                                                                                                                   <  2:810109/145‑1 (MQ=255)                                            aaacCCGCTGAAAAAATATCTACAGGATTACA‑AAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATACCCCTCTTCCCACGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc                                                                                                       <  2:469310/148‑1 (MQ=255)                                                    tGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                             >  1:660408/1‑150 (MQ=255)                                                     ggAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                           >  1:411317/1‑151 (MQ=255)                                                     ggAAAATTAGCTACA‑GATTAAGAAAAAAAAAAAGGCTCCCTTCGGCCTGCCCTGACAAAAAAGCCCTCCCCCCACAAAGAGGGCGGTTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                            <  1:519638/150‑1 (MQ=255)                                                       aaaaaTTTCTACAGAATTAAA‑AAAAAAAAAAGGCCCCCTTCGGCTTCCCCTAAAAAAATGCCCCTCTCCCCACGAAGGGGGCCCTTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                            <  1:712683/145‑1 (MQ=255)                                                       aaaaaTATCTACAGAATTAAG‑AAAAAAAAAGGGCTCCCTTCGCCTCGCCCTGAAAAAATCCCCCTCTCCCCACGAAGGGGCCCCTTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                              <  1:851778/143‑1 (MQ=255)                                                          aaTATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                      >  2:698856/1‑151 (MQ=255)                                                           atatCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                     >  2:599803/1‑151 (MQ=255)                                                                 aCAGGATTACG‑AAAAAAAAAGTTCTCCCTTCGCCTGCCCCTCAAAAAAGACCCCTCTCCCCAGCAGGGGGGCCCCAACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                     <  1:120355/145‑1 (MQ=255)                                                                  cAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                              >  1:393296/1‑151 (MQ=255)                                                                  cAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                              >  2:100752/1‑151 (MQ=255)                                                                   aGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                             >  2:61148/1‑151 (MQ=255)                                                                    ggATTACA‑AAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                  <  2:299514/145‑1 (MQ=255)                                                                        ttaaaaaaaaaaaaaGCCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACAAAGGGGGCCGCTACCTTTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                               <  1:6625/140‑1 (MQ=255)                                                                          acaaaaaaaaaaaGGCTCCCTTCGGCTTCCCCTGACAAAATCCCCCTCTCCCCACGAAGGGGCCCTCTACCCTTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                       <  1:104509/148‑1 (MQ=255)                                                                             aaaaaaaaaaGGCTCCCTTCGGCTTGCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCCCTACCTTTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                     <  1:198581/150‑1 (MQ=255)                                                                             aaaaaaaaaaGGCTCCCTTCGGCTTCCCCGGACAAAACACCCCTCTCCCCACAAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                                            <  2:586152/143‑1 (MQ=255)                                                                             aaaaaaaaaaGGCTCCCTCCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                     <  2:841482/150‑1 (MQ=255)                                                                               aaaaaaaaGCCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCCCCCCACGAAGAGGCCCCCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                  <  1:955953/149‑1 (MQ=255)                                                                                tAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                 >  1:218874/1‑151 (MQ=255)                                                                                    aaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                             >  2:128205/1‑151 (MQ=255)                                                                                    aaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                             >  1:45603/1‑151 (MQ=255)                                                                                      aGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                           >  1:412180/1‑151 (MQ=255)                                                                                         ctcccTTGGGCTGCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                                         <  2:90942/147‑1 (MQ=255)                                                                                         ctcccTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                           <  1:985524/145‑1 (MQ=255)                                                                                           aCCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                      >  2:176585/1‑151 (MQ=255)                                                                                             cTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                    >  2:656788/1‑151 (MQ=255)                                                                                                  gCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                               >  1:738175/1‑151 (MQ=255)                                                                                                    ttccccTAAAAAAATACCCCTCTCCCCGCGAGGGGGCCCGCAAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                 <  2:200026/144‑1 (MQ=255)                                                                                                    ttGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCTGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                             >  2:26260/1‑151 (MQ=255)                                                                                                      gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                            >  2:948321/1‑150 (MQ=255)                                                                                                         cTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                                        >  1:156856/1‑151 (MQ=255)                                                                                                         cTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgc                                         >  1:915257/1‑150 (MQ=255)                                                                                                                 aTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                 >  1:551372/1‑150 (MQ=255)                                                                                                                             ccACGAAGAGGCCCCCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                    <  2:146780/151‑1 (MQ=255)                                                                                                                               aCGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                  >  2:202119/1‑151 (MQ=255)                                                                                                                                  aaGGGGCCCCCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                 <  2:606699/149‑1 (MQ=255)                                                                                                                                  aaGAGGCCCCCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                 <  2:524459/149‑1 (MQ=255)                                                                                                                                     agGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc              >  2:287239/1‑149 (MQ=255)                                                                                                                                      ggcccGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagca             <  1:232228/146‑1 (MQ=255)                                                                                                                                          cgtaacccTTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa       <  2:167106/145‑1 (MQ=255)                                                                                                                                          cccTAACCTAAATCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa       <  2:419741/149‑1 (MQ=255)                                                                                                                                                  ttAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTa  <  2:913305/148‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACG‑AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTA  >  NC_000913/3584937‑3585228

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