Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,798	0	T	C	37.5%	20.6 / 10.6	33	intergenic (+405/+310)	eco/mqo	serine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (5/15);  new base C (2/10);  total (8/25)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTT  >  NC_000913/2304552‑2305035                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAAGCCGCATCCCACATCCAACGGCCgg                                                                                                                                                                                                                                             >  1:499304/1‑250 (MQ=16)        ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                                                                                                                                                                                                                                      <  1:306829/251‑1 (MQ=37)          gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCCACGTCCCAGCCGGGTG‑c                                                                                                                                                                                                                                   >  1:44412/1‑251 (MQ=255)                   cgcccagccccgTTGCCTGATGCGACGCGGGCGCGTCTTAACAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCAAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGa                                                                                                                                                                                                                          <  2:22565/240‑1 (MQ=11)                                 gcgggaggccaCGCTGGCGCGTCTTTACAGGACTCAACCGCTGTGACGTGCCCCACCCCGTAGGTCTGATAAGCCGGGTACGCCGCCCCCGACATCGAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGAGAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCtt                                                                                                                                                                                                            <  2:354055/241‑1 (MQ=11)                                 gccggATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTTTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTGACGCCGCATCCGACATCTAACGCCCAAGCCGTTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTTCACCGCTGTGAATTGCTCCACCCCGTAGGTTGGATAAGGCGGCTACGCCGCATCCGACATGCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCtt                                                                                                                                                                                                            <  1:365747/247‑1 (MQ=11)                                    tGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCCTCCGACGTCCCACGTCCGAGCCGGTTG‑CCTGGTGCGACGCTGGcgcggcctatt                                                                                                                                                                                                         >  1:18255/1‑244 (MQ=255)                                    ggctGCGCCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                                                                                         <  2:59438/248‑1 (MQ=21)                                        gCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAgg                                                                                                                                                                                                      <  1:129716/250‑1 (MQ=35)                                          gACGCTGGCGCGGCTTATCAGGCCTACACCGCTGAGAAGGGCCCCCCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTGACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                                                                                                                                                                   <  2:44412/251‑1 (MQ=255)                                                   gcgTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAACCCGTAGGTTGGATAAGGCGGTTTCGCCGCATCCGACATCCCAAGTCCGGGCCGGTTGGCCT‑ATGCGACCCTGGCTCGTCTTTTCCAGCCTCCACCGCtg                                                                                                                                                                                          >  1:173433/1‑251 (MQ=11)                                                                     acCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTTCGCCGCATCCCAAATCCAACGTCCGGGCCGGTTG‑CCTGGTGCGACGCTGGCGCGTCTTATCCGGCCTACACCGCTGTGAAGTGCCCCAcccccc                                                                                                                                                                        >  1:179578/1‑249 (MQ=14)                                                                             gAAGTGCTCCACCCCGTAGGTCGGATAAGTCGGTTACGCCGCATCCGTCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTGTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                                                                                                                 <  1:496044/250‑1 (MQ=255)                                                                               aGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCCTCCCACATCCAACGCCCGAGCCGGGTG‑CACGCTGCGCCGCTGGCGCGTCTTATCAGGCCCACACCGCCTTGGAGTGCTCCACCCCCTAGGGCgggt                                                                                                                                                               >  1:146835/1‑248 (MQ=11)                                                                                   ccgcacccctcggggcgggtagggCGGTAACGCCGAATCCGACATCTAACGCCCGAGCCGGGTGCGTGATGCGACGGTGGCGCGTCTTATCAGGCCTACGCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                                                                                                                          <  2:140630/229‑1 (MQ=11)                                                                                         acagtgggtggtaTAAGGCAGATACCCCGCACCCGGCATCTACCACCCCATCCCGTTACCTGATGCGACGCTCGCGCTGCTTATCAGGCCTACACCGCGGTGAAGCGCGCCACCCCGTAGGTGGGTGAAGGAGGTTACGCCGACCCCGACATCCAACCCCCGAGCCGGTTG‑CGTGATGCGACGCTGCCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc                                                                                                                                                    <  2:173433/239‑1 (MQ=2)                                                                                                  cggctAAGGCGGTTCCGCCGCATCCGACATCTAACTCCCCAGCCGGTTGCCTGATGCTACGCTGGCGCGTCGTAGCAGGCCCACACCGCTGTGAAGTGCCCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc                                                                                                                                           <  2:499304/247‑1 (MQ=14)                                                                                                  cGGAAAAGGCGGTCCCGCCGCAGCCGACAGCCAACGCCCGAGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATc                                                                                                                                            <  2:50083/250‑1 (MQ=11)                                                                                                         ggCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCCACATCCCACGCCCGAGCCCGTTG‑CCTGATGCGACGCCGGCGCGTCTTATCAGGCCTTCCACGCTTTGAAGGGCCCCACCCCGTAGGTCGGATAAGGCGGTTaacccccctccctcatct                                                                                                                                    >  1:482722/1‑234 (MQ=14)                                                                                                                                 caccGCCCGAGCCGTTTGCCAGATGGTGCGCTGGCGCGTCTTGTCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCTGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTG‑CCTGATGCGCCGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATg                                                                                                           <  2:81606/248‑1 (MQ=11)                                                                                                                                    cccccgagccgggtgctggATGCGTCGCTGGCGCGGCTTATCAGGCCTACACCGCGGTGAAGTGCTCCTCCCCGTAGGGTGGATAAGGCGGTGACGCCGCATCCGACATCCAACCCCCGAGCCGGTTG‑CCGGATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGa                                                                                                        <  2:211495/233‑1 (MQ=11)                                                                                                                                                  gCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGGTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCCCATCCAACATCTAAACGCCCA‑GCCGGGTGCCTGATGCGACGCCGGCGCGGCtt                                                                                          >  1:9478/1‑251 (MQ=11)                                                                                                                                                           gACGCTGGCGCGTCTTATCAGTCCTACACCGCTGTGAAGTGCTCCACCCCGTATTTTGGATAAGGTGGTTACGCCGCATCCGACATCCAACGTCCGAGCCTGTTG‑CCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTTAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                                                 <  1:538502/251‑1 (MQ=14)                                                                                                                                                                gggcgcgtcggagcgggccgaagccGCTGTGAAGTGAGCCACCCCGTAGGTGGGATAGGGCGGTTACGCCGCATCCGACTGCCAACGCCCGAGCCGGGTG‑CCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAcc                                                                            <  2:50997/228‑1 (MQ=11)                                                                                                                                                                      gTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGCCATCCAACGCCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtgtg                                                                      <  2:474238/251‑1 (MQ=255)                                                                                                                                                                          tagcAGGCCTACACAGCTGTGAAGGGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAGCGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                   <  1:355077/247‑1 (MQ=255)                                                                                                                                                                                           cctatcttGCTCCCCCCCGTAGTTGGGATAGGGCGGGTCCGCCCCATCCTACATCCGATGTCCGCGCCGGTTG‑CCGGATCCGTCTCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCGAA‑CGCCCAAGCCGGGTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCGGTGAAGTGCTCCACCCCGTAGGtt                                                 <  2:215463/244‑1 (MQ=11)                                                                                                                                                                                                aGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTTTGAAGTGCTTCACCCCGTAGGTTGGATa                                            <  1:466591/251‑1 (MQ=255)                                                                                                                                                                                                       cACCCCGGAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCGTACACCGCTGTGAAGTGCTCCACCCCGTAGGGCGGATAAGGCGGTTACGCCGCGTCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGGCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCTCCTCGTAGGTTGGATAAGGCGGt                                     <  1:484213/251‑1 (MQ=14)                                                                                                                                                                                                             ggagggggcaTAAGGCGGTTACGCCCCATCCGACATCCACCGCGCGACCCGGTAG‑CCTGACTCGACGCCGGCACGTCGTACCTGGGCTCCACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCGGATGCGACGCGGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgcc                               <  2:179578/242‑1 (MQ=11)                                                                                                                                                                                                                               gtACGCCGCACCCGACACCCAACGCCCGAGCCGGTTG‑CCTGACGCGACGCTGGCGCGACTTATCAGGCCTACACCGCTGTTAAGTGCGCCCCCCCTGAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCGAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAGCAGGCCTACACCGCTGTGAAGTGCTCCTCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGc             <  2:229078/250‑1 (MQ=11)                                                                                                                                                                                                                                   gccgcATCCGACATCAAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGa         <  1:201890/251‑1 (MQ=37)                                                                                                                                                                                                                                          ccGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACACTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGtt  <  1:449135/251‑1 (MQ=37)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                GCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTG‑CCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAA‑CGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTT  >  NC_000913/2304552‑2305035

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