Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,717,771	0	A	G	24.0%	30.3 / 11.2	25	G154G (GGA→GGG)	ghrB	glyoxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (8/11);  new base G (3/3);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGG  >  NC_000913/3717555‑3718016                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      cGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGAACGGCCTGGCGCTGGCCCAACGTgcgccct                                                                                                                                                                                                                     >  2:194999/1‑248 (MQ=255)          tACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTc                                                                                                                                                                                                             <  1:156353/250‑1 (MQ=255)                         tCTGCAGATGCACACGCCAACCGGATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGGTGTCTACCGCTCGTCGGGGTGTGGATGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTGACCATAAAACGCTGGGCATTGTCGGGAGGGGACGGATCGGCATGGCGCTGGCCCAACTTGCGCACGGTGGCTTCAGCAGGCACATCctc                                                                                                                                                                                              <  2:179663/250‑1 (MQ=255)                            gctgATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGGCGGATCGGCATGGCTCGGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCtata                                                                                                                                                                                          <  1:22969/251‑1 (MQ=255)                              tgATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGGCGGATCGGCATGGCGCTGGCACCACGTGCGCACTTTGGCTTTCACATGCCCCTCCTCTATAAc                                                                                                                                                                                        >  2:60500/1‑251 (MQ=255)                                aTGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAAcgc                                                                                                                                                                                      <  2:232909/251‑1 (MQ=255)                                  gCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAAcgcgc                                                                                                                                                                                    <  1:3671/251‑1 (MQ=255)                                          cAACCGTATTAAAAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATGGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATACCGCGCGCCGccac                                                                                                                                                                            <  2:13680/251‑1 (MQ=255)                                           aaCCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGATACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGGCGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGccacc                                                                                                                                                                           <  2:183728/251‑1 (MQ=255)                                                     aCAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCATTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACCCGCGCCGCCACCCTAAAGAAGc                                                                                                                                                                 >  2:189091/1‑251 (MQ=255)                                                       aGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGGCGGATCGGCATGGCGCTGGCACAACGGTCGCACTTTGGCTTCAACATGCCCATCCTATATAACGCGCGCCGCCACCATAAAGAAGc                                                                                                                                                                 >  2:13518/1‑249 (MQ=255)                                                       aGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGc                                                                                                                                                                 >  1:211526/1‑249 (MQ=255)                                                                 gCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGCGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGGGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCCTAAAGAAGCAGAAGAACGcct                                                                                                                                                     >  2:146570/1‑249 (MQ=255)                                                                       aCGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGCGGCAGCAGCACGGGTAAAAGCAGGCGAATGCACCGCGAGCATAGGCCCGGACTGCTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGACCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGACCGCTTCACCGc                                                                                                                                               >  1:30277/1‑251 (MQ=255)                                                                        cGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGcc                                                                                                                                              >  1:46921/1‑251 (MQ=255)                                                                             aTGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTa                                                                                                                                         <  1:240182/251‑1 (MQ=255)                                                                               ggCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGCCCGGGTAAAAGCAGGCGAATGGCCCGCGAGCATAGGACCGGACTGGTACGGCACTGACGTTCAACATAAAACACTGGGCATTGTCGGGCTGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCAATTTGGCTTCAACATGCCCATCCTCTATAACGAGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCACGCTACt                                                                                                                                       >  1:152951/1‑251 (MQ=255)                                                                               ggCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGGAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTGGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAATGCCCGCTACt                                                                                                                                       <  2:137303/251‑1 (MQ=255)                                                                                               aCCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATActct                                                                                                                       <  1:143680/251‑1 (MQ=255)                                                                                                     cgtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGGCGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACa                                                                                                                 <  1:207605/251‑1 (MQ=255)                                                                                                                           aaCGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTg                                                                                           >  1:268029/1‑251 (MQ=255)                                                                                                                                                                                   gACGTTCACCATAAACCACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATACCGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCaa                                   <  2:211526/251‑1 (MQ=255)                                                                                                                                                                                                 aaCACTGGGCATTGTCGGGGTGGGACGGATCGGCAGGGCGCTGGCACAATGTGCGCATTTTGGCTTCAACAAGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGATATTGCGATTTGGATACTCTGTTACAAGAGTCAGTTTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTGTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCg                     <  2:156202/251‑1 (MQ=255)                                                                                                                                                                                                                  ggaggggACGGATCGGCATGGCGCTGGCACACCGTGCGCACTTGGGCTTCAACATGCCCATCCTCTAAAACGCGCGCCGCCACAATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGcc    <  2:199712/247‑1 (MQ=255)                                                                                                                                                                                                                     tGGGGCGGATCGGCATGGCGATGGCACAACGTGCGCACTTTGGCTTAAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCgg  >  1:87457/1‑250 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      CGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAAAATGAAATCCTCCGCCATTTTCATTAATGCCGG  >  NC_000913/3717555‑3718016

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