Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	689,514	0	A	C	52.4%	17.9 / 25.1	21	G456G (GGT→GGG)	lnt	apolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (8/2);  new base C (0/11);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 2.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-04

AGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGA  >  NC_000913/689376‑689653                                                                                                                                           |                                                                                                                                            agcagcAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGcc                                                                                                                                  >  2:1146601/1‑150 (MQ=255) agcagcAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCa                                                                                                                                 <  1:2520649/151‑1 (MQ=255)                              cccACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgct                                                                                                    <  1:372800/151‑1 (MQ=255)                                      ggCCAGTTGCCGTTAGGTGCGTATGGTGTGAGTCCGGGGGTCGGCGTCACGTTAGTGGTTAACCCCCCGCGGGTGAACTGGGGGATCACCGCCTGAACCCCCCCCGGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCaa                                                                                              <  1:3471573/149‑1 (MQ=255)                                         gaGTTGCCGGTACGTGCGTATGGGGAAAGTGTGGCGCCCTGCGGTAAGTTCGTGGGTAGCCCCCCGCGTGTAAAAAAAAAAATCATCGCCTGAATCTCACCCTGCGCGCCACTCACCGCCGTCATGCCGTTGTTGGTGCTGCGCACCAgcg                                                                                         >  1:1167464/2‑149 (MQ=255)                                                 ggTACGTGCGTATGGTGTGGGTCCGGTGGTGGGCGTCACTTTATTGTTTACCCCCTCGGGGGTGAACGGCGGGTTCATCGCCTGATTCCCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcg                                                                                   <  1:2148985/149‑1 (MQ=255)                                                   tACGTGCGTAGGGTGTGAGCCCGGTGGGCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGGGGGATCTTCGCCTAAACCCCCCCCTGCGGCCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCa                                                                                <  1:1811039/150‑1 (MQ=255)                                                          gTATGGGGTGTGCCGGTTGGTCGGCGTCACGTTAGTGTTTAACCCCTCGCGGGTGACCTGGGGGTTCATCCCCTGAATCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGc                                                                              <  2:3827018/145‑1 (MQ=255)                                                                                 gCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                  >  1:2453288/1‑150 (MQ=255)                                                                                    tCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                              >  1:3562407/1‑151 (MQ=255)                                                                                    tCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgt                                               >  2:2917004/1‑150 (MQ=255)                                                                                     cACGTTAGTGTTTACCCCCTCGCGGGTGAACTGCGGGATCACCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGaa                                                 <  1:3256208/147‑1 (MQ=255)                                                                                         ttAGTGTTGAACGTCTCGGGGGTGATCTGCGGAACACCCCCCTGACCCCCACCCTGGGGGCAAAACTCCGCCTTAATGCCGTTGTTGGTCCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                            <  1:2113702/148‑1 (MQ=255)                                                                                                          gcgGGTAACCTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTa                        <  2:2276830/151‑1 (MQ=255)                                                                                                                  aCTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAc                  <  2:434510/149‑1 (MQ=255)                                                                                                                    tGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg              >  2:3929863/1‑151 (MQ=255)                                                                                                                        ggATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATc          <  1:1095336/151‑1 (MQ=255)                                                                                                                             atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg     >  1:3141713/1‑151 (MQ=255)                                                                                                                             atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg     >  1:694764/1‑151 (MQ=255)                                                                                                                                 ccTGAATCCCCCCCTGCGGTCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTCCCAAACCACGCATCGTTgg     <  1:3224962/147‑1 (MQ=255)                                                                                                                                    gAACCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGgaga  <  1:874674/147‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            AGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGA  >  NC_000913/689376‑689653

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