Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,329,894 | 0 | A | G | 45.1% | 70.6 / 12.3 | 51 | D188G (GAC→GGC) | sohB | inner membrane protein, S49 peptidase family protein |
Reads supporting (aligned to +/- strand): ref base A (7/21); new base G (23/0); total (30/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTG > NC_000913/1329745‑1330043 | gCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGa < 1:555534/150‑1 (MQ=255) ccGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑a < 1:401863/149‑1 (MQ=255) cTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTATTGGCTGG‑GTTGCGGGC‑aaaa > 2:721835/1‑148 (MQ=255) ggAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGt < 1:95088/150‑1 (MQ=255) aaGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCAATATTGGCCGT‑GTGGGGGGC‑CAAAATtttt > 1:419237/1‑149 (MQ=255) aGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTATGGCCCGT‑GTGGGGGGC‑AAAAATttttt > 1:236034/1‑148 (MQ=255) tggtggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGT‑GTGGCGGGG‑AAAAATTTTTTCG‑CAcc > 1:676960/1‑150 (MQ=255) tggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTATATAAAGGGCTCT‑GTGGGGGGG‑AAAAATGTTTTCC‑CAcccctt > 2:10368/1‑147 (MQ=255) tggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGCCTGT‑GTGGGGGGC‑AAAAATGTTTTCG‑CAcc > 1:779249/1‑147 (MQ=255) ggCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTg < 1:284274/151‑1 (MQ=255) catggtgcatggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAACACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTATGGCCCGT‑GTGGGGGGC‑AAAAATTTTTTCT‑CAcc > 1:440369/1‑144 (MQ=255) catggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTCATTAATGTCCCT‑GGGGGGGGC‑AAAAATGTTTCCGCCCCCG‑TTTgttttttgg > 2:38736/1‑143 (MQ=255) catggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑gg < 2:647144/151‑1 (MQ=255) ggtTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGtt < 1:224688/151‑1 (MQ=255) cTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTTT‑GTTGGGGGG‑AAAAATTTTTTCG‑CACCC‑CTTGCtt‑tt > 2:512672/1‑135 (MQ=255) cTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGCCTGT‑GTGGGGGGG‑AAAAATATTTTCC‑CACCC‑GTTTTTA‑TTTTGGGGTTCCTTTgcg > 1:291270/1‑149 (MQ=255) tGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTAAAGATGGCCCGT‑GTGGGGGGG‑AAAAATTTTTTCC‑CCCCG‑GTTGCTATTTTT‑GGGTTccttgggg > 1:655808/1‑144 (MQ=255) tCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGgc < 1:100563/149‑1 (MQ=255) gCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGT‑GTGGGGGGG‑AAAAATTTTTTCG‑CACCC‑GTTGCTATTTTT‑GGGTTCCTTTggggggggggg > 2:303850/1‑139 (MQ=255) cAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑aaa < 1:208758/151‑1 (MQ=255) aaCTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGT‑GTGGGGGGC‑AAAAATTTTTTCG‑CCCCC‑CTTTCTT‑TTTTGGGGTTCCTTTgcggggggggggcaa > 1:425345/1‑136 (MQ=255) tGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTAAATTATTGGCCGT‑GGGGGGGAC‑AAAAATTTTTTCC‑CCCCG‑GTTTTTT‑TTGT‑TGGTTGCAATTGG‑GT‑GGGGGCG‑G‑aaaaacc > 1:196607/1‑147 (MQ=255) cAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGGTACATGAAGGTCTCT‑GGGGGGGGA‑AAAAATTTTTTCG‑CACCG‑CTTTCTC‑TTTT‑GGGTTCCATTTGC‑GG‑GGGGGGG‑C‑gaaa > 2:992073/1‑145 (MQ=255) cGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACt < 1:930378/151‑1 (MQ=255) gTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTAATGGCTGTGGGGGGGGGA‑AAAATTTTTTCCC‑CCCCT‑TTTGTTT‑TTTG‑GGGTTTC‑TTTGC‑Gt‑gggggggggaaagacc > 2:1040749/1‑135 (MQ=255) gTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGATATCTTATGGGCTTT‑GTGGGGGGCGAAAATTGTTTCC‑‑CACCC‑TTTTCTA‑TTGT‑GGGGTTCTTTTGC‑GG‑GGGGGGG‑T‑GAAA‑GCCCCACtt > 2:382643/1‑151 (MQ=255) gcgcGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTaa < 2:173059/150‑1 (MQ=255) gcgATAAAAACATTCCTTTAACTGTGACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CAGCG‑TTTGCTA‑TAGT‑GGCGTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTa < 1:489360/147‑1 (MQ=255) aTAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TCTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCtt < 1:228599/151‑1 (MQ=255) tAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTAAATGATGGGCTGT‑GGGGCGGGC‑AAAAATGTTTTCC‑CCCCC‑GTTGCTA‑TTTT‑GGGGTTCATTTGG‑GT‑GGGGGGG‑CGAAAT‑CCCAACCTTTAACCCCt > 1:71286/1‑150 (MQ=255) aaaaaCATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTcc < 1:623226/151‑1 (MQ=255) cATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGGCTGT‑GGGGGGGGC‑AAAAATGTTTCTC‑CACCC‑TCTGCTA‑TTTT‑GGGTTTCCTTTGT‑GT‑GGGGGCGCC‑AAAT‑CCCCACCTTTaccacatcctgg > 1:345768/1‑138 (MQ=255) tCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTTAATTATAGCCCGG‑GTGGGGGGC‑AAAATTGTTTTCC‑CCCCC‑TTTTCTT‑TTTT‑TGGTTTCCTTGGCGGG‑GGGGGGG‑C‑AAAA‑CCCCACTTTTAACCCCTTCCGGaaa > 1:578055/1‑149 (MQ=255) cTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGGCTGT‑GTGGGGGGG‑AAAAATTTTTTCG‑CACCC‑CTTTCTA‑TTTT‑GGGGTTCCTTGGC‑GG‑GGGGGGG‑C‑AAATACCCCAATTTTAACCCCtcctcgaaa > 1:392306/1‑139 (MQ=255) aaCTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGat < 1:10369/151‑1 (MQ=255) gTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATGGCCCGG‑GCGGGGGAG‑AAAAATTTTTTTC‑CCCCG‑GTTTTTA‑TTGT‑GGGGTTCATTGGC‑GT‑GGGGGGG‑C‑AAAA‑CCCCCACCTTTTCCCCCTCCTTAAAAACAAAGaaa > 1:1041853/1‑146 (MQ=255) gTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTg < 1:394945/151‑1 (MQ=255) tACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATTATTGCCTGG‑GTGGGGGGA‑AAAAATGTTTTCC‑CCCCGCTTTTTTT‑TTTT‑GGGTTCCTTTGCG‑GG‑GGGGGGC‑C‑AAAG‑GCCCACCTTT‑ACCCCCTCCTTAAAAACAAAAAAAtt > 2:387363/1‑148 (MQ=255) cAAAGTCGCTGCCAGCGGCGGGTATATTATGGGCTGT‑GTGGCGGAC‑AAAAATGTTTTCC‑CACCC‑GTTGCTT‑TTTT‑GGGGTTCATTGGG‑GTGGGGGGGG‑C‑AACT‑ACCCAACTTAAACCCCTT‑CTTAAAAGCCAAGATATTGata > 2:722722/1‑143 (MQ=255) cGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCacac < 1:465242/150‑1 (MQ=255) cGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACt < 2:625823/151‑1 (MQ=255) cTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGGGGGG‑AAAAATTTTTTCC‑CACCC‑CTTGCTA‑TTTT‑GTGTTTCATTTGC‑GG‑GGGGGGG‑C‑AAAA‑ACCCCACTTTTATCCCCTCCCGAAAAAAACAGAAATTTATATctca > 2:796023/1‑137 (MQ=255) cAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCggg < 1:444462/151‑1 (MQ=255) gcggcgGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCa < 1:1007927/151‑1 (MQ=255) cggcggTTACATGATGGCCTGT‑GTGGGGGGC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTTCATTGGC‑GG‑GGGGGGG‑C‑AAAA‑GCCCAACTTTTAACCCGTCCTGAAAAGAAAAAATATTGATATCGAAATGCACACTGCCggg > 2:951926/1‑148 (MQ=255) cggTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGtat < 2:917352/151‑1 (MQ=255) ggTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGtata < 1:818598/151‑1 (MQ=255) ggTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGtat > 1:994385/1‑150 (MQ=255) ggCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTg > 1:44694/1‑151 (MQ=255) ggCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTg < 2:680997/151‑1 (MQ=255) cGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTggg > 1:784363/1‑151 (MQ=255) ggAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTggg > 1:810437/1‑150 (MQ=255) aC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTg > 2:183154/1‑150 (MQ=255) | GCTCCGTCTGGAAAGCCCTGGTGGCATGGTGCATGGTTACGGGCTGGCGGCTTCGCAACTGCAGCGTCTGCGCGATAAAAACATTCCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGT‑GTGGCGGAC‑AAAATTGTTTCCG‑CACCG‑TTTGCTA‑TTGT‑GGGTTCCATTGGC‑GT‑GGTGGCG‑C‑AAAT‑GCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTG > NC_000913/1329745‑1330043 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |