Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,953,550 | 0 | A | C | 56.7% | 21.5 / 42.2 | 30 | T25P (ACC→CCC) | ilvD | dihydroxyacid dehydratase |
Reads supporting (aligned to +/- strand): ref base A (1/12); new base C (17/0); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.50e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-06 |
CAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTG > NC_000913/3953408‑3953688 | cAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATg < 1:1431/151‑1 (MQ=255) taaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGAcc < 1:316547/151‑1 (MQ=255) cAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCCGGATGACCGACGc > 2:553854/1‑150 (MQ=255) aaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGAttt < 2:615222/151‑1 (MQ=255) ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGGGCGCTGTGGCGCGCCCCCCGGAGGACCGACCCCGAGTTCGGTAAgg > 1:389320/1‑150 (MQ=255) ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGGGCGCCGTTGGGCGCCCCCCGGCTGACCGACGCCgcgttctggaagg > 1:297305/1‑139 (MQ=255) ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCCCCGGCGCCCTGGCGCGCCCCCGGCAAGAACGCCCCCGCGTTTGGTAAg > 1:614294/1‑150 (MQ=255) gcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCCCGTGCGCTGTGGCGCGCCCCCGGGATGACCGACGCCGATTTTGGTTAGGCGGtt > 2:633961/1‑151 (MQ=255) cGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAAc < 1:552997/151‑1 (MQ=255) gacaCTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGGGGCGCGCCCCCGGGAGGAACGACGCCGGTTTTGGTTAGCCGGTTATCGCGGGTGGGGAATCGTTc > 1:153656/1‑150 (MQ=255) acacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGtt < 1:654612/148‑1 (MQ=255) cacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGGGCGCCCCCCGGCTGACCCACCCCGATTTCGGGGAGCCGGTTATTGCGGTTTTGGAATCGTTCAc > 2:46673/1‑150 (MQ=255) gAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTGGGGCGCCCCCCGGAAGACCGCCCCCCCTTTCTGTAAGCCGAGTATCGCGGGTGTGGAATCCTTCcccccaatt > 1:72306/1‑146 (MQ=255) aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTAc < 2:488108/151‑1 (MQ=255) gTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCCTTGGCCGTGGCGCGCCCCCCGAAAGGACGACCCCGCGTTCGGTAAGGCGATTATCGCGGGTGGGGACTCGCTCCCCC‑CATTTTTACCGGGGCa > 1:452898/1‑151 (MQ=255) ccTAAGTACCGTTCCGCCGCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGCTCACGTc < 1:416237/150‑1 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTGGGCCGTGGGGGGCCCCCGGAATGAACGACGCCGATTTTGGGGAGCCGGTTATTTCGGGTGTGGAATCGGTTACCC‑CAATTTTACCGGGGCGCGTc > 1:657779/1‑150 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTGCGCTGTGGCGCGCCCCCCGGATGAACCACCCCCCTTTTGGTAAGCCGATTTTTGCGGTTGTGAACTCCTTTCCCC‑CATTTTTACCGGGGCACGTc > 2:532925/1‑150 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCCGGAAGACCCACGCCCGTTTCGGTTAGCCGGGTATCGCGGGTGTGGAATCGTTCTCCC‑CATTTTTAACGGGGCACGTc > 1:207465/1‑150 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCCGCGCTGCGGCGCGCCCCCCGGACGACCGACGCCGCTTTCGGTAAGCCGATTATCGCGGCTGTGGACACCTTCACCC‑AAATTGTACCggggtacacc > 2:386450/1‑143 (MQ=255) ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGAGGCTCGGGCGCTGGGGCGCGCCCCCCGGATGACCGACGCCGGCTTTGGGGAGCCGGTTATCGCGGGTGTGAAATCGCTCCCCCAAATTTTTCCCgggcccggc > 1:294062/1‑144 (MQ=255) cTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCa < 2:14850/151‑1 (MQ=255) gTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGTGGCGCGCCCCCGGAATGACCGACGCCCGTTTTGGTAAGCCGGTTTTCCCGGGTGGGAAATCGTTCACCC‑CACTTTTTCCgggggcccgccccct > 1:181914/1‑138 (MQ=255) tACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTgc < 1:368730/151‑1 (MQ=255) ccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACt < 2:75130/151‑1 (MQ=255) ccacTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCGGAATAACCGACCCCGCTTTTTGGAAGACGGTTATCGCGGGTGTGGAATCCTTCACCC‑AATTTTTATCCGGGCACCTCCCTCTGCGCGATCTCggggaa > 2:485104/1‑145 (MQ=255) cacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTg < 2:27413/151‑1 (MQ=255) acTCATGGTCGTAATATGGCGGGGGGTCGGGCGCTGTGGCGCGCCCCCCGGAAGAACGAAGCCGCTTTTGGGGAGGCGGTTTTCGCGGGTGTGGAATCGGTCACCC‑CCCTTTTTTCGGGTCCCGGCCATCCGCGCGGtctccgggaactcg > 1:55695/1‑141 (MQ=255) ggCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGaa > 2:455836/1‑151 (MQ=255) cTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTg < 2:292309/151‑1 (MQ=255) | CAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTG > NC_000913/3953408‑3953688 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |