Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,384	0	A	C	52.8%	7.0 / 30.4	36	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (12/5);  new base C (0/19);  total (12/24)
Fisher's exact test for biased strand distribution p-value = 4.94e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.52e-02

GTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCAT  >  NC_000913/3760240‑3760531                                                                                                                                                    |                                                                                                                                                     gTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGC‑GGCCCCCCCAACAACTTTTCGCCCGCGCAA‑CTCACC‑AGa                                                                                                                                                >  1:201579/1‑150 (MQ=255)   ggCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACt                                                                                                                                              <  2:255966/150‑1 (MQ=255)              ccgcGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGCTCCTCTCTTTTTTACCAACAATAATTCCTTCCTGC‑GGCCCGCCTAACAACCTTTCGCCGGAGAAA‑ATCACC‑CGACTGACGCccgccg                                                                                                                                   >  1:619309/1‑150 (MQ=255)                gcgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGATCATCTCTTTTTTACCCACCATAATTCCCGCCTGC‑GGCCCGCCTAACAACTTTTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCaa                                                                                                                                >  1:160048/1‑151 (MQ=255)                        taTTTTCCCCGCGGGTAATCCACGCTTCAGCGGGGGGCTTGG‑CAGGCGGGCGATCACCTCTTTTTTACCAACAAAAATTCCTCCCGGC‑GGCCCGCCTAACAACTTGTCGCCGGAAAAA‑CCCCCC‑AGACTGGCTCCCGCCGCAATCAACTcc                                                                                                                        <  2:500722/150‑1 (MQ=255)                               tCCCCGCTTAATGCCCGCTTCAGCGGG‑GGCTTTGCCAGGCGGCCGATCATCTTTTTTTTCCCAACAATAATCCCGGCCTCC‑GGCCCTCCTAACACCTTCTCCCCGGAGAAA‑CTCCCC‑ATCCTAACGCCCGCCGCAATCAACTCCTGCTGCa                                                                                                                 <  2:632210/151‑1 (MQ=255)                                 cgcgcgTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGGTCCTCTCTTTTTTACCAACAATAATTCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCAtt                                                                                                               >  1:246990/1‑151 (MQ=255)                                   cgcgTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGGTCATCTCTTTTTTTCCAACAATAATTCCTGCCCGC‑GGCCCCCCTAACAACTTGTCGCCTGAGAAA‑ATCACC‑AGACTGACGCCCGCCGCAATCAACTCCTCCGGCATTgg                                                                                                             >  2:620951/1‑151 (MQ=255)                                              cGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGATCATCTCTTTTTTACCCACAATAATTCCTGCCCGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑CGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa                                                                                                  >  1:573953/1‑151 (MQ=255)                                              cGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGC‑GGCCCGCCTAACAACTTTTCGCCGGAGAAA‑CTCACC‑CGACTGACGCCCCCCCCCATCCACTCCTGCGGGCTTGGCCCTTTTGGCa                                                                                                  >  2:333726/1‑151 (MQ=255)                                                                tG‑CGGGGGGGCGATCATCTTTTTTTTCCAAAAAAAAATCCCGGCCGGG‑GCCCCGCCTAAAAATTTGTCGCCGGAGAAA‑CCCACC‑GGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  2:505538/151‑1 (MQ=255)                                                                tG‑CGGGGGGGCGATCATCTTTTTTTACCCAAAAAAATTTCCTGCCTGC‑GCCCCGCCTAAAAACTTGTCCCCGGAAAAC‑CTCCCC‑AAACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  2:173982/151‑1 (MQ=255)                                                                tG‑CGGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGC‑GCCCCGCCTAACACCTTGTCGCCGGAAAAC‑CCCCCC‑AGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  1:411356/151‑1 (MQ=255)                                                                tG‑CAGGGGGGCGATCATCTCTTTTTTACCACAAAAAATTCCGGCCTGG‑GGCCCGCCAAACAACTGGTCGCGGGAGAAA‑CCCCCC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  2:163804/151‑1 (MQ=255)                                                                tG‑CAGGCGGGCGATCATCTTTTTTTTCCCAAAAAAAATCCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAC‑CCCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  2:610957/151‑1 (MQ=255)                                                                tG‑CAGGCGGGCGATCATCTCTTTTTTACCAACAATAATCCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CCCCCC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  2:472144/151‑1 (MQ=255)                                                                 tgaaGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCGGCCGGC‑GCCCCGCCTAACACCTGGTCGCGGGAGAAC‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  1:516811/148‑1 (MQ=255)                                                                 g‑CAGGCGGGCGATCATCTCTTTTTTCCAAACAATAATCCCTGCCGGC‑GCCCCGCCTAACAACTTGTCGCGGGAGAAC‑CTCCCC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  1:83858/150‑1 (MQ=255)                                                                 g‑CAGGCGGGCGACCTTCTCTTTTTTCCCAACAAAATTCCCGGCCGGC‑GCCCCGCCTAACAACTTGCCGCCGGGGAAA‑CTCCCC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  2:131498/150‑1 (MQ=255)                                                                  gcgggggggcgTTCATCTCTTTTTCACCAAAAATAATCCCTGCCTGC‑GGCCCGCCTAACATCTTGTCCCGGGAAAAC‑CCCCCC‑AGACTCACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  1:235919/144‑1 (MQ=255)                                                                  gcgggcggccgGTCATCTCTTTTTTCCCAAAAATAATTCCTGCCGGC‑GGCCCTCCTAACAACTTGTCGCCGGAGAAC‑CTCCCC‑AGCGAAACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  1:97035/147‑1 (MQ=255)                                                                  gaaGGCGGGCTATCATTTTTTTTTTCCAAACAATAATCCCTCCCTGC‑GCCCCGCCTAACAATTTGCCGCGGGAAAAA‑CTCACC‑AGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                                <  1:620753/148‑1 (MQ=255)                                                                       cggccgATCATCTCTTTTTTACCAACAATAATCCCCGCCTGC‑GGCCCGCCTAACAACTTGTCGCGGGAAAAA‑CTCCCC‑AGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGAcc                                                                           <  1:271749/150‑1 (MQ=255)                                                                             atcatcTCTTTTTTACCAACAATAATTCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAg                                                                     >  1:463422/1‑150 (MQ=255)                                                                              ccatcTCTTTTTTCCCAACAATAATTCCGGCCGGC‑GCCCCGCCTAACAACTTGTCGCGGGAGAAA‑CTCCCC‑AGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGcc                                                                   <  1:527619/150‑1 (MQ=255)                                                                                ttctcTTTTTTCCAAAAAATAATCCCTCCCTGC‑GGCCCCCCTAACAACTTGTCGCCGGAGAAA‑CTCCCC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCAc                                                                 <  2:59769/150‑1 (MQ=255)                                                                                        ttacccAAAAATAA‑TCCCGCCCGCGGGCCCGCCTAACATCTTGTCGCGGGAAAAC‑CTCCCC‑AGCCGAGCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTaaa                                                         <  1:179518/147‑1 (MQ=255)                                                                                          tACAAACAAAATTCCCTGCCCGC‑GCCCCG‑CTAACCACCTTACCCCGGGGAAACCCCCCC‑ACACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGCCCAGCGAGCCACTGCCTAAATc                                                       <  2:159992/151‑1 (MQ=255)                                                                                                caaTAATTCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCa                                                  >  1:102115/1‑150 (MQ=255)                                                                                                caaTAATTCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACCGGCTAAGATCGACCCGCGAGCCACTGCCTAAATCAGTCa                                                  >  2:635443/1‑150 (MQ=255)                                                                                                        cgtgccggc‑GCCCCGCCAAACACTTGGCCGCGGGAAAAA‑CCCCCC‑AGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGaa                                         <  1:541739/149‑1 (MQ=255)                                                                                                        cctgccggC‑GGCCCCCCAAACAACTTGTC‑CCGGAAAAA‑CCCCCCAAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGaa                                         <  1:272523/151‑1 (MQ=255)                                                                                                        ccggccggc‑GGCCCCCCTAACACCTTGTCGCGGGAGAAC‑CACCCC‑AGACGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCGCTACGGGaa                                         <  2:609881/148‑1 (MQ=255)                                                                                                                   ccccGCCAAACACCTTGTCCCGGGAGAAA‑CCCCCA‑AAACTCACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGctc                               <  2:160088/150‑1 (MQ=255)                                                                                                                        ccTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTg                           >  1:276265/1‑150 (MQ=255)                                                                                                                                                 cTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCAt  >  2:315085/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTG‑CAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGC‑GGCCCGCCTAACAACTTGTCGCCGGAGAAA‑CTCACC‑AGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCAT  >  NC_000913/3760240‑3760531

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