| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,567,738 | 0 | A | T | 42.4% | 40.9 / 17.7 | 34 | V362V (GTT→GTA) | yddW | liprotein, glycosyl hydrolase homolog |
| Reads supporting (aligned to +/- strand): ref base A (5/14); new base T (14/0); total (20/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.45e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCC > NC_000913/1567589‑1567887 | ccGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTaa < 1:795314/151‑1 (MQ=255) ttCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑CCTGGGCGGGTTTACGTCAtac > 2:392624/1‑149 (MQ=255) cGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑CCTGGTCGGGTTTAAGACATcccc > 1:381722/1‑149 (MQ=255) gggCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGGCGGGTTTAAGACATCCG‑cc > 1:779697/1‑151 (MQ=255) acaGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑CCTGGGCGGGTTAAAGACATCcc‑cccc‑cc > 1:695573/1‑151 (MQ=255) aGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGGGGGGTTTAAGACAAACG‑cccc‑cc > 2:847908/1‑149 (MQ=255) cATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGG‑CCTGG‑CGGTTTTACGAAAACCC‑CCCC‑CCC‑ttt > 1:263295/1‑151 (MQ=255) aGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAaca < 1:151929/151‑1 (MQ=255) gAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCa < 1:62266/151‑1 (MQ=255) gAGCTGC‑TTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGT‑CCGGGCGGGTTTAA‑GACATACG‑CCCC‑CCC‑TTTTTTCAacacacca > 2:49255/1‑146 (MQ=255) ctTTTTCAGTTCCGGTACGCCGCCGTCAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑GCTGGGGGGGTTTACGCCAAACC‑CCCA‑CCC‑CTTTGACCACACGCCCCAAcgc > 1:777006/1‑151 (MQ=255) tttCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTGCCTGGCCGGGTTCAC‑ACAAACGCCCCA‑CCC‑TTTTCCCaaaccccacaa > 2:641302/1‑136 (MQ=255) tCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCgg < 2:195029/150‑1 (MQ=255) tCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTCCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑CCTGGGCGGGTTTATGAAAACCC‑CCCA‑CCA‑CTTTTCCCacacaaccaaaagcgcgg > 1:685009/1‑136 (MQ=255) cAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑CCTGGTGGGGTTTACGAAACCCC‑CCCC‑CCC‑TTTTTCCCacacccccaaaaacgcgg > 2:733809/1‑135 (MQ=255) ccGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGGTTTACGACAACCC‑CCCC‑CCC‑TTTTTCCCACACGGCATAACGCGCGGCAcccc > 1:52320/1‑150 (MQ=255) tAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑CCTGGGCGGGTTTAAGACAACCC‑CCCC‑CCC‑CTTTTTCCACACCTCATAAAGCGCGGCCCTCCCTGGGAACgg > 1:491724/1‑145 (MQ=255) gTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCCCGGCGGGTTTAAAACAACCG‑CCCC‑CCACTTTTTCCAA‑ACATCATAAAACGCGGCACCCCCCGGGAACGGcccatgaaaataggg > 2:625739/1‑137 (MQ=255) ttGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGG‑CCTGGTCGGTTTTAAGACAACCG‑CCCACCCA‑TTTTTCC‑ACACCCCATAAAGCGCGGCACTCCCCGCGAACGGCCAGGAAAATTTGGGGGCAAAGTAAttc > 2:780816/1‑147 (MQ=255) tGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAg < 1:41765/150‑1 (MQ=255) gAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGc < 1:172511/150‑1 (MQ=255) aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAAt < 2:776970/151‑1 (MQ=255) aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGA‑CCTGGTCGGTTTTACCACCTTCC‑CCCC‑CCA‑CTTTTCCCACACGCCATAACACGCGGGCCCCCCTGAGAACGGCCAGTAAAATTTGGGGGGAATATaaactcacaa > 1:103804/1‑142 (MQ=255) cccACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGtt < 1:158201/151‑1 (MQ=255) ataGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCAt < 1:436425/151‑1 (MQ=255) ataGAAGGCGATACCGATATACAGGC‑GGGG‑CCCGGTCGGGTTTAAGACAACCC‑CCCC‑CCC‑CTTTTTCACCACGCCATAAAGCGCGGCGCCCCCTGTGAACGGCCCGCAAAATTTGGGGGGAATGTAATCCACCAATCCTTTTTTCAccc > 2:179144/1‑149 (MQ=255) gATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGgtg < 2:363002/151‑1 (MQ=255) gATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCg < 2:25756/151‑1 (MQ=255) ataCAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAg < 2:46761/151‑1 (MQ=255) cAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGAt < 2:429085/151‑1 (MQ=255) c‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGt > 1:820221/1‑151 (MQ=255) gT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCAt < 1:610528/151‑1 (MQ=255) cGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGATGcc < 1:899246/151‑1 (MQ=255) tAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGcc > 2:370278/1‑151 (MQ=255) tAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGcc > 1:151277/1‑151 (MQ=255) | CCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGT‑CCTGGTCGGTTTAACGACATCCG‑CCCA‑CCA‑TTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCC > NC_000913/1567589‑1567887 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |