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

GCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689395‑689661                                                                                                                        |                                                                                                                                                    gCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgt                                                                                                                      >  2:612438/1‑151 (MQ=255)           cccACAGCGGCCAGTTGCCGTTCCGGGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgct                                                                                                            <  2:168817/151‑1 (MQ=255)                   ggctggTTGCGGGTACGTGCGTATGGTGTGAGTCCGGTGGTGGGCGTCACTTTAGTGGTTAACCCCTCGCGGTTAAACGGCGGGATCATCGCTTGAGTCTCCCCCTGCGGACCAATCACCGCCGTAATGCCGTTGTTGGtgctg                                                                                                           <  1:111009/139‑1 (MQ=255)                           gCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcg                                                                                             >  2:685769/1‑150 (MQ=255)                            ccGGTACGTGCGTAGGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgc                                                                                            <  2:436001/150‑1 (MQ=255)                               gTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGcc                                                                                         >  2:94465/1‑150 (MQ=255)                                                       gtggtCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCa                                                                <  2:195403/151‑1 (MQ=255)                                                                      tgaGTGGTTACCCCCTCGCGGGTGAACTGCGTGTTCATCGCCTAATTCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                        <  2:534853/142‑1 (MQ=255)                                                                       tagggtttAACCCCTCGCGGGTGAACTGCGGGATCACCGCCTGAATCCCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGc                                                  <  1:94428/143‑1 (MQ=255)                                                                               aaCCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCaa                                        <  2:287708/151‑1 (MQ=255)                                                                                   catCGCGGGGGAACGGCGGGACCATCGCCTGAACCCCCCCCTGGGGGCAATTCTCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGa                                    <  2:62963/149‑1 (MQ=255)                                                                                     tCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGGGTTGCCATGGACCAATAGAtt                                  >  1:649707/1‑151 (MQ=255)                                                                                      cgcgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 >  2:648317/1‑151 (MQ=255)                                                                                      cgcgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 >  2:807558/1‑151 (MQ=255)                                                                                      cgcgGGTGAACTGCGGGATCATCCCCTTAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAtt                                  <  2:219207/150‑1 (MQ=255)                                                                                         gggaaaaCTGCGGGATCTTCCCCTGAATCTCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 <  1:637482/143‑1 (MQ=255)                                                                                              aaCGGCGGGATCCTCGCCTAAATCGCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAcc                         <  1:555573/151‑1 (MQ=255)                                                                                                         catcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTg              >  1:116692/1‑151 (MQ=255)                                                                                                                        aCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  2:28308/1‑148 (MQ=255)                                                                                                                        |                                                                                                                                                    GCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689395‑689661

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: ‑