Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,969,794	0	A	C	52.6%	22.2 / 31.1	19	G397G (GGT→GGG)	tap	methyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (9/0);  new base C (0/10);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.08e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGA‑CACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTG  >  NC_000913/1969647‑1969929                                                                                                                                                      |                                                                                                                                        ggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGA‑CACGATTGACTGTCTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCt                                                                                                                                      >  2:74534/1‑151 (MQ=255)                     cGATCATGGGCGCGGGGGGTTTATTCCCCAGTTTCGACCCCTGCGGGA‑CACGATTGACTGACTCTTCCATCCGCCCTTTGCTCTCTTTTGCTTCCTGCGCGCTCCGGCTGGCAAGATTGCGCCC‑TTCCCCTGCCACTACCGCAAAACCAc                                                                                                                   <  1:304071/150‑1 (MQ=255)                           tGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGA‑CACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                            >  1:30051/1‑151 (MQ=255)                           gggTCGCGGCGGCGTTTTTCCCCA‑TTTCAGACCCCGCTTGACCGCGATTGCCTCCTCCTCCGATCCCCCCTTTGTTCTCTTTTGCTGCCTGCGCGCTGCGGCTGGGAGGATTGGGCCC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                            <  1:355940/150‑1 (MQ=255)                                                ccAGTTTCGAACCCTGCTGGA‑CACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTc                                                                                       >  2:353297/1‑151 (MQ=255)                                                 catttccaaccccTGCTGGA‑CACGATTGACTGCCTCTTCGATCAGCCTTTTGATCTTTTTGGCGCCCTGCGCGCTCCGGCTGGAAAGATTGCGCCC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                      <  2:356955/141‑1 (MQ=255)                                                    tttCGACCCCGGCTGGA‑CCCGTTTGACTGCCTCTTCGATCACCCTTTTGATCTTTTTTGCTGCCTGCGCGCTCGGGCTGGCAAGATTGCGCAC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCAc                                                                                    <  2:69147/150‑1 (MQ=255)                                                        gAACCCTGCTGGA‑CACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                               >  2:28672/1‑151 (MQ=255)                                                              tgctgGA‑CACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                          >  2:2849/1‑150 (MQ=255)                                                               gctgGA‑CACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                          >  1:367450/1‑149 (MQ=255)                                                                           ttGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAatat                                                             >  2:63555/1‑151 (MQ=255)                                                                                            cgcccttttGATCTTTTTTGCTGCCTGCGCGCTCCGGCTGGAAAGTTTGCCCTC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAt                                             <  2:342699/144‑1 (MQ=25)                                                                                               cccTTTTATTTTTTTTCCTGCGGGCGCGCTGCGGCTGGCAAG‑TTGCCCCCTTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTgg                                                          <  1:181840/134‑1 (MQ=18)                                                                                                  tttgtcctctTTTCCCCCCTCCGCGCACCGGCTGGAAGGATTGCCCCT‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTTTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc                                           <  1:45055/140‑1 (MQ=21)                                                                                                   ttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                      >  2:261889/1‑150 (MQ=255)                                                                                                                               ggCTGGCAAGTTTGCGCCC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt           <  2:151005/149‑1 (MQ=255)                                                                                                                               ggCGGGCAAGATTGCGCCC‑TCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt          <  2:94484/150‑1 (MQ=255)                                                                                                                                  ggaaaaGTTTGCCCCC‑TTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGa      <  2:271712/147‑1 (MQ=255)                                                                                                                                      aaGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtg  >  1:164677/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                        GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGA‑CACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCAC‑TTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTG  >  NC_000913/1969647‑1969929

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