Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,086	0	A	C	44.2%	71.2 / 10.4	43	Y580D (TAC→GAC)	ggt	gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (20/4);  new base C (0/19);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.65e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCG  >  NC_000913/3584936‑3585232                                                                                                                                                         |                                                                                                                                                   attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAACAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTa                                                                                                                                                    <  1:269949/151‑1 (MQ=255)    attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTAcc                                                                                                                                                  <  2:366360/150‑1 (MQ=255)      tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGATTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCccg                                                                                                                                               >  2:613736/1‑151 (MQ=255)        aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCccgcc                                                                                                                                             <  1:85808/151‑1 (MQ=255)        aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCccgcc                                                                                                                                             >  1:207321/1‑151 (MQ=255)          aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCccgccg                                                                                                                                            <  1:460111/150‑1 (MQ=255)           aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGtt                                                                                                                                          >  1:542677/1‑151 (MQ=255)            aGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGtt                                                                                                                                          >  1:256259/1‑150 (MQ=255)            aGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCACTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGtt                                                                                                                                          >  2:267685/1‑150 (MQ=255)                  tggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAatcat                                                                                                                                   >  1:102046/1‑151 (MQ=255)                  tggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAatcat                                                                                                                                   >  1:438243/1‑151 (MQ=255)                          ttctttctTGGGAAACGGAAAACCCGCTGAAAAAATATCTACAGGATAACGAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAAGACCCCTCTCCCCACAAAGGGGCCCGCAA‑CCGTA‑AGTCCCCCGCCGTTAAATCATCCACCGa                                                                                                                            <  2:395981/150‑1 (MQ=255)                           tcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                          >  1:136260/1‑151 (MQ=255)                               tctcTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                      >  1:469586/1‑151 (MQ=255)                                    gggACCGACAAACCCGGGGGAAAAATATCCAAAGGATAACAAAAAAAAAAAGGCTCCCTCCGCTTTCCCCTCACAAAACCCCCCTCTCCCAGCGAGGGGGCCCGTTA‑ACCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                 <  2:759431/151‑1 (MQ=255)                                    ggaaaCGGAAAAACCGGTGGAAGAATATCTACAGGATTACAAAAAAAAAAAGCCTCCCTTCGGTTTCCCCTAAAAAAATGCCCCTCTCCCCACGAAGGGGCCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                  <  1:298879/147‑1 (MQ=255)                                             aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCg                                                                                                         >  2:346739/1‑150 (MQ=255)                                                    cTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCATCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                 >  2:373667/1‑151 (MQ=255)                                                               tCTACAGGATTAAAAAAAAAAAAAGGCTCCCTCCGGTTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTA‑CCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                        <  1:106125/149‑1 (MQ=255)                                                                  aCAGAATTACAAAAAAAAAAGGTCTCCCTTCGGCTTCCCCGAAAAAAATGCCCCTCTCCCCACGAGGAGGCCCGCTA‑ACCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                                   <  1:228939/151‑1 (MQ=255)                                                                     ggATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGACAAAATCGCCCTCTCCCCACGAAGGGGGCCGCTA‑ACCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                  <  1:594341/149‑1 (MQ=255)                                                                        ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCTGGCCCAACCATa                                                                              >  2:523142/1‑150 (MQ=255)                                                                                  aaaaaGGCTCCCTTCGGCTCCCCCTAAAAAATCACCCCTC‑CCCCAACAAGGGGGCCCCCACCCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                   <  2:286117/151‑1 (MQ=255)                                                                                   aaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                   >  2:625791/1‑150 (MQ=255)                                                                                      aGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                               >  1:79153/1‑151 (MQ=255)                                                                                           cccTTCGGCTTGCCTAAAAAAAAGACCCCCCTCCCCCCGAAGGGGGCCGCTA‑ACCTT‑AGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                               <  2:81005/145‑1 (MQ=255)                                                                                           cccTTCGGCTTGCCCTGACAAAATACCCCTCTCCCCACGAAGGGGCCCGCAA‑CCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                          <  1:1210/150‑1 (MQ=255)                                                                                           cccTCCGCCTTCCCCTAAAAAAATACCCCTCTCCCCACAAGGAGGCCCGCTA‑CCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCa                                                           <  1:792058/149‑1 (MQ=255)                                                                                               tCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGcc                                                      >  1:379440/1‑151 (MQ=255)                                                                                                 ggCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                    >  1:728423/1‑151 (MQ=255)                                                                                                    ttGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                 >  1:290874/1‑151 (MQ=255)                                                                                                    ttGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                 >  1:58128/1‑151 (MQ=255)                                                                                                     tccccTGAAAAAATGCCCCTCTCCCCACGAAGGGGCCCGCTA‑CCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                    <  2:835597/145‑1 (MQ=255)                                                                                                       cccTGAAAAAACACCCCCCTCCCCACAAAGGGGCCCGCTA‑CCTTT‑AGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                       <  2:438137/142‑1 (MQ=255)                                                                                                              aaaaacgcctcctcccccCGAAGAGGCCCGC‑A‑ACCTTGATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgc                                             <  2:788246/131‑1 (MQ=255)                                                                                                                            ccgaCGAAGGGGGCAGCTT‑ACCTT‑CGCCCACCGCCGTTAAATCATCCACCGAGCGCGGGTCGTATGCGCCGTACAACTCACCGGCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGGCCgttt                           <  2:474411/146‑1 (MQ=255)                                                                                                                            cccACGAAGGGGCCCGCTA‑CCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                         <  1:346836/151‑1 (MQ=255)                                                                                                                                 gAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGGCGTTAAATCATCCACCGAGCACGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCCTAATGCTTTGTGTACTGCCCATCGCCTCTTTCATCGCCACTTTCTGACCTTTTGCTTcc                     >  2:585642/1‑150 (MQ=255)                                                                                                                                      gccccGCTA‑CCCTT‑AGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCca                    <  2:593599/143‑1 (MQ=255)                                                                                                                                         ccccTA‑CCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg            <  1:816275/148‑1 (MQ=255)                                                                                                                                           gCTA‑CCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa           <  1:669270/150‑1 (MQ=255)                                                                                                                                              a‑aCCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt       >  2:365171/1‑151 (MQ=255)                                                                                                                                               accttt‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg        <  1:272298/146‑1 (MQ=255)                                                                                                                                                      aGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCg  <  1:625944/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                   ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTA‑ACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCG  >  NC_000913/3584936‑3585232

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