Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 436,504 | 0 | A | C | 42.9% | 44.2 / 12.2 | 36 | T291P (ACC→CCC) | thiL | thiamine monophosphate kinase |
Reads supporting (aligned to +/- strand): ref base A (16/4); new base C (0/15); total (16/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.90e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGG > NC_000913/436364‑436651 | cgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCgg > 1:476751/1‑151 (MQ=255) ttCTCGCCATGTTGATCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGATCTGAACC‑GTGGAGCGCTGCATGTAGGTCTCGGACACCTGGGCGTACCGTTTCCC‑TGTATCGGGCa < 1:143786/150‑1 (MQ=255) ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCa > 1:620591/1‑150 (MQ=255) aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgc > 2:100775/1‑151 (MQ=255) gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgcc < 2:705857/150‑1 (MQ=255) ttGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGAccgccg < 1:518469/150‑1 (MQ=255) ggggCGCTCTGGGGGGGTAAATATTCCGTTTTGTTTTTCTCTGTCCGGAAATTGACCG‑GTGGCGCGGTGGGTGTGGCTCCCGTCCCCCGGGGCG‑CCCGTTTTCCTTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt < 2:126826/148‑1 (MQ=255) cgcTTTCGGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGACCC‑GTGGCGCGCTGGATGTGGCTCTCGGACACGTGGGCGTACCGTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg < 2:162554/150‑1 (MQ=255) gcTCTCTACCGGTAAAGAGTACGAGTTGTGTGTCACTGTGCCGGAACTGGACC‑GTTGCGCGGTGGATGTGGCTCTCGGACACCTGGGCGTACCATTTACT‑TGTATCGGGCACATGGCCACCGCGATCGAGGGGGTATGTTTGATTCgaggc > 1:78604/1‑147 (MQ=255) gggcgtggAAGTTTAGGAGTTGTTTTTCACTGTG‑CGGAACCGCACCGGGGGCGCGCTGGATGTGGCTCGCGGCCACCGGGGCGTCCTTTTTACC‑TTTAGCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg < 1:156201/139‑1 (MQ=255) gggggTGAAGTTTAGGAGTTGTTTTTCACGGTGCCGGACCTGACCC‑GTGGCGGGGGGGATGTGCCTCTCGGCCACCGGGGCGTCCCTTTTCCT‑TGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg < 1:210473/143‑1 (MQ=255) gAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt > 2:203421/1‑150 (MQ=255) gAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATCTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt > 2:438668/1‑150 (MQ=255) gATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt > 2:601221/1‑151 (MQ=255) ttcgggtttGTTTTTCACTGTGCCGGAACGGACCC‑GTGGCGCGCTGGAGGTGCCTCTCGGACCCCGGGGCGTCCCTTTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt < 2:711743/143‑1 (MQ=255) tgttTTTCACTGTGCCGAAACTGAACC‑GGGGCGCGTTGGATGTGGCTCTCGGACCCGTGGGCGTCCCTTTTCCC‑TGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa < 2:262518/150‑1 (MQ=255) tttCACTGTGCCGGAACTGACCC‑GTGGCGCGCTGGATGTGGCTCTCGGCCACCTGGGCGTCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg < 2:78260/150‑1 (MQ=255) cggaaccataaCC‑GGGGCACGCGGCTTGTGGCTCTCGGACACCGGGGATGCCCGTTTCCT‑TGTACCGGGAAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg < 2:123296/136‑1 (MQ=255) ccGGACCTGACCC‑GTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCTTACCTTTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa < 2:476841/146‑1 (MQ=255) ggAACTGACCG‑GTGGCGCGGTGGTTGTGGCTCTCGGACCCCGGGGCGTCCCGTTTCCC‑TTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt < 2:142538/151‑1 (MQ=255) aaaCTGACCG‑GTGGCGCGCTGGATCTGCCTTTCGTACCCCTGGGCTTCCCTTTTCCT‑TGTATGGGCCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt < 1:438576/147‑1 (MQ=255) aaCTGACCC‑GGGGCGCGCTGGTTGTGGCTCTCGGACCCCTGGGCGTACCTTTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc < 1:682532/151‑1 (MQ=255) aaCTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc > 1:252828/1‑151 (MQ=255) aaCC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc > 2:254399/1‑151 (MQ=255) gTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCT‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa < 2:595812/149‑1 (MQ=255) tGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa > 1:379767/1‑150 (MQ=255) gtggtgggggCTCCCGGACCCTGGGGCGTCCCTTTTCCT‑TGTAGGGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGa < 2:521351/139‑1 (MQ=255) ggttggggCTCGCGGACACGGGGGCGTCCCTTTTACC‑TGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt < 2:269716/140‑1 (MQ=255) tctcGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg > 2:602548/1‑151 (MQ=255) ctgggACCCCGGGGCGTCCCTTTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg < 2:294504/147‑1 (MQ=255) ctcGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc > 2:431227/1‑151 (MQ=255) acCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa > 1:702642/1‑148 (MQ=255) cgtGGGCTTCCCTTTTCCC‑TTTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa < 2:466203/148‑1 (MQ=255) cTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg > 2:198277/1‑151 (MQ=255) tttACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTgg > 1:52956/1‑151 (MQ=255) ttACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGt > 1:662437/1‑148 (MQ=255) | CGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGG > NC_000913/436364‑436651 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |