Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,079,468	0	A	G	66.7%	‑1.3 / 13.1	18	D55G (GAT→GGT)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (0/6);  new base G (12/0);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 5.39e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-04

TGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCTG‑TTAAT‑GGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079324‑1079616                                                                                                                                                 |                                                                                                                                                        tGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAg                                                                                                                                                   <  2:337820/151‑1 (MQ=255)                                   gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGTGGGTATGAGCGGGTGGGCGT‑TTAAT‑GGGGTGGCggggggtggt                                                                                                               >  2:406385/1‑145 (MQ=255)                                    atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTAGGGTATGGGGGGCT‑GGCTG‑GTTAT‑GGGGGTGGCGGgggcgggtt                                                                                                             >  2:386627/1‑146 (MQ=255)                                    atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGCGGCG‑GGGTG‑TTTAT‑GGGGGTGCCggggggcgttt                                                                                                             >  2:363537/1‑143 (MQ=255)                                    atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTAGGGGGC‑GGGGG‑GTTAT‑GGGGGTGGCggggggtggtt                                                                                                             >  2:8335/1‑143 (MQ=255)                                      attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATGAGCGGGT‑GGGTG‑TTTAT‑GGGGGTGCCGGGCGCTGgtttt                                                                                                           >  2:146067/1‑151 (MQ=255)                                              ggTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGAATTGGAGGGG‑GGCGG‑GTGAT‑GGGGTTTCCGGGCGCCGGTTTTTTTTTCgg                                                                                                   >  1:204036/1‑151 (MQ=255)                                                           gCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGGT‑GGCTGTTTATT‑GGGGTTG‑CGGGGGCCGGTTTTTTTTTCGGGGTTTTCGAAAgg                                                                                      >  2:57216/1‑150 (MQ=255)                                                                     aaaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCTG‑TTAAT‑GGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGcc                                                                             <  2:293457/150‑1 (MQ=255)                                                                      aaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCTG‑TTAAT‑GGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAt                                                                           <  1:57216/151‑1 (MQ=255)                                                                           tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCTG‑TTAAT‑GGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGcc                                                                      <  1:257837/151‑1 (MQ=255)                                                                                  acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATGAGCGGGT‑GGCTG‑TTTAT‑GGGGGTGGCGGGGGCTGTTTTTTTTTTCGGGGTTTTCGAAAGCTGGGTCGCCATTGGCCTGACAtt                                                               >  2:241807/1‑151 (MQ=255)                                                                                       atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGGGGT‑GGGTG‑TTTAT‑GGGGGTGCCGGGGGCTGTTTTTTTTTTCGGGGTTTTCGAAAGCTGGGACGCCCTTGGCCTTGCATTTgggg                                                          >  2:303367/1‑149 (MQ=255)                                                                                                 gTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGGGGCC‑GGCGG‑GTTAT‑GGGGTTGTCGGGCGCCGGTTTTTTTTCTGGGGGTTCCGCAAGCCGGATTGGCCCTGGCCTGACATTATGCGCGGGGGttta                                                >  1:148544/1‑149 (MQ=255)                                                                                                 gTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGCC‑GGGTG‑TTTATGGGGGTTGCCGGG‑GCTGGTTTTTTTTCCCGGGTTTCCGAAAGCTGGGTCGCCATATGGCTGACATTTAGCGCGGGGGTTaa                                                >  2:188952/1‑151 (MQ=255)                                                                                                       gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATTAGCGGCT‑GGGTGTTTAAT‑GGGG‑TGGCGGGGGCTGTTTTTTTTTTCGGGGTTTTCGAAAACTGGGTTGCCCTTGGGCTGACATTTGGGGCGGGGGTTAAATggga                                          >  2:152925/1‑149 (MQ=255)                                                                                                             ggCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCTG‑TTAAT‑GGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtgg                                     <  1:365178/150‑1 (MQ=255)                                                                                                                                                gATATGAGCGGCT‑GGCTG‑TTAAT‑GGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaac  <  2:142647/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                        TGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCTG‑TTAAT‑GGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079324‑1079616

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