| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,538,109 | 0 | A | C | 49.1% | 71.6 / 21.9 | 58 | G219G (GGT→GGG) | nanM | N‑acetylneuraminic acid mutarotase |
| Reads supporting (aligned to +/- strand): ref base A (23/6); new base C (0/28); total (23/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATT > NC_000913/4537962‑4538256 | ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCAcc < 1:28048/150‑1 (MQ=255)ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCt < 1:451244/151‑1 (MQ=255) cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt < 2:764885/150‑1 (MQ=255) aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc > 1:769975/1‑151 (MQ=255) cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAAc > 2:394909/1‑151 (MQ=255) ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAAccg > 2:100370/1‑151 (MQ=255) ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAAccg < 2:756527/151‑1 (MQ=255) aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAAccgcc > 1:24070/1‑151 (MQ=255) tACGCCATCTGGAGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGCGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACcgccgc > 2:133758/1‑150 (MQ=255) ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTAAAATTTTTCCCGGTGAAATCAATTTAAAATCCGGCATCCGTCCGCATTCCTGGTTTGG‑C‑TTCCCCATTAATAAGCCAGGTTTTATCCCCTTTATTC‑ACAACCGCCGCACCa < 2:676880/150‑1 (MQ=255) caaTCTGTGGAGGAGACGGGAGCAAGTTTTTTCCATTTTAATTTAATCCGGGTAAAATCAAGTTAAAATCCGCCTTCCTTTCGCAATCCTGTTTTGG‑C‑TTCCCCATAAAAAAGCCAGGTTTAATCCCCTTTATTC‑ACAACCGCCGCACCa < 2:805422/148‑1 (MQ=255) tGGTGATGAGACGGGAGCAAGCTTACTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGtt > 2:330103/1‑151 (MQ=255) ggTGATGAGACGGAAGCAACCTTTTCCCATTTTAATTTATTCCCGGTAAAATAAATTTAAAATACGGCATCCTTCCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGTTTTTATCCCTTTTATTC‑ACAACCGCCGCACCAGCCg < 2:849763/148‑1 (MQ=255) gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTcc > 2:550432/1‑151 (MQ=255) gaCGGGAGCAAGTTTATTCCATTTAAAATTATTCCCGGTAAAATCAAGTCAAAATCCGCCATCCGTCCGCAACCCTGTTTTGG‑T‑TTCGCCATAAAAAAGCCAGGTTTTATCACCTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCa < 1:210106/150‑1 (MQ=255) gggAGCAAGTTTATCCCATTTAAATTTATTCCCGGTAAAATAAATTTCAAATACGGCACCCTTTCGCAATCCTGTTTTGC‑T‑TTCGCCATAAAAAAGCCAGTTTTTATCACCTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGc < 1:868940/151‑1 (MQ=255) gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGc > 1:79588/1‑151 (MQ=255) gggAGAAAGCTTATTCCATTTAAATTTATTCCCGGTAAAATAAATTTAAAATCCGCCCCCCGTCCGCAATCCGGGTTTGG‑C‑TCCCCCATTAATAACCCGGGTTTTACCCCTTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAggg < 2:25796/150‑1 (MQ=255) ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt > 2:505137/1‑151 (MQ=255) tttaAAATTATTCCCGGGAAAATAAATTTAAAATACGCCCCCCGTTCGCAATCTGGTTTTGG‑T‑TTCCCCATTAATAAGCCAGTTTTTACCCCTTTTATCC‑AAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc < 2:450450/149‑1 (MQ=255) tattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg > 2:726609/1‑151 (MQ=255) ccggaaaaataaatTTAAAATACGCCATCCTTTCGCACCCCTGTTTTGG‑T‑TTCGCAATTAATAAGCGAGTTTTTATCCCTTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:270944/140‑1 (MQ=255) gtgaaaaCCAATTTAAAATCCGGC‑CCCCTTCCCTATCCTTGGTTGG‑CTTTCCCCATAAAAAAGCAAGTTTTTCTCCCTTTTTTCC‑CCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg < 2:584989/145‑1 (MQ=255) ggTAAAAACAAGTTCAAATACGGCATCCTTCCGCACCCCTGTTTTGG‑C‑TCCGCCTTAAATAAGCCGTTTTTTATCCCCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt < 2:875707/147‑1 (MQ=255) aaaTCAAGTCAAAATCCGGACTCGTTCCGCATTCCTGGTTGGC‑T‑TTCCCCATTAATAGGCCGGTTTTTATC‑CCTTTTTTCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt < 2:460826/151‑1 (MQ=255) aaaCCAATTTAAAATACGGCATCCTTCCCCAACCCGGTTTTGG‑C‑TCCCCCATTAAAAACCCGGTTTTTACCCCTTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt < 2:103037/151‑1 (MQ=255) aaGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa > 1:294896/1‑150 (MQ=255) ttCAAATACGGCATCCGTCCGAAATCCTGGTTTGG‑T‑TTCGCCATTAATAAGCCAGGTTTTATCCCCTTTATTC‑AAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg < 1:423760/149‑1 (MQ=255) aaaTAGGCCATCCTTCCGCAATCCTGTTTTGG‑C‑TCCCCCATTAATAAGCCAGTTTTTACCCCTTTTTTCC‑AAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa < 2:355792/150‑1 (MQ=255) aaaTACGGCACCCGTTCGCAATCCTGTTTTGG‑C‑TCCGCCATTAAAAAGCCAGGTTTTATCCCCTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa < 1:283955/150‑1 (MQ=255) aaTACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCCCCATAAAAAAGCCAGTTTTTATCCCCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca < 1:154284/151‑1 (MQ=255) gccATCCGTTCCCAACCCTGGTTTGG‑C‑TTCGCCATAAAAAACCAAGGTTTTACCCCTTTTATTC‑ACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa < 1:778070/147‑1 (MQ=255) cGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg > 2:712298/1‑151 (MQ=255) ttCGCAACCCTGGTTTGC‑C‑TTCGCCATAAATAAGCCGGTTTTTATCCCTTTTATTC‑ACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga < 2:379392/150‑1 (MQ=255) ctcacccTTGTTTGGCC‑TCCCCCATCAATAACCACGGTTTCTTCCCGTCTATCC‑ACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgca < 2:159043/144‑3 (MQ=255) aacccTGTTTTGG‑C‑TTCCCCATTAATAAGCCAGTTTTTACCCCCTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat < 1:565555/148‑1 (MQ=255) aaTCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat > 1:835465/1‑151 (MQ=255) ggTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt > 1:472511/1‑151 (MQ=255) ggTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt > 2:891121/1‑151 (MQ=255) ttGG‑C‑TTCGCCATTAAAAAGCCGGTTTTTATCCCCTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct < 1:870678/150‑1 (MQ=255) ttggttcccccATAAAAAACCCGGTTTTTACCCCTTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct < 2:95922/141‑1 (MQ=255) ttCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt > 1:825111/1‑151 (MQ=255) tccGCCATTAATAAGCCGGTTTTTATCCCTTTTTTCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg < 2:812859/144‑1 (MQ=255) gCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGTCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg > 2:26750/1‑151 (MQ=255) gCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg > 2:288237/1‑151 (MQ=255) ccATTAAAAAGCCAGTTTTTACCCCCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt < 1:803051/147‑1 (MQ=255) ccATAAATAAGCCAGTTTTTATCCCCTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt < 1:880039/151‑1 (MQ=255) ataaaGCCAGGTTTAACCCCCTTTATCC‑ACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa < 1:161706/147‑1 (MQ=255) aataaGCCAGTTTTTACCCCTTTTATCC‑CCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa < 2:718741/149‑1 (MQ=255) taaGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt > 2:621154/1‑150 (MQ=255) agacccgggTTTTACCCCTTTAACTC‑ACCCCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAACTTCTGCTTTTTTGTCaaaa < 2:65140/142‑1 (MQ=255) aaGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt > 2:771536/1‑149 (MQ=255) aaGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt > 1:624070/1‑149 (MQ=255) aGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt > 2:64761/1‑151 (MQ=255) ccggtttttATCCCCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAGGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa < 1:508334/146‑1 (MQ=255) caagtttttATCCCCTTTTTCC‑ACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa < 1:1566/146‑1 (MQ=255) gtttttATCCCCTTTATCC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAg < 1:746477/147‑1 (MQ=255) aTCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtt > 1:99414/1‑151 (MQ=255) | CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTC‑ACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATT > NC_000913/4537962‑4538256 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |