Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	689,514	0	A	C	42.6%	89.9 / 14.5	55	G456G (GGT→GGG)	lnt	apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/4);  new base C (0/23);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 3.23e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTC  >  NC_000913/689367‑689658                                                                                                                                                     |                                                                                                                                                 catcaACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCAccc                                                                                                                                               <  2:314109/151‑1 (MQ=255)    caACACGTCAGCAAAGCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGc                                                                                                                                            <  1:506265/151‑1 (MQ=255)        acaGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGc                                                                                                                                        >  1:211695/1‑151 (MQ=255)                aaaaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAcc                                                                                                                                <  2:827339/151‑1 (MQ=255)                        acaagcccttcaccccccACAGCGGCCAGTTCCGGGTACGGGCTTATGGTGT‑GTGTCCGGTGGCCGGCGCCCCTTTAGTGGTAACCCCCCGGCGGGTGACCGGGGGGATCTCCGCCTGAACCTCCCCCTGCGGGCCAATCACCGCCGTAAt                                                                                                                        <  1:399504/136‑1 (MQ=255)                          aaTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGc                                                                                                                      >  2:272590/1‑151 (MQ=255)                              ccGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgtt                                                                                                                  >  2:739148/1‑151 (MQ=255)                                gTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGGGGGCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttg                                                                                                                 >  1:861499/1‑150 (MQ=255)                                    gCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGGGGGCGGCGTCACGTTAGAGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGt                                                                                                            >  1:976663/1‑151 (MQ=255)                                      aCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgc                                                                                                          >  2:412910/1‑151 (MQ=255)                                      aCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGGGGGCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgc                                                                                                          >  1:764834/1‑151 (MQ=255)                                       cccACAGCGGCGTGTTTCCGGTCCGTGGTGATGGTGG‑GAGTCCGGTGGTCGTGGCCACGTTAGTGGTTAACCCCGCGGGGTTGAACTGGGGGAGCACCGGTTGAATCTCACCCGGGGGCCAAACCACCGCCGTAATGCCGTTGTTGGtg                                                                                                           <  1:798836/149‑1 (MQ=255)                                       cccACAGCGGCCATTGGCGGGTACGTGCGTATGGTGT‑GGGTCCGGTGGTCGGCGTCAGGTTAGTGTTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTAAATCTCCCCCGGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGtgct                                                                                                         <  2:883091/151‑1 (MQ=255)                                         cacaGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCCCGCGGGTGAACTGCGGGATCATCGCCTGATTCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgctgc                                                                                                       <  1:76985/151‑1 (MQ=255)                                          acaGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcg                                                                                                      >  2:27395/1‑151 (MQ=255)                                            aGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcgc                                                                                                     >  1:184168/1‑150 (MQ=255)                                                gCCAGTTGCCGGTACGTGCGTATGGTGT‑GAATCTGGGGGTCCGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCaa                                                                                                   >  1:918009/1‑148 (MQ=255)                                                     ttGCGGGTACGTGCGTATGGTGG‑GGGTCCGGTGGTCGGTGTAACGTTAGTGGTAAACCCCTCGCGGGTGACCTGCGGGATCATCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCa                                                                                                    <  2:491598/142‑1 (MQ=255)                                                         cGGTACGTGCGTATGGTGTAGATCTCGGTGGTCGCCGTATCATTA‑AAAAAAACACATCGCGGGTGAAATGCGGGATCATCGACTGAATCCCACCCTGAGGGACAATCACCGCCGTAATGCCGTTGTTGGTGCTCCGCAACCGTGGgcgc                                                                                         >  1:586699/1‑149 (MQ=255)                                                          ggTACGGGCGTATGGTGT‑GAGTCCGGTGGGCGGCGCCCCGTTAGTGGTTACCCCCTCGGGGGTGACCTGGGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                               <  1:963543/142‑1 (MQ=255)                                                                 gCGTATGGTGT‑GAGTCGGGTGGTCGGCGTCACGTTAGTGGTTACCACCTCGGGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCGGGGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCa                                                                               <  2:585706/151‑1 (MQ=255)                                                                   gTATGGTGT‑GAGTCCGGTGGTCGGGGTCACGTTAGTGTTTAACCCCTCGCGGGTGAGCGGCGGGATCATCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAg                                                                              <  1:639490/150‑1 (MQ=255)                                                                    tATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgc                                                                             <  2:908535/150‑1 (MQ=255)                                                                       ggTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcac                                                                         <  1:336591/151‑1 (MQ=255)                                                                         tgt‑gAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgc                                                                       >  1:823971/1‑151 (MQ=255)                                                                          gt‑gAGTCCGGTGGTCGGCGCCACGTTAGTGGTTAACCCCCCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgc                                                                       <  1:792834/150‑1 (MQ=255)                                                                                      ggtGGGCGTCACGTTAGTGGTTACCCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATc                                                           <  2:1056818/151‑1 (MQ=255)                                                                                         cGGCGTCACGTTAGTGGTTAACCCCTCGGGGGTGAACTGCGGGATCATCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTg                                                         <  1:1152811/150‑1 (MQ=255)                                                                                          ggCGTCACGTTAGTGGTAAACCCCCCGGGGGTGAACGGCGGGATCCTCGCCTAAACCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCt                                                          <  1:654473/148‑1 (MQ=255)                                                                                          gctgTCACTTTAGTGGTTAACCCCTCGGGGGTGACCGGCGGGATCACCGCCTGAACCCCCCCCGGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacg                                                                        <  1:1082211/131‑1 (MQ=255)                                                                                               cccGTTGGTGGTTAACCCCTCGGGGGTAAACTGGGGGATCACCGCCTGAATCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgt                                                    <  1:69419/147‑1 (MQ=255)                                                                                                aCGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                 >  1:775304/1‑151 (MQ=255)                                                                                                    tAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGcc                                              >  2:935869/1‑150 (MQ=255)                                                                                                      gTGGTAACCACCTCGGGGGTGAACTGGGGGATCACCGCCTGAATCTCCCCCGGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCa                                             <  1:630006/149‑1 (MQ=255)                                                                                                       tgctcaccacCTCGGGGGTAACCTGTGGGATCAGCGCCTGAATTTCTCCCTGCGGTCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATg                                           <  2:617864/143‑1 (MQ=255)                                                                                                        ggTTAACACCCCGGGGGTGACCGGCGGGATCACCGCCTGAATCTCCCCCGGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGa                                         <  1:1114483/151‑1 (MQ=255)                                                                                                                 cTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                    >  2:1139861/1‑147 (MQ=255)                                                                                                                    gcgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                              >  2:407281/1‑150 (MQ=255)                                                                                                                     cgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAc                            >  1:991571/1‑151 (MQ=255)                                                                                                                        gTGAACTGGGGGATCATCGCCTGAATCTCCCCTGGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCaa                         <  2:294338/151‑1 (MQ=255)                                                                                                                           aaCTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAcc                      >  2:183408/1‑151 (MQ=255)                                                                                                                             cggggggATCATCCCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCa                     <  1:445704/146‑1 (MQ=255)                                                                                                                              tGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg                   >  2:998853/1‑151 (MQ=255)                                                                                                                               gggggATAATCGCCTAATTCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg                   <  1:866080/148‑1 (MQ=255)                                                                                                                               gCGGGATCATCCCCTAAATCCCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGc                  <  2:900324/151‑1 (MQ=255)                                                                                                                                 gggATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                >  2:411459/1‑151 (MQ=255)                                                                                                                                     tcatcGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGt             <  2:976697/150‑1 (MQ=255)                                                                                                                                       atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg          >  1:433129/1‑151 (MQ=255)                                                                                                                                          gCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg          >  2:600605/1‑148 (MQ=255)                                                                                                                                          gCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg          >  1:362196/1‑148 (MQ=255)                                                                                                                                           ccTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg          >  2:1050184/1‑147 (MQ=255)                                                                                                                                           ccTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg          >  1:734699/1‑147 (MQ=255)                                                                                                                                           ccTGAATCCCCCCCGGCGGGCCAGTCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg          <  2:67714/147‑1 (MQ=255)                                                                                                                                             tGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt      >  1:53089/1‑149 (MQ=255)                                                                                                                                                 tctcACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTc  >  2:142936/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTC  >  NC_000913/689367‑689658

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