| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,097,299 | 0 | A | G | 50.0% | 17.5 / 45.3 | 42 | intergenic (+470/+266) | ycdU/serX | putative inner membrane protein/tRNA‑Ser |
| Reads supporting (aligned to +/- strand): ref base A (5/16); new base G (3/18); total (8/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.97e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01 | |||||||||||
GCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCA > NC_000913/1097152‑1097446 | gCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGccc < 2:4682/151‑1 (MQ=14) acacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTc < 1:99482/151‑1 (MQ=14) acacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTc < 3:68898/151‑1 (MQ=14) cacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCg < 2:1624816/151‑1 (MQ=11) cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTg < 2:25563/149‑1 (MQ=11) cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTg < 2:1058228/149‑1 (MQ=11) cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑Tg‑ggggggggtgcgcacccc > 1:810979/1‑140 (MQ=11) cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCGCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACt < 1:1357518/147‑1 (MQ=17) gTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑Gg‑gggggggggttggaccccc > 2:981232/1‑139 (MQ=1) ttCACGGTCTGTCCAACTGGCTGTCGCCATTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACtc < 1:562378/146‑1 (MQ=18) tCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACtcgctt < 2:594949/150‑1 (MQ=9) cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTGAGtcccct > 2:48806/1‑144 (MQ=9) aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGG‑GGGGGGGGTTTGCtcgcttcgc > 2:1718930/1‑151 (MQ=11) cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTGCGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTtcgctc < 2:336188/151‑1 (MQ=18) cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGTGGT‑GA‑GGGGGGGATTGACTCGCTtcgctc > 1:2656/1‑151 (MQ=18) cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGAAATATACG‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTtcgctc < 1:1222970/150‑1 (MQ=16) ggTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTGACTCCCCccgctc > 2:1533363/1‑150 (MQ=255) gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTtcgctcg < 2:1008746/151‑1 (MQ=9) gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTGAcacccccccctccc > 4:560079/1‑138 (MQ=255) gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTtcgctcgc < 3:18647/151‑1 (MQ=33) gTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGcc < 4:432720/149‑1 (MQ=9) gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTg < 2:866453/151‑1 (MQ=33) tCCAACTGGCTGTCGCCTGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTtcgctcgc < 1:1040810/147‑1 (MQ=9) ccAACTGGCTGTGGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGc < 2:115679/151‑1 (MQ=9) aaCTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg < 3:528189/151‑1 (MQ=9) aaCTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCggg < 2:1130624/151‑1 (MQ=14) cTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GA‑GGGGGGGGTTGACTCGCTTCGCCCGCCCCGCCGGcc > 2:979253/1‑150 (MQ=255) ggCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa < 2:1409130/149‑1 (MQ=33) gCTGTCGCCAGTTGTCGAACCCCGGACGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTACGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa < 2:1221839/148‑1 (MQ=18) gTCGCCAGTTGTCGGGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc < 1:1282882/150‑1 (MQ=9) gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCg < 2:1017470/150‑1 (MQ=33) gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGc < 3:272871/151‑1 (MQ=33) cGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc < 1:1699087/148‑1 (MQ=9) cGCCAGTTGCCGAACCCCGGTCGCGGCTTCTCGCCCCCCCGGTGTGTGCAATATTCGAAAAAAAAGCCCTTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTc < 1:473732/151‑1 (MQ=11) ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc < 2:669107/151‑1 (MQ=33) aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt < 2:152099/151‑1 (MQ=9) aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt < 1:1547391/151‑1 (MQ=9) aGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt < 1:1043576/151‑1 (MQ=9) cGAGCCCCGTTCGGGGCTTCTCATCCCCCCGGTGTTTGCAATATACGAAAAAAAAGCCCGTTCTTTCGTTCGAGCTCTTCTTTAAATATGGCGGT‑TAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc < 3:114423/151‑1 (MQ=11) cGAGCCCCGCTGGGGGGTTCCGGTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGGCCGTTCTTTCGGACGAGCTCTTCTTTAAAAATGGCGGT‑GAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGt < 1:701590/149‑1 (MQ=11) cGAACCCCGGTCGGGGCTTCTCATCCCCCCGTTGTGTGAAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTCCTTTAAATATGGCGGT‑AA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCa < 1:467533/151‑1 (MQ=255) cGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc < 2:350876/150‑1 (MQ=33) | GCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCA > NC_000913/1097152‑1097446 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |