Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	689,514	0	A	C	45.0%	53.3 / 10.6	40	G456G (GGT→GGG)	lnt	apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (21/1);  new base C (0/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.68e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGT  >  NC_000913/689370‑689657                                                                                                                                                  |                                                                                                                                                caACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAGATCTCGGTGGTCGCCGTATCATTA‑AAAAAAACAACTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGc                                                                                                                                           >  1:1043573/1‑151 (MQ=255)           gcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATc                                                                                                                                  >  2:988623/1‑150 (MQ=255)             aaaaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAcc                                                                                                                               >  2:1035286/1‑151 (MQ=255)               aaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAccgc                                                                                                                             >  2:1038048/1‑151 (MQ=255)                              tCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttgtt                                                                                                              >  2:237267/1‑151 (MQ=255)                                      cacaGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑AGGTCCCGTGGTCCGCCTCACCTTTACGGTTAACACCTCGCGGGTTAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTACTGCCGTTGTTGGtgctgc                                                                                                      >  1:535294/1‑151 (MQ=255)                                            ggCCAGTTGCGGTTACGTGCGTATGGTGT‑GGGTCCGGTGGTGGGCTTCCCGTTAGTGGTTACCCCCCCGCGGGTAAACTGGGGGATCTTCGCCTAAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCaa                                                                                                  <  1:1121331/149‑1 (MQ=255)                                            ggCCAGTTGCCGGTACGTGCGTATGGTGT‑AGATCCCGTGGTGCGCCTCACCATAAAGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTAACCCTGCGGGGCACTCACCGAGGCAATGCCGGTGTTGGTGCTGAGCAAc                                                                                                 >  1:913025/1‑150 (MQ=255)                                                aGTGGCCGGTACGGGCGTAGGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTAAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                              <  2:136741/149‑1 (MQ=255)                                                       gggaggggCGTTTGGTGT‑GAGCCGGGTGGTCGGCGTCGCGTTGGTGGTTAACCCCTCGCGGGGGAACTGCGGGATCACCGCCTGAATCCCCCCCGGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcgc                                                                                      <  1:768076/143‑1 (MQ=255)                                                       ggTACTTGCGTATGGTGG‑GGGTCCGGTTGTCGGCTACCCGTTAGTGGTTACCCCCGCGCGGGGGAACGGGGGGATCACCGCCTGATCCTCCCCCTGCGGGCCAACCTCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcg                                                                                       <  1:246476/149‑1 (MQ=255)                                                        gTAGGTGTGTAGGGTGG‑GGGTCCGTTGGCGGGCGGCGGGTTGGTGTTAACCCCCTCGCGGGTGAACTGGGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCa                                                                                    <  1:618232/151‑1 (MQ=255)                                                            gTGCGTATGGTGT‑AAGTCTGGGGGTCCGCCGCACGTTTGAGGATAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTc                                                                                >  1:847856/1‑151 (MQ=255)                                                             tGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                               >  2:62028/1‑151 (MQ=255)                                                              gCGTATGGTGG‑GGGTCCGGTGGTCGGCGCCACGTTAGTGGTTACCTCGGCGCGGGTGAACTGCGGGAACAGCGCCTGAACCCCCCCTTGCGGGCCAAGCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCa                                                                              <  1:1131425/151‑1 (MQ=255)                                                                gTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgc                                                                            >  2:1130805/1‑151 (MQ=255)                                                                    ggTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcac                                                                        <  2:344643/151‑1 (MQ=255)                                                                       gt‑gtgtCCGGTGGCGGGCGTCACGTTAGTGGTTACCCCCTCGCGGGGGAACTGCGGGATCACCGCTTGACTCTCCCCCTGCGGGCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgca                                                                     <  1:546153/151‑1 (MQ=255)                                                                         tgcgTCGGGTGGTGGGCGTCACGTTAGTGGTGACCCCCTCGCGGGTGAACTGCGGGATCATCCCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgca                                                                     <  1:845332/147‑1 (MQ=255)                                                                            gTCGGGGGGTCGGCGCCACGTTAGTGTTTAACGCCTCGGGGGTAACTTGGGGGATCATCGCCGGAACCTCCCCCGGGGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                 <  1:106515/151‑1 (MQ=255)                                                                            gTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                 >  2:146707/1‑151 (MQ=255)                                                                             tCCGGTGGTGGGCGTCACGTTAGTGTTTAACCCCTCGGGGGTGAACTGGGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                 <  2:91609/150‑1 (MQ=255)                                                                             tCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcgc                                                                >  2:457971/1‑151 (MQ=255)                                                                                ggtggtCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGcc                                                             >  1:339181/1‑151 (MQ=255)                                                                                      cGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                       >  1:688783/1‑151 (MQ=255)                                                                                       gggtcccctttAGTGGTTAACCCCTCGCGGGTGACCGGGGGGATCATCCCCAGAACCCCACCCTGCGGGCCAAGCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTg                                                        <  1:449993/140‑1 (MQ=255)                                                                                       ggCGTCACGTTAGTGGTTAACCCCCCGCGGGTGAACTGCGGGATCATCGCCTGAATCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                       <  2:331115/150‑1 (MQ=255)                                                                                               gTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTg                                               <  2:924927/150‑1 (MQ=255)                                                                                               gTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTg                                               >  1:237465/1‑150 (MQ=255)                                                                                                                  cgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAc                           >  2:205348/1‑151 (MQ=255)                                                                                                                   gggTGAACTGCGGGATCATCGCCTGAAACTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                             <  2:578209/148‑1 (MQ=255)                                                                                                                   gggTGAACTGCGGGATCATCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGagggacc                          <  2:904045/151‑7 (MQ=255)                                                                                                                    ggagAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAcc                          >  1:114911/4‑150 (MQ=255)                                                                                                                    ggagAACTGCGGGATCATCGCCCGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGCTGTTGGTGCTGCGCAACAGAGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAc                           >  1:446665/4‑149 (MQ=255)                                                                                                                            gggggATCCGCGCCTGATTCCCCCCCTGGGGTCAATTCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCAc                   <  2:705000/147‑1 (MQ=255)                                                                                                                              ggtatcatcGCCTGAATCTCCCCCTGCGGGCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg                  <  2:163663/145‑1 (MQ=255)                                                                                                                                    atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg         >  2:514072/1‑151 (MQ=255)                                                                                                                                        ccTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg         >  2:680048/1‑147 (MQ=255)                                                                                                                                            attctcCCCCGGCGGGCAAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGtt           <  1:1034792/139‑1 (MQ=255)                                                                                                                                               ctcACCCTGCGCGCCAATCACAGCCGTAATGCCGATGTTGGTGCTGCGCAACAGTGGGCGCGCAAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGt  >  1:1019169/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                CAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGT  >  NC_000913/689370‑689657

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