Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,421,443 | 0 | A | C | 50.0% | 39.0 / 16.0 | 38 | Y40S (TAC→TCC) | folX | D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase |
Reads supporting (aligned to +/- strand): ref base A (14/5); new base C (0/19); total (14/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.41e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC > NC_000913/2421298‑2421592 | ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTAcccc < 2:551356/150‑1 (MQ=255) aaGTAAAACAACGGAGACCCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCccgccg < 1:507544/151‑1 (MQ=255) aaGTAAAACAACGGAGAACCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGATGATCCACTACCccgccg > 1:406422/1‑151 (MQ=255) tAAAACAACGGAGAACCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGAt > 2:579559/1‑150 (MQ=255) aaaaCAACGGAGAACCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGGAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATa > 2:55064/1‑150 (MQ=255) ggaaacccTGAAAGGCCACAC‑CCTGCCGTTTTTTTTCTTAAAAAGACCCTTC‑GTTTGCGTACTTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTCAGGA‑TATTGTAATCAATGTGACAACCCACTCCCCCGCCGATAAAgcgcgc < 2:571810/144‑1 (MQ=255) acccTGCAGTGCCACAC‑CCTGCCGTTATTATTCGTAAAAAGACCCTTC‑TTTGGCGTCCGTTTATGGAAATAAAGAAAGAAGAAATTACCAACCGTCAGAA‑TATTTTTATAAAGGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTa < 1:523803/148‑1 (MQ=255) aaCCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt > 2:408957/1‑149 (MQ=255) ccTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGCATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACACCGCCGATAAAGCGCGCACTAGCg > 2:237845/1‑151 (MQ=255) aaggccacaC‑CCTGCCGTTATTATTGGAAAAAAGACCCTCC‑TTTTGCGTACTTTTACCGAAATAAAGGAAGAAAAAATAAACACCCGGAAGAA‑TTTTTTTAACAGTGGGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGa < 1:71814/142‑1 (MQ=255) agggaccaacccggCCGTTATTATTCTAAAAAAAACCCTTC‑TTTTGCGTTCGTTTACGGAAATAAAGAAAAAAAAAATAAACACCCGTAGGGA‑TTTTTTTACAAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGa < 1:579574/137‑1 (MQ=255) aatgCACAACCCTGCCGTTATTATTAATA‑AAAAAACCTTCGGTTTGCGTCCGTTTATAGAAATAAAGAAAAAAAAAATAAACACCGGCCAGGA‑TATTTTTAAAAATGTGACGATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAg < 2:411014/146‑1 (MQ=255) cacaC‑CCTGCCGTTATTATTCTAAAAAAGACCTTTC‑TTTTGCTTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTAAGGA‑TTTTTTAACAAAGGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat < 1:353709/148‑1 (MQ=255) ccTCCCGTTTTTTTTCGTATAAAAACCTTTC‑TTTTGGGTATTTTTAGAGAAATAAAGAAAGAAGAAATAACCACCCGTCAGGA‑TTTTGTTAAAAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa < 1:413911/145‑1 (MQ=255) cTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa > 1:48551/1‑150 (MQ=255) ttataCGAAAAAAGACCCTCC‑GTTTGGTTGCTTTAATCGGAATTAGGGAAGAAAAAAAAAACCCCCTTCGTTA‑TTTTTTTACCGATGCGGCGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg < 1:387023/146‑1 (MQ=255) attCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc > 2:372259/1‑151 (MQ=255) ataAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg > 2:225660/1‑151 (MQ=255) aaaaaGATCCTTT‑TTTTGCATACGTTAATCGGAATAAGGAAAAAAGAAATAAACAACCTTGGGAT‑TTTTGTTAAAAATGGGCCGTCCCACTACGCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt < 1:426251/142‑1 (MQ=255) taAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt > 1:177353/1‑151 (MQ=255) aaacccTTC‑GTTTGCTTCCTTTTATCGAAATAAAGAAAAAAAAAATAACCACCCGTAAGAA‑TTTTTTTATAATTGTGACTACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc < 1:113030/147‑1 (MQ=255) aGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc > 1:317741/1‑151 (MQ=255) cTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa > 2:10676/1‑150 (MQ=255) ttgtgcgtttttttATCGAAATTAAGGAAAAAAAAATAAACACCCGCCAGAT‑TATTGAAAACAATGTGCCAATCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa < 1:135962/134‑1 (MQ=255) tttgggatgtttttTCGGAATGAAGGAAGAAAAAAAAAACACCCG‑AAGGATTGTGGTAATCATTGTGATGATCTACGCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaaa < 2:393576/137‑1 (MQ=255) tttGCGTATTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCTTCGGAT‑TTTTGTTATAATTGGGCCGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa < 2:114389/144‑1 (MQ=255) ttGCTTCCGTTTATGGAAATAAAGGAAAAAAAAATAAACACCCTTCAGAA‑TTTTTTAATAAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt < 1:555789/150‑1 (MQ=255) gtgttcttttATCGAAATAAAGAAAGAAAAAATAAACACCCGTAAGGT‑TTTTGTTAACAATGTGCCATCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt < 2:471743/144‑1 (MQ=255) tttATCGAAATAAAGAAAAAAAAAATAAACACCCGTAAGTA‑TTTTTTTATAAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg < 2:317734/147‑1 (MQ=255) ttacggaaaTAAAGAAAAAAAAAATAAACACCCGTCAGAT‑TATGTTTAAAAATGTGCCCACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg < 1:93633/143‑1 (MQ=255) tgaaaaTAAAGGAAAAAAAAATTAACAACCGACATGA‑TTTTTTTATCAATGTGACGACCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt < 1:417981/144‑1 (MQ=255) aaTTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt > 1:264510/1‑151 (MQ=255) taaagaaagaAAAAATTAACACCCGTAAGAA‑TTTTTTTATCAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa < 1:147516/148‑1 (MQ=255) aaGAAAAAAAAAATAACCACCCGTGAGGA‑TATTTTTATAAAGGGGCCCACCCACCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc < 1:225665/151‑1 (MQ=255) aaaTAAACATCCGACAGAT‑TATTGTTAAAAATGAGACGCTCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct < 1:576582/149‑1 (MQ=255) aaTTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTCCTCTTTg > 1:60113/1‑151 (MQ=255) tttttttATAAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt < 2:391857/139‑1 (MQ=255) tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc > 2:379552/1‑151 (MQ=255) | CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAA‑CCTGCCGCTATTATTCGTATAAAGAACCTTC‑GTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGA‑TATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC > NC_000913/2421298‑2421592 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |