| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 41.1% | 74.7 / 19.8 | 56 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (30/3); new base A (0/23); total (30/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.21e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01 | |||||||||||
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACT > NC_000913/436367‑436659 | tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGc > 1:404322/1‑150 (MQ=255) ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGc < 2:436043/149‑1 (MQ=255) tcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAAt > 2:933979/1‑150 (MQ=255) tcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTAGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATg > 2:925593/1‑151 (MQ=255) cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑aaa < 2:488523/148‑1 (MQ=255) tGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGAccgcc > 2:779750/1‑151 (MQ=255) tGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGAccgcc > 1:266578/1‑151 (MQ=255) tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGat > 2:397242/1‑151 (MQ=255) gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGat > 1:1338000/1‑150 (MQ=255) cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCg > 2:1354312/1‑150 (MQ=255) gAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGa > 1:142270/1‑149 (MQ=255) ggCGCTGCGCTGGGCGCTCTCGGGGGGTGAAG‑ATTACGAGTGGTGTTTCCCTGT‑GCCGGAACTGACCCGTGGCGCGGTGGATGTGGCTCGCGGCCCCGGGGGCGTCCCTTTTCCCTTTATG‑GGGA‑AAATGACCGCCGATATCGAAGGGc < 1:780311/149‑1 (MQ=255) tgcgctgGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGttt > 2:1366706/1‑151 (MQ=255) gcgctgGGCGCTCTCGGGCGGTGAAG‑ATTACGATTTGTGTTTCACTGT‑GCCGGAACTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATC‑GGGA‑AAATGACCGCCGATATCGAAGGGCttt < 1:960277/146‑1 (MQ=255) ctgGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt > 2:153616/1‑151 (MQ=255) tgGGCGCTCTGGGGGGGGAAAG‑ATTAGGGTTTGTTTTCCCTTGT‑GCCGGACCTAACCGGTGGCGCGTTGTTTGTGGCTCCGGGACGCCGGGGCGTCCCTTTTCCCTGTATC‑GGGA‑AAAGGACCGCCGATATCGAAGGGCTTTGtttt < 1:605655/147‑1 (MQ=255) tgGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt > 1:607031/1‑150 (MQ=255) gggCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg > 1:688908/1‑151 (MQ=255) gggCGCTCT‑GGGGGGGGAAGAATTCCGCGTTGGGTTGCACTGT‑CCGGGAACTGACCCGTGGGGCGCTGCATGGGCCCCTCGGCCAGCGGGGCGTCCTGTTTACTTTTTTC‑GGCA‑AAATGACCGCCGATATCGAAGGGCTTTGttt < 1:423386/145‑1 (MQ=255) tctctGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg > 2:932913/1‑151 (MQ=255) tctGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACggtg > 2:1294550/1‑149 (MQ=255) gggggTAAAG‑TTTACGAGTTGTGTTCCACTGT‑GCGGGACCTGACCCGTGGCGGGCTGTATTGGGCTCTCGAACCCCGGGGCGTCCCTTTTCCTTTTTTG‑GGGA‑AAATGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc < 1:1312918/144‑1 (MQ=255) ggCGGTGAAT‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTACCTGTAGG‑GGGA‑AAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg < 1:1069699/145‑1 (MQ=255) ggCGGTGAAG‑ATTGCGGTTGGTTTTTCCTTGT‑CCCGGACCTGACCGGGGGCGCGTTGGTTGTGGCTCCCGGACCCGTGGGCGTTCTTTTTCCCTGTACG‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa < 2:1325037/149‑1 (MQ=255) ggTGAAG‑ATTACGAGTTGTTTTTCACTGG‑GCCGAAACTGACCGGGGGCGCGTTGGTGGTGGCTCTCGGCCCCCGGGGCGTCCCTTTTCCCTGTATG‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt < 1:19525/150‑1 (MQ=255) aaG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc > 1:1151363/1‑151 (MQ=255) aG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc > 2:921964/1‑150 (MQ=255) ttaggtgTTGTGTTTC‑GTGTGGCCGGACCGGCCCCGTGCGGCGCTGGAGGTGCCTCTGGGCCACGTGGGGGTCCTTTTTCCCTTTATG‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc < 2:1286215/141‑1 (MQ=255) ttACGAGTTGTGTTTCACTGT‑GCCGAAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGCCACCTGGGCGTCCCTTTTACCTGTATC‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt < 1:1178424/150‑1 (MQ=255) aGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg > 1:942560/1‑151 (MQ=255) aGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGAGGTGGCTCTCGGACAACTGGGCGTACCGATTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTCTTCGTGACGGCGAACCTGTTACATTAGACTg > 2:236984/1‑151 (MQ=255) ttGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg > 2:831473/1‑150 (MQ=255) tgtgTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa > 1:134051/1‑150 (MQ=255) ttgttcccttgt‑gCCGAAACTGACCCGGGGGGCGTTGGAGGTGCCTCCCGGCCACCTGGGCGTCCTTTTTACTTTTACC‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa < 2:539879/140‑1 (MQ=255) ttCACTGT‑GCCGAAACTGACCCGTGGCGCGCGGGATGTGGCTCTGGGACCCTGGGGCGTCCCTTTTACCTGTATC‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa < 1:832052/150‑1 (MQ=255) cattgt‑gCTGGAACGCACCCGGGGCGCGTTGGAGGTGGCTCTCGGACCCGTGGGCTTGCTTTTTCCTTTTACC‑GGGA‑AAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata < 2:92254/147‑1 (MQ=255) cACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat > 2:881199/1‑151 (MQ=255) tgt‑gCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc > 1:772809/1‑151 (MQ=255) gt‑gCCGGAACTGAACCGTGGCGCGTTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTCCCTTTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat < 1:981343/145‑1 (MQ=255) gctgaaacgaaaCCGTGGCGCGTTGGTTGGGGCTTCCGGCCCCTGGGGTGAGCCTTTCCCCTGTATC‑GGCA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa < 1:1151170/137‑1 (MQ=255) ccGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt > 2:427244/1‑151 (MQ=255) ggAACTGAGCCGTGGCTCGCTGGTTGGGTCTCTCGGACCCCGGGGCGTCCCTTTTTCTGGTGGC‑GGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt < 2:840274/151‑1 (MQ=255) gAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg > 2:256781/1‑151 (MQ=255) gacccGGGGCGCGTTGGATGTGGCCCTCGGACACCTGGGCGTCCCGTTTCCCTGTATC‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa < 1:1007442/146‑1 (MQ=255) ggCTCGGTGGAGGTGCCTCCCGACCACCTGGGTGT‑CCTTTTTCCTTTTTCGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa < 2:144482/149‑1 (MQ=255) tcTGGTTGTGCCTCCCGGACCCCTGGGCGTCCCTTTTCCTTGTACC‑GGAA‑AAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt < 1:703970/148‑1 (MQ=255) ggttgtgCCTCTCGGACCCGGGGGCGCCCCTTTTCCTTTTAGC‑GGAA‑AAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt < 1:375865/145‑1 (MQ=255) cGGAAACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt > 2:640725/1‑151 (MQ=255) cTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg > 2:279132/1‑151 (MQ=255) tGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg > 1:1219028/1‑150 (MQ=255) ggggGTCCCGTTTCCCTGTATC‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGa < 1:1000442/150‑1 (MQ=255) ggggtccGTTTTCCTTGTTC‑GGGCAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc < 2:1152692/134‑1 (MQ=255) gtcccGTTTCCCTGTATC‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTaa < 1:332705/147‑1 (MQ=255) ccGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg > 1:263009/1‑151 (MQ=255) ttttCCCTGTAGG‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc < 1:969943/146‑1 (MQ=255) ttACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc > 1:987125/1‑151 (MQ=255) ccTTTATG‑GGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCa < 1:1104121/149‑1 (MQ=255) tATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt > 1:1069814/1‑151 (MQ=255) | TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAG‑ATTACGAGTTGTGTTTCACTGT‑GCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATC‑GGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACT > NC_000913/436367‑436659 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |