Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,158	0	C	T	45.5%	‑3.2 / 25.1	22	intergenic (+329/+407)	ycdU/serX	uncharacterized protein YcdU/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base C (3/9);  new base T (8/2);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTT  >  NC_000913/1096865‑1097444                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                               gTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGATGTGGGCGGGGGCACCCGGTGTGTTCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCGGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCACCCGGTGTGTGCAATATACGaac                                                                                                                                                                                                                                                                                          >  2:223426/1‑299 (MQ=255)           ctctTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCg                                                                                                                                                                                                                                                                                >  1:387671/1‑300 (MQ=255)                                    gAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGCCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTAATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCACGTACTTTCGTCCGGGCTCTTCTTTaa                                                                                                                                                                                                                                                       >  1:183150/1‑300 (MQ=255)                                                             gATTTCCCCTCGCCCTAACAGACGTGGGCGGGTGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAc                                                                                                                                                                                                                              <  1:277657/300‑1 (MQ=255)                                                              aTTTCCCCTCGCCCTAACACACGTGTGCGAGGGCACCCGGTGTGTGCAATATACNAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCGTTAAATAGGGCGGGGAGGGTGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTTGCGGGCGCCAGTTTTCGACCCCCGGCCGGGGCTTCTCATCCCCCCTGTGGGTGCAATACACGAAAAAAANGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACt                                                                                                                                                                                                                             <  2:681059/300‑1 (MQ=14)                                                              aTTTCCCCTCGCCCTAACAAACTTGGGCGGGGACACCCGGCGTGTGCAATATACGAAAAAAAAGCCTATACTTTCGTACGAACTCGTCTTTAAAGATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCGCACTGCGTTCACGGTCTGTCCAACCGGCTTTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACt                                                                                                                                                                                                                             <  2:21109/300‑1 (MQ=14)                                                                       cGCCCTAACAGACGTGGGCGGGGGCACCCGGTGGGTGCAATAGACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcttcgct                                                                                                                                                                                                                    <  1:332633/300‑1 (MQ=255)                                                                            tAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACCCGCTTCGCTCGCC‑CTGCGGTCAGCCCGCTCACTGCGTTCACGGTCCGTCCCACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAttagcccaccccatccgcca                                                                                                                                                                                                               >  1:705180/1‑282 (MQ=11)                                                                            tAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTAAACCCCCCAACCCGCCACAATATCCA‑CCCCCTCCCCTCATCCACATCCCCTCCCCCCAGCTACCGCCACTTGCCGACCCAAGGTCGGGGCTTCTCATCCACCAGGTGTGTGCAATATACGAAAAAAAAGAAAGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAATCGCTTAGCTCGccc                                                                                                                                                                                                               >  1:705215/1‑300 (MQ=9)                                                                                        gCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCTTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGATTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc                                                                                                                                                                                                   <  1:180753/300‑1 (MQ=38)                                                                                                                   cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAGTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAc                                                                                                                                                                        <  1:63481/300‑1 (MQ=17)                                                                                                                                           aCGAGCTCTTCTTTAAATATGGCGGTTAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCATTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAAccc                                                                                                                                                <  1:709018/300‑1 (MQ=17)                                                                                                                                                                           gggATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGTTCTGTCCAACTGGCTTTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGc                                                                                                                <  1:427894/300‑1 (MQ=17)                                                                                                                                                                                     cgcttcgctCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTTTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCACGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGaa                                                                                                      >  1:550843/1‑300 (MQ=17)                                                                                                                                                                                                                                                             gaaccccGGTCGGGGCTTCTCACCCCCCCTTTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTTGTTAAATATGGTGGTGAGGGGGGGGTTGACTTGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcttcgct                               <  1:445739/297‑1 (MQ=11)                                                                                                                                                                                                                                                              aCCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTAGGGGCTTCTCATCCCCCCGGTGTGTGTAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGaaccacctcgtgc                              >  1:166494/1‑288 (MQ=12)                                                                                                                                                                                                                                                                            cTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCAGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCACTCGCCCTGCgatccgc                >  1:111779/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                 cATCCCCCCGGTGTGTGCAATATACGAAAAAAACGCCCGTACTTTCGTATGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCCCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCAGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAACCCGTACTTTCGTACGAGCTCTTATTTCAATATGACGGTGAGGGGGGGATTAACCCGCTTCACGCGCCCTGCGAGCAGccctcc           >  1:623328/1‑297 (MQ=11)                                                                                                                                                                                                                                                                                 cATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCACGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTAGCTCGCCATGCGGGCAGCCCGCt           >  1:622918/1‑300 (MQ=16)                                                                                                                                                                                                                                                                                    ccccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc        <  1:398186/300‑1 (MQ=17)                                                                                                                                                                                                                                                                                          tgtgggcgCAATACACGAAAAAAAATCGCGTGCGTTCGTACAAGCTCTGCTGTGAGGATAGGGGGGAGGGGAAGATGGACTTGCTTCGCGCGCCCTTCTGGCGGCCCGCTCACTACTTCCACAATCTGTCCGATCGGCGGGCGCTAGTCTCCGACCCCCGGTCGTGGCTTCTTGCCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGtt  <  2:87855/293‑1 (MQ=2)                                                                                                                                                                                                                                                                                          ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTAACTCACCTCGCTCGCCCTACTGGCAGCCCTCTCACTGCGtt  >  1:315231/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                               GTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCC‑CTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTT  >  NC_000913/1096865‑1097444

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