Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	689,514	0	A	C	51.3%	44.9 / 56.7	37	G456G (GGT→GGG)	lnt	apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (13/5);  new base C (0/19);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 2.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-06

CATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAG  >  NC_000913/689367‑689652                                                                                                                                                     |                                                                                                                                           catcaACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCAc                                                                                                                                           <  1:143889/149‑1 (MQ=255)     aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCg                                                                                                                                     >  1:535102/1‑151 (MQ=255)         cagcagCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGcc                                                                                                                                 >  1:412696/1‑151 (MQ=255)            cagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAAt                                                                                                                              <  1:426793/151‑1 (MQ=255)                  aaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAccgc                                                                                                                        <  2:596049/151‑1 (MQ=255)                      aaaCAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGt                                                                                                                     >  1:521098/1‑150 (MQ=255)                          aaTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGc                                                                                                                >  1:310610/1‑151 (MQ=255)                                       cccACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgct                                                                                                   >  1:516160/1‑151 (MQ=255)                                               ggcaggttccgggTAGGTGCTAAGGGTGGGAGTCCGGGGGTGGGCGTCCCGTTATTGGTTAA‑ACCCCCGGGGGTGAACGGGGGGATCATCGCCGGAATCTCCCCGGGCGGGCCAAGCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACa                                                                                           <  2:751426/140‑1 (MQ=255)                                                     ttGCCGGTACGGGCGTATGTTGTGAGCCCGGGGGTCGGCGTCACGTTAGTGGTAAA‑CACCGCGCGGGTGAACTGGGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTg                                                                                        <  2:738072/148‑1 (MQ=255)                                                      ggcggttccGTGCGTTGGGTGGGGGTCCGGTGGCCGGGGCCACGTTAGTGGTTAC‑CCCCCCGCGGGTAAACTGCGGGATCACCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcg                                                                                    <  2:380540/143‑1 (MQ=255)                                                       gcggttACGTGCGTAGGGGGGAGGTCGGGTGGTGGGGTCCACGTTAGTGTTTAC‑CCCCTCGCGGGTGAACTGCGGGTTCACCCCTTGAATCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgc                                                                                     <  2:323244/144‑1 (MQ=255)                                                       gCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgc                                                                                   >  1:766886/1‑151 (MQ=255)                                                                       ggTGGGAGTCCGGTGGTGGGCGTCACGTTAGTGGTTAA‑CCCCTCGCGGGTGAACGGCGGGTTCATCGCCTGAACCTCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgcgc                                                                     <  1:467856/149‑1 (MQ=255)                                                                               tCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcgc                                                           >  2:663530/1‑151 (MQ=255)                                                                                   gtggtCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGcc                                                        >  2:611209/1‑150 (MQ=255)                                                                                     ggtGGGCGTCGCTTTAG‑GGTTAACCCCCCCGCGGTTAAACTGGGGGATCCTCGCCTGAACCCCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATc                                                     <  2:588851/151‑1 (MQ=255)                                                                                         ggCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                  >  1:432063/1‑150 (MQ=255)                                                                                         ggCGTCACGTTAGTGGTTAA‑CACCCCGCGGGTGAACTGCGGGATCATCGCCGGAACCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTg                                                   <  1:777299/149‑1 (MQ=255)                                                                                              cACTTTAGTGGTTAA‑CCCCTCGCGGGTGAACTGGGGGATCATCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtgt                                            <  1:295051/151‑1 (MQ=255)                                                                                              cACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCCTCTGGAAggg                                             >  2:79309/1‑148 (MQ=255)                                                                                               aCGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtgt                                            <  1:29421/150‑1 (MQ=255)                                                                                               aCGTTAGTGATTAA‑CACCTCGCGGGTGAACGGGGGTACAATCGCTTGAATCTCCCCGCGCGGGCCAATTGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                             <  1:649021/149‑1 (MQ=255)                                                                                                        gTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                  >  1:675628/1‑151 (MQ=255)                                                                                                                ccggaggggTGAACGGGGGGATCCTCCCCTGACCCTCCCCGGCCGGGCCAAACACCGCCGTAATGCCGTTATTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                              <  2:635687/142‑1 (MQ=255)                                                                                                                ccTCGCGGGTAAACTGCGGGATCACCCCCTGAACCTCCCCCTGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAg                            <  1:176419/150‑1 (MQ=255)                                                                                                                ccTCGCGGGCGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCACGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAg                            >  2:412192/1‑150 (MQ=255)                                                                                                                    gcgGGAAAACGGCGGGCTCAGCGCCTGAACCTCCCCCTGCGGGCCCATCACCGCCGTAATGCCGTTGGTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGAACGCATTCGCGCCATCTGGAAGTGTTGCCATGGGCCAATAGAtt                         <  2:460293/149‑1 (MQ=255)                                                                                                                           agcggggggacaagcgcgTGAATCCCCCCCGGGGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAc                 <  2:570528/136‑1 (MQ=255)                                                                                                                             cTGGGGGATCATCGCCGGAATCTCCCCCGGCGGGCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAcc                <  2:445994/149‑1 (MQ=255)                                                                                                                               gggggTTCATCGCTTGAATCTCCCCCTGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCa               <  2:727529/146‑1 (MQ=255)                                                                                                                                 gggATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCACCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCAc              <  1:512864/147‑1 (MQ=255)                                                                                                                                       atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg    >  1:179005/1‑151 (MQ=255)                                                                                                                                       agccccTGAACCCCCCCCGGTGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTg     <  1:614987/146‑1 (MQ=255)                                                                                                                                        tcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg    <  1:111966/150‑1 (MQ=255)                                                                                                                                         cGCCTGAATCTCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGgag  <  2:520806/151‑1 (MQ=255)                                                                                                                                          gCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGt       <  2:766519/145‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                           CATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAA‑CACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAG  >  NC_000913/689367‑689652

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