Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,376,505 | 0 | A | C | 42.9% | 86.7 / 19.6 | 56 | T158P (ACC→CCC) | ycjR | putative TIM alpha/beta barrel enzyme |
Reads supporting (aligned to +/- strand): ref base A (24/8); new base C (0/24); total (24/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.75e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT > NC_000913/1376355‑1376650 | cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAaca < 1:1059868/151‑1 (MQ=255) ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAacac < 2:393594/149‑1 (MQ=255) tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACAcc < 2:60230/149‑1 (MQ=255) aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACAccc < 2:574436/149‑1 (MQ=255) cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGc < 2:458886/151‑1 (MQ=255) gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGcc > 2:71456/1‑151 (MQ=255) ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCg > 1:1009979/1‑151 (MQ=255) gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGAt > 1:299850/1‑151 (MQ=255) ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCccg < 1:433921/151‑1 (MQ=255) gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAC‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATAATATGATCAACACCCTCTCCGATGCGCGCCGTTAc > 2:539386/1‑151 (MQ=255) tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc > 2:270199/1‑149 (MQ=255) ggACGGCGCCCGAAAAAGGGTGAGTGATCCCCGGCGCGTACGGAAACAGGCCGCCGCGCGTCCGGGAC‑CGGTGTTTTATCCCAACCTTTTAAC‑CCC‑C‑TATAAGGATCATATAAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGaa < 1:246679/151‑1 (MQ=255) ggCGACCGAAAAGGGGGGGTTGTTCCCCGGCGCGTACGGAAACGGGTCCCCGCGCGTACCGAAC‑CTGTGGTGTATCTCGACCCGTTAAC‑CCT‑C‑TATCAGGTTCATTTGATCAACACCCCCGCCGATGCCCGCCGTTACAtcgtc < 2:61414/144‑1 (MQ=255) ggCGACCAAAAAGTGGTGAGTGTTCCCGTGCGCTTACTGAAACAGGTCGCCGCGCGTACCGGAC‑CGGTGGTGTTTCTC‑AACCCTTAAACCCG‑C‑TACCAGGATCATATGACCAACCCCCTCGCAGATCCCCGCCGTTACATCGTCGAAAACg < 2:494400/151‑1 (MQ=255) aaaaggggTGAGTGATTCCCTGCGCTTCCGGAAACGGGCCCCCGCGC‑TTCCGGAACCCGTGTTGTATCCCG‑ACCCTTAAA‑CCGCC‑TATCAGAATAATATAAAAACCCCCCGGGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa < 1:295993/144‑1 (MQ=255) aaaaTGGTGAGTGCTCCCCGGCGCGTACGGAAACGGGTCGCCGCCCGTACCGGAA‑CGGTGGTGTTCCTCGACCCTTAAAC‑CCG‑C‑TACCAGGATAATATGACCACCCCCCTCGCCAAGGCCCGCCGTTACATCGTCGAAAACGATCTTa < 1:739458/148‑1 (MQ=255) aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc > 1:441991/1‑151 (MQ=255) aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc > 2:478845/1‑151 (MQ=255) aaaaGGGTGAGTTATCCCCGGCGCTTCTTGAACCAGGTCCCGGCGCGTCCCGAAA‑CGGGGGTGTATCT‑GAAACTTTAAACCCG‑T‑TATCGGTATCATAGGATCAAACCCCCCGCGGAGCCCCGCCGTTACATCGTCGAAAACGATCTTaa < 2:316349/149‑1 (MQ=255) tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt > 2:168262/1‑151 (MQ=255) ggggagtgaTTCCCGGCGGGTCCGGGAAAAGGTCGCGGCGCGTACCGGAC‑CCGTGTTGTATCCCGAACCTTTAAA‑CCG‑T‑TATCAGTTTAATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa < 1:187148/148‑1 (MQ=255) gaTTCCCGGCGCGTACGGAAGCAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAC‑CCG‑T‑TATCAGGATCATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt < 1:462221/150‑1 (MQ=255) aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt > 1:807054/1‑151 (MQ=255) cccTGCCCGTACGGGAACGGTCCGCCGGGCGTCCCGAAC‑CCGGGGTGTACCTCAAACC‑TTAAA‑CCC‑CTTATAAGGACAATAGGAACACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt < 2:861185/148‑1 (MQ=255) gcgcTTACGGGAACGGGTCGCGGCGCGTCCGGGAC‑CGGTGGTGTATCTCGACCCTTTAAC‑CCG‑T‑TATCAGGTTCTTATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg < 1:521170/146‑1 (MQ=255) gcgTACTGGAACAGGTCGCCGCGCGTCCCGGAA‑CCGTGGTGTATCTCGACCCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt < 1:303794/151‑1 (MQ=255) ggAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac > 1:94164/1‑151 (MQ=255) aCAGTTCCCCGCGAGTCCCGGTG‑CTGTGGTCTATCCCGCACCATTAAG‑CCG‑C‑TAACGGCATCATTTTACCAGCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt < 1:200742/150‑1 (MQ=255) ggTCGCCGCGCGTACCGGAA‑CCGTGGTGTACCCCGAACCTTTAAC‑CCG‑C‑TACCAGTATCATATGACCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa < 2:883575/149‑1 (MQ=255) cgcgcgTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga > 1:179052/1‑150 (MQ=255) gcgTCCCGGAA‑CCGGGGTGTATCCCGACCCGTAAAC‑CCG‑C‑TACCAGAATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa < 1:630470/151‑1 (MQ=255) gtcccGGAC‑CGGGGGTGTTTCCCAAACCTTTAAC‑CCG‑A‑TATAAGGATCATATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga < 2:1027272/145‑1 (MQ=255) aCCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt > 2:451707/1‑150 (MQ=255) cGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc > 1:867865/1‑151 (MQ=255) ggaccggtggtgTACCCCGACCCGTAAAA‑CCG‑C‑TATAGGGTTCATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag < 1:644204/139‑1 (MQ=255) aa‑CCGTGGTTTATCCCGACCCGTAAAC‑CCG‑C‑TATCAGAATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc < 1:515169/148‑1 (MQ=255) tttttcccaaaCCGTTAAC‑CCC‑T‑TATCAGAATAATATAACCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg < 2:247383/140‑1 (MQ=255) tgtATCTCGACCCTTTAAC‑CCC‑C‑TAGCAGAATAATATAAACAACACCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg < 1:510171/146‑1 (MQ=255) tgtACCCCGAACCTTTAAC‑CCT‑T‑TATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg < 1:295879/150‑1 (MQ=255) tctcGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg > 2:175361/1‑151 (MQ=255) cccgaCCATTAAACCCG‑C‑TATCAGAATAATATAACCAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg < 2:616037/135‑1 (MQ=255) caacccGTTAAC‑CCG‑C‑AATCAGAATCATATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg < 1:451718/144‑1 (MQ=255) cGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt > 2:692078/1‑151 (MQ=255) gAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg > 1:58826/1‑151 (MQ=255) cctttAAC‑CCG‑T‑TACCGGGATCATATGACCACCCCCCCCGCCGAGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc < 1:10284/147‑1 (MQ=255) ccGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc > 1:264112/1‑150 (MQ=255) ccGTAAAC‑CCG‑C‑TACCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa < 1:432980/151‑1 (MQ=255) tAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg > 2:1050152/1‑151 (MQ=255) cccttaccAGGATCATATAACAACCCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg < 1:665333/140‑1 (MQ=255) c‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg > 2:887756/1‑151 (MQ=255) tATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc > 1:594812/1‑151 (MQ=255) aTCAGGATAATATAAACAACCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt < 2:994172/148‑1 (MQ=255) aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg > 2:905641/1‑151 (MQ=255) aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg > 2:942239/1‑151 (MQ=255) ggATAAAATGACCAACCCCCTCGCGGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt < 2:222139/144‑1 (MQ=255) cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg > 1:975080/1‑151 (MQ=255) tGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc < 1:409078/150‑1 (MQ=255) tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat > 1:549800/1‑151 (MQ=255) | CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCG‑C‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT > NC_000913/1376355‑1376650 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |