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

GGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCT  >  NC_000913/1969646‑1969939                                                                                                                                                       |                                                                                                                                                  gggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CAcc                                                                                                                                                 >  2:307275/1‑151 (MQ=255)     cACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGcc                                                                                                                                             >  1:338574/1‑151 (MQ=255)                            ggttCGCGGGGGCGTTATTCCCCAGTTCCAACCCCGGCGGGACGCGTTGGGGTGACTTTCCGACCACCCCTTGGCTCTTTTTTGCTCCCGGCGCATTCCGGCT‑GCAAGTTTTGGCCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                       <  1:465689/147‑1 (MQ=255)                                 gcggcggctTTATTCCCCATTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCGGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGctc                                                                                                                   <  1:544594/149‑1 (MQ=255)                                 gcggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCtctc                                                                                                                 >  1:250299/1‑151 (MQ=255)                                  cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                                >  2:392513/1‑151 (MQ=255)                                             ttCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGACCAGCCTTTTGATCTTTTTTGCTGCCTGCGCGCTCCGGCTGGAAAG‑ATTGCGCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                                      <  2:259459/150‑1 (MQ=255)                                                   aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                >  2:141497/1‑150 (MQ=255)                                                             ccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                     >  1:424595/1‑151 (MQ=255)                                                               ggctgAACCCGTTTGACTGCCTCTCCGACCGCCCTTTTGTTCCCTTTCGCCCCCGCCGCGCTCCGGTGGACAAG‑GTTGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                     <  1:432972/148‑1 (MQ=255)                                                                   ggACACGTTTGACTGCCTCTTCGACCACCCTTTTGTTCTTTTTTGCTCCCTGCGCGCTCCGGCTGGAAAG‑TTTGCGCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg                                                                                 <  2:199922/149‑1 (MQ=255)                                                                           ttGCCTTACTCTCCGACCCCCCCTTTGTTCTTTTTTCCTGCCTGCGCGCTCCGGCGGGCAAG‑TTGGCCCACTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                                         <  1:508227/149‑1 (MQ=21)                                                                           ttGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAatat                                                                       >  1:27149/1‑151 (MQ=255)                                                                           ttGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAatat                                                                       >  2:232635/1‑151 (MQ=255)                                                                                      ttCGACCAGCCTTTTGTTCTTTTTTGCGGCCGGCGCGCTACGGCGGGAAAG‑ATTGCGCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCt                                                                 <  2:52359/146‑1 (MQ=25)                                                                                            cAGCCTTTTGTTCTTTTTTGCTCCCGGCGCGCTCGGGCGGGAAAG‑ATTGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCaa                                                        <  2:577218/149‑1 (MQ=21)                                                                                               cccTTTGATCTTTTTTGCTCCCTGCGCGCTGCGGCGGGCAAG‑ATTGCGCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc                                                     <  2:289443/149‑1 (MQ=25)                                                                                                         ttttttGCTCCCTGCGCGCTCCGGCTGGCAAG‑TTTGCGCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                             <  2:225040/145‑1 (MQ=37)                                                                                                           ttttGCTGCCTGCGCGCTCCGGCAGGAAAG‑TTGGCCCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                        <  2:478052/150‑1 (MQ=37)                                                                                                           ttttGCCCCCTGCGGGCGCCGTCGGGAAAG‑TTTGCGCACTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg                                          <  2:188705/148‑1 (MQ=21)                                                                                                            tttCCTGCCTCCGGGCTGCGGCTGCAAAG‑ATGGCGCTCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                       <  2:499842/150‑1 (MQ=25)                                                                                                                ctccctgcGCGCTACGGCTGGAAAG‑ATTGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATaa                                   <  2:20539/147‑1 (MQ=37)                                                                                                                   cctgcGCGCTCGGGCTGCCGAG‑TGTGCGCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATaa                                   <  1:456326/147‑1 (MQ=25)                                                                                                                            accggcggGCAAG‑‑TTGCCCCCTTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                      <  2:209516/149‑1 (MQ=25)                                                                                                                               ggCTGGAAAG‑TTTGCCCCTTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt                    <  2:485963/150‑1 (MQ=25)                                                                                                                               ggCGGGCAAG‑ATTGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtt                        <  2:321773/146‑1 (MQ=255)                                                                                                                                            tGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGc       >  2:288318/1‑151 (MQ=255)                                                                                                                                                 aCTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtct  >  2:127758/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                  GGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAG‑ATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCT  >  NC_000913/1969646‑1969939

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