Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,554,013 | 0 | A | C | 38.6% | 51.9 / 11.7 | 44 | W553G (TGG→GGG) | maeA | malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme |
Reads supporting (aligned to +/- strand): ref base A (9/18); new base C (17/0); total (26/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.94e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTT‑CACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCGGAAGCAATA > NC_000913/1553780‑1554248 | cgcAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTCCGGCGGTAATCGCGGTGTTCGGCTTGCCACAAAt > 1:254202/1‑239 (MQ=255) gtTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTGGTCTCGGGGTTCTGCTTCCACGAAAtattttt > 2:1470757/1‑233 (MQ=255) aCGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTGGTCGCGGGATTCTGCTTGCCCGAAAATaatgttcaagt > 2:841316/1‑228 (MQ=255) acTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCCTGCCCGAAAATATCGGC‑AAT‑GGCCTGTTGCAGGGGTTCCGCAGG‑GGGTTTTAcc > 1:1559225/1‑238 (MQ=255) acTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGGTTTCCCGAAATTTTCTTC‑CAA‑GGGCCGTTGGAGGGTTTCGGCGGA‑GGGTTTCAcc > 1:697298/1‑238 (MQ=255) ttCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGACGCGGCATTCGGCCTGCCAGAAATTATTGCC‑AAT‑GGCGCGGTGCAGGGCGTTGGCCGG‑GGTGTTTTACTCCACGCCtc > 2:1248081/1‑238 (MQ=255) cATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCTGCTTGCCCGAAATTTTCGTC‑AAA‑GGCCTGTTGCAGGGCCCCGGCAGA‑GGTTTTTTCCGCCACCCCTtg > 1:467624/1‑238 (MQ=255) cctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCGTGGCCGGAAATATTGGC‑AAA‑GGGCCGTTGCCGGGGGTCGGCAGAGGGTTTTTACCGCC‑CGCCCTGCTTCTgtgc > 2:1167834/1‑237 (MQ=255) tCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAATTATCGTC‑AAT‑GGGCTGTTGCAgg > 1:1247913/1‑191 (MQ=255) cATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTTGCGGGATTCGGCCTGCCCGAAATTATTGTC‑AATGGGCCTGtt > 2:1557200/1‑172 (MQ=255) cAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGACGGGCGGGATTGGGGCTGCCACGAAATATTGTG‑AAT‑GGGTCGTTGGATGGCGTCCGCCGA‑GGGTTTTACTGGCAACCCCTGCTGTTGCGCCCTTCTGCCACTCGCaacgccactttg > 2:1118955/1‑229 (MQ=255) ttCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTTGCCGAAATTATCGTCAAAT‑AGCC‑GTTTCAGGGGTTCGGCCGA‑GGGTTTCACCGCCCCGCCCTGCTGCTTCGCCCCTTTTTCCACCGCCAAAGCAAATTTGCGGgggaagt > 2:1220937/1‑225 (MQ=255) ccAGTATGGATGAGATCTAAAAAGGGGAGAGGGAAATAGCCCGGTGGCCTTCCCTACCGGGCGCAGGCTTAGATGGGGGTACGGCGGTAGGCGCGGTATTCGGCTTGCCAGAAATTACCGTC‑AAT‑GGCCTGTTGCGGGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCgg < 2:787140/239‑1 (MQ=255) ttagggATGAGATCGAAAAAAGGGGGAGGTAAATACCCCGGTAGCCGTCACTACCGGGCGCAGGTTGAGAGGGAGGGAGGGCGTTAGTGGCGTTTTGCGGCTGCCAAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGTTTCCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAg < 2:1006960/235‑1 (MQ=255) aTGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTC‑aa > 1:341354/1‑115 (MQ=255) gggAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGTTTCGGCGTGGCCGAAATTATTGTC‑AAT‑gggc > 2:471558/1‑103 (MQ=255) agagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGTTTCGGCTTGCCCGAAATTATCATC‑CAT‑GGCCCGTTTCAGGGGTTCGGCAGAGGGTTTTCA‑CGCCACCCCCTGCTGTTTCCCCATTTTGTACAACGCCAACGCAAATTCGCGGGGGGCTTTCTGAATAACATTCCATTCCTGCAGTAACAATCCTTCGCCCTTCAGcc > 2:485945/1‑237 (MQ=255) tagccccgcccccGGTAGCCTTCACTCCCGGGCGCAGGTTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAG‑GGCCTGTTGCAGGGTTTCGGCAGA‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCa < 1:263388/229‑1 (MQ=255) agGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTTCCCGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGGTTTTACCGCCCCCCCTTGCTTCTGCGCCCTTTTTCCAACCGCCAACGCAAATTCGCGGGGGACTTTCTGTATATTTTTTAGTTCCGGCAGTACCATAACTTCCCCGTTCCGCACCa > 2:194072/1‑239 (MQ=255) ggAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAT‑GGCCTGTTGc < 1:248988/104‑1 (MQ=255) agcctgcgCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTTGGCTTGCCCGAAATTTTCGTC‑AAA‑GGGCCGGTTCAGGGGTTTGGCCGA‑GGGTTTTACCCCCACGCCCTTCTGCTGCGCCCCTTTTTCAACCCCAAAAGCAATTGCGCGGGGGACTTTCTTAATATCTTTTCGTTCCGGCAGGAACATACCCTCGCCCGTTACCACCaaa > 2:263388/8‑238 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGaaa < 1:36266/82‑1 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGaaa < 2:341354/82‑1 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑aa < 1:283866/92‑1 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑aa < 1:1557200/92‑1 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAT‑GGc < 2:1247913/96‑1 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTc < 1:841316/229‑1 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGGACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATc < 1:471558/87‑1 (MQ=255) aaaaTAGCCCGGTAGCCTTCACTCCCGGGCGCAGGTTTAGATGGAGGTACGGCGGTAGTCGCGGTTTTCGGTTTGCCAAAAATTATCGTC‑AAT‑GGCCTGTTGCGGGGTTTCGGCAGA‑GGTTTTCACCGCCACGCCTtgctgctg < 2:1504003/143‑1 (MQ=255) aaTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGtt < 1:1133901/201‑1 (MQ=255) aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTa < 1:206153/82‑1 (MQ=255) aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTa > 2:206153/1‑82 (MQ=255) aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGCTTGCCCGAAAATATCGTC‑AAA‑GGCCCGTTGCAGGGCTTCGGCAGA‑GGTTTTTACCGCCACGCCCTTCTGCTGCGCCCTTTTTCCAAACGCCAACGCACATGCGCGGGGGACTTTCTCAATATCTTTCAGTTTCGGCAGTCCCCTACCTTCGCCCGTCTGCACCCATgggg > 2:294139/1‑237 (MQ=255) tAGCCCGGTAGCCTTCACTACCGGGCGCCGGCTTAGATGGAGGTACGGCGGGTGTCGTCGGGTTCCGCTTGCCAGAAAATATCGTT‑CAT‑GGCCCGGTTGGGGGGTTCGGCGGA‑GGCTTTTCCCGCCCCCCCCTCCGGCTCTGCCACTTTTCCCACCCCAACCGCCAATTCGCGGGGGGCCTTCTTGATAGTTTTCTGCTCCCGCAGAAACATACCCTAGACTGTCAGCACCCATAGTgg > 2:857480/1‑238 (MQ=255) tAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCCGAAAATATCGTC‑AAT‑GGCCCGTTTCCGGGCTTTGGCCGA‑GGGTTTTACCGCCCCGCCTTGCTGCTGCGCCCCTTTTCCCACCGCCAACCCAATTTCGCGGGGGGCTTTTTGTATATCTTTTTGTTCCGGCAGTAc > 2:1133901/1‑201 (MQ=255) cACTCCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGTTTTCACCGCCa < 1:139054/112‑1 (MQ=255) cgcAGGCTTAGATGGAGGTACGGCGGTAGACGCCGGATTCGGCTTTCCCGAAAATTACGTC‑AAA‑GGCCTGTTGCCGGGGGTCGGCAGA‑GGGTTTTACCGCCACGCCCTTCTGCTGCGCCACTTTTCCAACCCCCAAAGCAACTTCGCGGGGGACTTTCTGAATATCTTTCAGTTCCTGCCGCAACAAAACTTCTCCGTTTATCACCCATGGTGAAAAATGCGCCAGCGGTC‑TACTTgg > 2:1328960/1‑237 (MQ=255) cgcAGGCTTAGATGGAGGTACGGCGGTAATCGCGGGATTCGGGTTGCCCGAAAATATTGTC‑AAT‑GGGCTGTTGCCGGGGTTCGGCAGA‑GGGTTTCACCCCCCCGCCTtgctg > 2:139054/1‑112 (MQ=255) ggaGGCTAAGATGGAGGTGCGGCGGTATCCGCGTTATTCGGCTGGCAAAAAATTATCGTC‑AAG‑GGCCTGTTGCAGGGCTTCGGCAGA‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGGCAGCGTTT‑CACTTGCCg < 1:1167834/237‑1 (MQ=255) aGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAATATTGGC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGGTTTTACCCCCCCGCCCTTCTGCTGCGCCATTTTTCCAACCGCCAACGCAATTGCGCGGGGGGCTTTCTGAATATCTTTTAGTTCCGGCAGTACCATACCTTCTCCGTTTCGCACCAATGGTGAATAATTCGCCAGCGTTTGTACTT‑CCGGCATCAGCAt > 2:175274/1‑235 (MQ=255) aGGTACGGCGGTAGGCGCGGTATTTGGCCTTCCCGAAAATATCGTT‑AAT‑GGGCTGTTTCAGGGGGTTGGGAGA‑GGGTTTTACCGCCACGCCTTTCTGCTGCGCCATTTTTCCAACCGCAAAAGCAATTTCGCGGGAGACTTTCTGAATATCTTTTACGTCCGGCAGGAACCTACCCTCGCCGCTTAGCACCCAAGGGGAATACTTCGCCCGCCGTT‑CACTTGCCCAcctatacctct > 2:880456/1‑227 (MQ=255) ggTGCGGCGGAGGCCGCGGTATTCGTCTGACAAGAAATTATCGCC‑AAG‑GGCCTGTTGCAGGGCTTCGGCAGG‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTT‑CACTTGCCGACATCAGCATCTCa < 1:352123/238‑1 (MQ=255) ccGGCTTGCCAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGCCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTT‑CACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCGGa < 1:294139/238‑1 (MQ=255) gCCAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTT‑CACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGACGGAAGCAATa > 2:113458/1‑239 (MQ=255) | CGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAT‑GGCCTGTTGCAGGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTT‑CACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGCGCCGGAAGCAATA > NC_000913/1553780‑1554248 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |