Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,384	0	A	C	44.9%	28.2 / 14.9	49	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (14/13);  new base C (2/20);  total (16/33)
Fisher's exact test for biased strand distribution p-value = 1.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCT  >  NC_000913/3760236‑3760527                                                                                                                                                       |                                                                                                                                                cAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGGGAAACT‑CAc                                                                                                                                                >  1:436976/1‑150 (MQ=255) cAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CAcc                                                                                                                                               <  2:194337/151‑1 (MQ=255)     gTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATTATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCCGAGAAAAc‑caccaaaa                                                                                                                                           >  1:324122/1‑148 (MQ=255)      tGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGAc                                                                                                                                           <  2:1100065/150‑1 (MQ=255)                    gcgAGGGTCCTTTTATCCGAGCGTAATGCACGCTTCAGCGGGTTGCCTTGCAGGCGGGGGATCACCCTTTTTTTACC‑CCCACTAATTCCTGCCTTGGGGCCGCCCAGCACACTGTTCCCGGAGAAACT‑CACCCCACCGGCGCCCGCCGCaa                                                                                                                            >  2:49445/1‑151 (MQ=255)                            ctttttATCCGCGCTTAATGCCCGCTTCAGCGGGTGGCTTTGCAGGGGGGCGTTCATCTTTTTTTTACC‑AACAATAATCCCTCCCTGCGCCCCGCCTAACACCTTGTCGCCGGAAAACCT‑CACCAGACTGCCGCCCGCCGCAATCAACTcc                                                                                                                    <  1:140604/149‑1 (MQ=255)                            ctttttATCCGCGCGTAATGCACGCTTCAGCGGTTGGCTTTGCAGCCGGGCGATCATCTCTTTTTCCCC‑AACAATATTCCCTCCCTGCGGCCCGCCTAACACCTTGTCGCGGGAAAAACT‑CCCCAGACGGACGCCCGCCGCAATCAACTcc                                                                                                                    <  2:290388/149‑1 (MQ=255)                            cATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGCTCCTCTCTTTTTTACC‑CACCATAATTCCTTCCTGCCGGCCGCCTAACAACTTTTCGCCGGAGAAACT‑CCCCAGCCTGACGCCCGCCGCAATCAACTc                                                                                                                     >  2:609565/1‑150 (MQ=255)                                   cccGCGGGTAATCCCCGCTTCGGCGGGTGGTTTGGCGGGCGGGCGATCATCTCTTTTTACCC‑AACAAAAATCCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACC‑CCCAAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCa                                                                                                             <  1:1178898/150‑1 (MQ=255)                                             aTGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCCCCTAACCACTTTTCGCCGGAGAAAAT‑CACCCGACTGACGCCCCCCGCAATCCACTCCCGCGGCATTGGCTCttt                                                                                                   >  1:56191/1‑151 (MQ=255)                                               gCACGTTTCAGGGGGGGGTTTTGCGGGGGGGCGATCATCTTTTTTTTCCC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCttt                                                                                                   <  1:899240/149‑1 (MQ=255)                                                    cTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACT‑CCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCcaa                                                                                            >  1:374556/1‑151 (MQ=255)                                                              tGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCCGACTGACGCCCCCCGCAATCAACCCCTGCGGCCTTGGCTCTTTCGGCAAACCGTACTgg                                                                                   >  1:27027/1‑150 (MQ=255)                                                                gCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACCACTTGTCGCCCGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACCCCTGCGGCATTTGCTCTTTCGGCAAACCGTACTTGCTa                                                                                >  2:1077561/1‑151 (MQ=255)                                                                gCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATCCCTGCCTGCGGCCCGCCAAAAAACTGGTCGCGGGAGAACCT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCt                                                                                 <  2:109796/150‑1 (MQ=255)                                                                gCTTGGCAGGGGGGCGATCATCTCTTTTTTCCC‑AACAATATTCCCTGCCGGCGCCCCGCAAAACAACTTGCCGCCGGAGAACCT‑CACCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACTGTACTGGCt                                                                                 <  2:1189932/150‑1 (MQ=255)                                                                gCTTGGCAGGCGGGCGATCATCTTTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCCCCTAACAAGTTGTGGCCGGAGAAACT‑CCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCt                                                                                 <  1:599924/150‑1 (MQ=255)                                                                 cTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCCTCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTCGCTCCTTCCGCAAAACGTACTGGCCaa                                                                               >  2:380930/1‑151 (MQ=255)                                                                 cTTTGCAGGCGGGCGATCATCCCTTTTTTACC‑AACAAAAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAAACAACTCCTGCGGCATTGGCTCTTTCCGCAAAACGGACTGGCTaa                                                                               >  2:97902/1‑151 (MQ=255)                                                                   ttGCAGGCGGGCGATCATCTCTTTTTTACC‑CACAATAATTCCTGCCTGCGGCCCGCCCAACAACTTGTCGCCCGGGAAACT‑CACCAGACTGACGCCCGCCGCAAACAACTCCTGCGGCATTTGCTCGTTCGGCAAACCGTACTGGCTCAGa                                                                             >  2:1050867/1‑151 (MQ=255)                                                                    tgagggggggcgATCATCT‑TTTTTTACCAAAAAATAATCCCTGCCGGCGGCCCGCCTAAAAACTTGTCGCCGGAGAAACT‑CACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                            <  1:609792/144‑1 (MQ=255)                                                                    tGCGGGCGGGCGATCATTTTTTTTTTACC‑AAAAATAATCCCTGCCTGGGGCCCGCCTAACAACTTTTCGCCGGAGAAACC‑CACCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                            <  1:972958/151‑1 (MQ=255)                                                                    tGCAGGGGGGGGATCATCTTTTTTTTCCC‑AACAAAAATCCCTGCCTGCGGCCCCCCTAACAACTTGCCGCGGGAAAAACC‑CCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                            <  1:300386/151‑1 (MQ=255)                                                                     gcgggcgggcgATCATCTTTTTTTTACC‑AAAAATAATCCCTCCCTGCGCCCCGCCTAACAACTTGTCGCCGGAAAAACC‑CCCCAGACTACCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                            <  2:187989/147‑1 (MQ=255)                                                                     gCAGGGGGGCGATCATCTTTTTTTCACC‑AAAAAAAATCCCTGCCTGCGGCCCGCCTAACAACTGGTCGCCGGAGAAACC‑CCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                            <  1:1236156/150‑1 (MQ=255)                                                                     gCAGGGGGGCATTCTTCTTTTTTTTCCC‑AACAATAATTCCTGCCTGCGGCCCGCTTAACAACTTGTCGCCGGAAAAACT‑CCCCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAt                                                                            <  1:536843/150‑1 (MQ=255)                                                                         gcggccgACAATCTTTTTTTTCCC‑AACAAAATTCCCTCCCTGCGCCCCCCCTAAAAACTTGTCCCCGGAAAACCC‑CCCCAGCCTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCg                                                                          <  2:626858/148‑1 (MQ=255)                                                                          cgggcgATCATCTTTTTTTTCCC‑AACAAAAATCCCTGCCTGCGCCCCGCCTAACAACTTGTCGCCGGAAAAACC‑CCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGAcc                                                                       <  1:768881/150‑1 (MQ=255)                                                                                      cgttttttCCC‑AAAAAAAATCCCTGCCTGCGCCCCCCCAAAAAACTTCCCGCCGGAAAAACC‑CACCCGGCGGGCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGc                                                          <  1:340971/149‑1 (MQ=255)                                                                                             tccc‑AAAAATAATCCCTACCGGCGGCCCCCAAAACATCTTGTCGCAGGGGAACCC‑CCCCGGCAGGACACCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATc                                                   <  2:1244083/149‑1 (MQ=255)                                                                                             tccc‑AAAAAAAATCCCTCCCTGCGGCCCGCCAAACAATTTGCCGCCGGAAAAACC‑CCCCAGCCTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATc                                                   <  1:117487/149‑1 (MQ=255)                                                                                             tccc‑AAAAAAAACTCCCGCCCGCGGCCCCCCTAACACTTTGCGGCGGGAAAAACC‑CCCCAGCCGGGCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATc                                                   <  1:165796/149‑1 (MQ=255)                                                                                             aaCC‑AAAAAAAATTCCTGCCTGCGCCCCGCATAACAATTGGCCGCGGAAAAAACT‑CCCCAGTCTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATc                                                   <  1:1015913/150‑1 (MQ=255)                                                                                                  cacaaTAATCCCGCCCAGCGCGCCGCATAA‑AATCTGTCCGCGGTGAAACCCCACCAGACGCACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGta                                               <  2:282686/150‑2 (MQ=255)                                                                                                    caaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCGATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCa                                              >  2:884114/1‑150 (MQ=255)                                                                                                    caaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCa                                              >  2:981890/1‑150 (MQ=255)                                                                                                    caaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCa                                              >  1:331395/1‑150 (MQ=255)                                                                                                            cctgccggCGCCCCGCCTAAAAACTTGTCGCCGGAAAACCC‑CACCAAACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGaa                                     <  2:469502/151‑1 (MQ=255)                                                                                                            cctcccggCGGCCCGCCAACCAATTTGTCGCGGGAAAAACC‑CCCCAGACTGCCGCCCGCCGCAATCAACTTCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGaa                                     <  1:1043792/151‑1 (MQ=255)                                                                                                             cggccggcGCCCCCCCAACCAACTGGTCGCGGGAAAACCT‑CACCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAAc                                    <  2:436837/149‑1 (MQ=255)                                                                                                             ccccctgcGGCCCCCCTAAAACCTTGCCGCCGGAAAACCC‑CCCCAACCAGACTCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAAc                                    <  2:368725/148‑1 (MQ=255)                                                                                                                      gCCCGCCGAACAACTGGTCGCCGGAGAAACT‑CCCCAGCCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGGGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGct                            <  2:168835/150‑1 (MQ=255)                                                                                                                       cccccccaaaaccTTGCCCCCGAAGAACCC‑CCCCAGACGGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGctc                           <  1:1101030/138‑1 (MQ=255)                                                                                                                       cccGCCTAACAACTTGCCGCCGGAAAACCT‑CCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGctc                           <  1:1092707/150‑1 (MQ=255)                                                                                                                          gCCCAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCGTACTGGCTAAGATCGACCCGCGAGCCACTGCCTAAATCAGTCACTACCGGAACCTCCAGCTCTTTg                       >  2:222695/1‑151 (MQ=255)                                                                                                                           ccaaacactttGTCCCGGGAAAACCC‑CCCCAGCCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGc                      <  1:274525/142‑1 (MQ=255)                                                                                                                           ccTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTg                       >  2:593281/1‑150 (MQ=255)                                                                                                                                   cTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgcgc               >  1:801690/1‑150 (MQ=255)                                                                                                                                       ccGCCGGAAAACCC‑CCCCAGCCTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGccacc           <  1:832528/149‑1 (MQ=255)                                                                                                                                                aaaccccACCATACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCAGCt  <  2:958628/145‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                CAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACT‑CACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCT  >  NC_000913/3760236‑3760527

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