Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,256,061 | 0 | C | T | 33.3% | 27.1 / 10.3 | 30 | E136K (GAA→AAA) | plsB | glycerol‑3‑phosphate 1‑O‑acyltransferase |
Reads supporting (aligned to +/- strand): ref base C (3/17); new base T (10/0); total (13/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.52e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATTAAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAAC > NC_000913/4255812‑4256302 | tattaACCAGACCCTGACGGCCAGGCAAACGTGGGCTGCCGGCGGCCAGACGTGGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCGGAGCGATAGTTTTACCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑c < 2:570946/250‑1 (MQ=255) attaAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCGCCTTTTT‑ccc > 1:187169/1‑249 (MQ=255) taAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCGGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACg < 1:407162/250‑1 (MQ=255) acggacccggACGAGCAGGCAAACGGGGGCCTCCGGCGCCCATGCGTTGACGGGCAAAGGGCATCCGCGCCACGCGCCCAAGTTTCTGAGCGATAGTTTTATCCGGGCC‑GTGTCCACCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCGCCTTTTT‑CACGcccc < 2:552507/242‑1 (MQ=255) gATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTT‑TACGcccccgc > 1:521545/1‑248 (MQ=255) acgggggCCTCCGGCCGCCAGCCGGTGACGGGCAAAGTGCAAACGCGCCGCGGGCGCCAGTTGCGGAGCGATAGTTTTATCCGTGCA‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAG‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCg < 2:511305/247‑1 (MQ=255) gTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCATCTCGCCTTTTTTCGCGCCCC‑gctggacgacaaacaaccccccc > 2:263845/1‑229 (MQ=255) ggCCTAAGGCAGCCAGACGTTGACGGGAAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTACTTTTCCCTTTT‑TTCGCCCC‑cgcgcgcggcccaaaacccacccaacc > 1:547949/1‑232 (MQ=255) gCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCCCGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑g < 1:306537/251‑1 (MQ=255) aGCCAGACGTGGCCGGACAAAGTGCATCCGCGCCGCGCGCGCCAGTTTTTGAGCGATATTTTTA‑CCGGGCCCGTGTCCTTCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACGGTCGCTCCCGAGCCGCAGTACAGCGAAAAAG‑TTCTGGACTCCGTGAAGCAGACGCGGCGGCGGTTTCACTTCGCCTTTTT‑CACTCCCC‑GGCGCGCGACCAAACTTC‑ACCGACACT‑GGCACCATCt < 1:661302/251‑1 (MQ=255) cgcgCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTT‑TCCGCCCC‑CGCGCGCGACCAAACATCCACCG‑CACT‑GGCACCATCTTCACCACCAGAATTTGGGT‑GGT‑ACGGTg > 1:530069/1‑251 (MQ=255) aCGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCCTTTT‑TAAGCCCC‑CGCGCGCGACCAAAAAAC‑ACCGACAAT‑GGGACCCTCTTCACATCCCGA‑CTTggggtggctacgtggcaaa > 1:524468/1‑234 (MQ=255) ggTTTCTGGGCTATAGTTTTACCCGTGCC‑G‑GTTCACCCCGCCATACGGCGCAGCGAAACTGACGGCGAAAAACGCACAAAACTTTGGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGTTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTc < 2:684802/250‑1 (MQ=255) aGTTTCTGAGCGATAGTTTTATCTGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTc < 1:635383/251‑1 (MQ=255) gAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCCCTTT‑CACGCCCC‑GGCGCGCGACCCAACACC‑ACCCACACT‑GGCACCATATGCCCATCCACA‑TTTTGGGTGCT‑ACGGTGCAAATCGAGAAAGTCGTgggac > 1:267239/1‑248 (MQ=255) gaaaGTTTTACCCGTGCG‑GTGTCCACCC‑GCCATACGGCGCACCGAAACGGCCGGCGAGAAACACAAAAAACTGGCGCGACCGAGCCACAGTCCAGCAAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGAAAGTCGTGGAACAGCt < 2:227868/248‑1 (MQ=255) tCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTT‑TACGCCCC‑CGCGCGCGGCCCAACATC‑ACCGCAACC‑GGCGCCCTCTGCACATCCCGA‑TTTTGGGTGCT‑ACTGGGCAAAACGAGATAGTCGTGGGACAGCTTAATATActc > 1:213452/1‑251 (MQ=255) gtgcc‑gtgTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtcttc < 1:134411/251‑1 (MQ=255) cgcgcgtgttcccaccgcccaccGGCGCAGCGAAACGGACGGCGAGAAACGCACAAACCTGGCTACCCCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtctt < 2:225497/229‑1 (MQ=255) ggcctggctccacccggcaggcggtgcAGCGAAACTGGCGGCGAGGACAGTAAAAAGCGGTCGCGGCCTAGCCACAGTCCAGC‑AAAAATATTCTATCCGCCGTTGAGCATACGCGGCGGCGGGTTCACTTCGCCTTTTG‑CTCCCCCC‑GGCGCGCGACCAAACATC‑ACCGGCACT‑GGCACCATCTGCACATCCAGA‑TTTGGGGTGCC‑ACGGTGCAAATTGAGAGAGTCGTGGATCGGCTTAATAGACtcttc < 2:660434/226‑1 (MQ=255) ggtatcagccc‑GCCATCGGGCCCAGCGAAACTGAGGGCGAGAAACGCAAAAACCTGTCGCACCCGGGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTc < 2:1609/242‑1 (MQ=255) gtgtTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCTCCCTTTT‑TACGCCCC‑GGCGCGCGGCCAAACATC‑ACCCACACC‑GGCACCCTCTGCACAACCAGA‑TTTTGGTTTCTTAGCTTGCAAA‑CAGGATAGTCGGGGAACAGCTTAAAATACtctttcttcc > 1:643962/1‑245 (MQ=255) gCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTTGCCTTTTT‑TACGCCCC‑CGCGCGCGCCCAAACATC‑ACCGCCACC‑GGGCCCCTCTGCACAACCCCA‑TTTGGGGTGCT‑TCGGTGCACATCGAGGTAGTAGTGGAGAAGCTTAATAGaattctcctttggcgtgtgat > 2:251662/1‑230 (MQ=255) ggCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTAACTTCTCCTTTTT‑TACCCCCC‑GGCGCCCGGCCAACCATC‑ACCCACCCT‑GGCCCCCCCTTCACATCCCGC‑TTTTGGGTGCT‑ACGGGGCAAATACAGATAGTCGTGGGACAGCTTTATAGACTCTTCTTTTCGCGTGTAAAAGGGGCAAACACGCGGAccgccctg > 2:262788/1‑248 (MQ=255) acaGTACAGCTAAAAAT‑TTCTGTACGCCTTTAAGCATACGCATCGGCGGGTTCACTGCGCTTTTTT‑CACGCCCC‑TGCGCGCGTCCAAATAAC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGATACAAATCTGGATAGTCg < 2:640706/155‑1 (MQ=255) caGTACAGCGACAAAT‑TTCTGTACGCCGTTACGCATACGCAGCGGCGGGTTCACTTTTACTTTTT‑TAATCCCC‑GGGGCGCGACCACACATC‑ACCGCCCCC‑GGGAACCTCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGGTAGTAGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAGAACACGCGGCCCGCCGTGGATGAACACAAAGCGCGGCGGGAGCGTggcgcggg > 2:442667/1‑244 (MQ=255) gtccGCCGTTAAGCATACGCGGCGGGGGGTTCCCTGCGCTTTTTT‑CACGCCCC‑GGCGCGCGAACAACCATC‑CCCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGAAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAAc < 1:253173/245‑1 (MQ=255) gccacttaaagagcccgcgagggggggTTCACTTCGCTTTTGG‑CACGCCCG‑GGCGCGCGCCAAAACACC‑CCCGACCCT‑GGCGCCATCGGCACTTCCAGT‑TTGGGGGTGCT‑GGGGTCAAAAGGGATATAGTGGGGGAACAGCTTAAGAGACTTGTTTTTGGGCGGGTAATGGGTAAACCCCCGGGCCCCGCCGTGAATGAACACAGAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTc < 2:312791/227‑1 (MQ=255) cGCCGTTAAGCATACGCAGCGGCGGGTTCACTTTGCCTTTTT‑TAACCCCC‑GGCGCGCGACCACACACC‑ACATACACTCGGC‑CCATCTGCCCATCCAAACTTTGGTTTCAT‑ACGGGCCAAAT‑GAGCTCGACTTGGGACAGCTTAAATGACTCTTCTTTTTGCGGGTTATAGGTGAAAACACCCGGCCCGCCCTGCATGGACACATAGAGCGGCAGTAGCGGGGCGGCGATTTTCAGCGGccataaccgggcc > 1:316277/1‑239 (MQ=255) gCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACCCCCC‑GGGGCGCGACCAAACCAC‑ACCCACACT‑GGCCACACCTGCACATCCAGC‑TTTGGGGTGTT‑ACGGGGCAAATAGAGATAATCGTGGAACAGCTTAATAGACTATTCTTTCGGCGGGGAAAAGGGGAACACACGCGGCCCGCCGCGAAAGAAAACAAAACGCGGCGGTAGCGTGCCGTCGATTTTCAGCGGCCCTCACggggc > 1:680599/1‑247 (MQ=255) gcggcggcgGTTTCCCTTC‑CCTTTTT‑CCCCCCCCCGGCGCCCGACCAAACACC‑CCCGACACG‑GGCCCCATCTGCACATCCAGA‑TTTGGTTTGCT‑ACGGTGCAAATGGAGATAGTCGTGGAACAGCTTAATAGACTCTTTTTTCGTCGTATAATAGGTGAACACCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGcc < 1:449172/248‑1 (MQ=255) cGCCTTTTC‑CCCGCCCC‑GGCGCGCGACAAAACATC‑CCCGACCCT‑GCCACCATCTGCACATCCAGA‑TTTGGTTTGCT‑ACGTTGCAAATCGAGATAGCCGTGGAACAGCTTAATAGACTCTTCTTCCGGCTTGTAAGAGGTGAACACACGCGTCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAAc < 1:58789/251‑1 (MQ=255) | TATTAAACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCC‑GTGTTCATCC‑GCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAAT‑TTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GGCGCGCGACCAAACATC‑ACCGACACT‑GGCACCATCTGCACATCCAGA‑TTTGGGTTGCT‑ACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAAC > NC_000913/4255812‑4256302 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |