| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,538,109 | 0 | A | C | 42.6% | 84.0 / 16.6 | 54 | G219G (GGT→GGG) | nanM | N‑acetylneuraminic acid mutarotase |
| Reads supporting (aligned to +/- strand): ref base A (26/5); new base C (0/23); total (26/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.29e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTT > NC_000913/4537961‑4538257 | tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCa < 2:80317/149‑1 (MQ=255) cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑ttt > 1:279836/1‑150 (MQ=255) gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTAt < 2:135603/151‑1 (MQ=255) cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTAtt > 2:228531/1‑151 (MQ=255) aaaaCCGCCAGGTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAATTTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATACGCCAGGTTTTATCACC‑TTTATT‑c > 2:215467/1‑151 (MQ=255) aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑ca‑c > 1:672474/1‑150 (MQ=255) aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑ca‑c > 1:624463/1‑150 (MQ=255) gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAAcc > 1:308539/1‑151 (MQ=255) gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTCATT‑CA‑CAAcc < 2:468020/151‑1 (MQ=255) ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAAccg < 2:660350/151‑1 (MQ=255) gCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAAccgccg < 2:357898/151‑1 (MQ=255) cTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACcgccgc > 1:299445/1‑151 (MQ=255) ggTGAGGAGGAGGGAGCAAGTTTATTCTTTTTTAATTTATTACGGGTAAAATAAAGTTAAAATA‑CGGCTTCCGTTCTCATCCTTGTTTTGGCTT‑CGCCATTAATAAGCCGGGTTTTACCCCT‑TTAATC‑CC‑CAACCGCCGCACCAGcc < 2:898659/147‑1 (MQ=255) tgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTAc > 1:178023/1‑151 (MQ=255) gaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTAcc > 1:63579/1‑149 (MQ=255) aCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAgg > 1:605235/1‑151 (MQ=255) gggAGCAAGCTTTTTCCATTTTAATTTATTACCGGTAAAATAAAGTTAAAATC‑CGGCATCCGTTCGAAATCCTGTTTTGGTTT‑CGCCATTAATAAGCCAGGTTTTATCCCC‑TTTATC‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGc < 1:810568/151‑1 (MQ=255) gggAGAAAGCTTTTTCCATTTAAAATTATTCCCGGGAAAATAAAGTAAAAATC‑CGCCTCCCGTTCCAAACCCTGGTTTGGCTT‑CGCAATTAAAAACCCGGGTTTTACCCCC‑TTTATC‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGc < 2:1009816/151‑1 (MQ=255) aatctttttcttttttAATTTTTTACCGGTAAAAACAA‑TTCAAACACCGGCACCCGTCCGCACTCCTGTTTTGGTTT‑CGCCATAAATAAGCCGGGTTTTACCCCT‑TTTATT‑CA‑AACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc < 2:29729/138‑1 (MQ=255) ttCCTTTTAAAATTATTCCCGGGAAAATCAATTTAAAATC‑CGGCTTCCTTTCGCAACCCTGTTTTGGCTT‑CCCAATAAATAACCCGGGTTTTATCCCT‑TTTATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg < 1:506161/150‑1 (MQ=255) tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt > 1:62553/1‑150 (MQ=255) cctttttAAATTTTTCCGGGGAAAAAAAAGTTAAATTC‑CGGCACCCTTTCGAAATCTTGGTTTGGTTT‑CCCAATAAAAAAGCATGTTTTTACCCCT‑TTTATC‑CA‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa < 1:1048694/148‑1 (MQ=255) ttattcCCGGTAAAATAAAGTAAAAATA‑GGGCATCCGTCCGAAATCCTGGTTTGGCTT‑CGCCATAAATACGCAGGTTTTTACCCCT‑TTTATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 2:450416/150‑1 (MQ=255) tattaCCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg > 1:99628/1‑151 (MQ=255) aCCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt > 1:337201/1‑151 (MQ=255) ccGGTAAAATCAATTTAAAATA‑CGGCACCTGTTCGCAATCGTGGTTTGGCTC‑CCCAATAAAAAGGCCGGTTTTTATCCCT‑TTTATA‑CC‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg < 1:228556/148‑1 (MQ=255) aatcaaatttaaaaTA‑CGGCTTCTTTCCGCATTCCTGTTTTGGCTT‑CCCCATTAATAGGCCAGTTTTTATCCCT‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa < 2:370226/139‑1 (MQ=255) aaaaGTTAAAACA‑CGGCGCCCTTCCACATTCCTGTTTTTGTTC‑CGCATTTAAAAAGCCGGGTTTTTCCCCT‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa < 2:292776/148‑1 (MQ=255) aaGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa > 1:463060/1‑150 (MQ=255) aaGTCCAAAAC‑CGCCCTCCTTTCCCAACCTTGTTTTGCCTC‑CCCCATAAAAAAGCAGGTTTTTATCCCT‑TTTTTC‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa < 2:1013389/149‑1 (MQ=255) aaaTT‑CGGCATCCGTCCGCAATCGTTTTGTTGCTC‑CCCCATTAAT‑AGCCCTGTTTTAACCCCTTTTATC‑CG‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac < 2:650677/148‑1 (MQ=255) gCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt > 2:626171/1‑151 (MQ=255) gCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt > 1:653219/1‑151 (MQ=255) cATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt > 1:323216/1‑151 (MQ=255) ctttccGCAAT‑CTGGTTTGGCTTCCCCCATAAATACGCC‑GGTTTTTCCCCC‑TTTTTTCCA‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt < 2:489972/142‑1 (MQ=255) ccttccGCAATCCTGTTTTGGCTC‑CGCCATTAAAAAGCCGGTTTTTATCCCT‑TTTATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt < 1:930106/142‑1 (MQ=255) ttCGCAATCTTGTTTTGGTTT‑CGCCATTAATAAGCCAGTTTTTATCCCC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga < 1:29790/150‑1 (MQ=255) tccccactCCTGTTTTGCCTC‑CCCTATAAATAACCCTGTTTTTTACCCT‑TTTTTC‑CA‑AACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg < 2:636991/142‑1 (MQ=255) gCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga > 2:247271/1‑150 (MQ=255) tCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata > 2:196015/1‑150 (MQ=255) ccTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCCCC‑TTTATC‑CA‑CAACCGCCGCACCAACCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata < 1:713037/149‑1 (MQ=255) tGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct > 2:477982/1‑151 (MQ=255) ttGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg > 1:541121/1‑151 (MQ=255) cgcttttAAAAAGCCAGTTTTTATCACC‑TTTATT‑CC‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt < 2:913760/146‑1 (MQ=255) taaatacGCCAGTTTTTATCCCC‑TTTTTC‑CA‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg < 2:588756/137‑1 (MQ=255) taataaGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa > 2:744236/1‑150 (MQ=255) aaaaaGCCAGGTTTTATCCCC‑TTTATA‑CA‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa < 1:1040123/146‑1 (MQ=255) aaaatcccggTTTTTTTCCC‑TTTTTT‑CCCCAACCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg < 2:787850/137‑1 (MQ=255) ccggtttttACCCCT‑TTTATC‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa < 1:12183/146‑1 (MQ=255) cAGGTTTTACCCCT‑TTTATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAg < 1:268485/151‑1 (MQ=255) aGGTTTTATCACC‑TTTATT‑CA‑GAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa > 2:622415/1‑149 (MQ=255) aGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc > 2:59215/1‑151 (MQ=255) ggTTTTATCCCT‑TTTTTC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa < 1:474372/144‑1 (MQ=255) ttATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg > 2:991815/1‑150 (MQ=255) tCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGCTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAACATAGTGAGCATTGAttt > 2:894449/1‑151 (MQ=255) | TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCATCCGTTCGCAATCCTGGTTTGGCTT‑CGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTT > NC_000913/4537961‑4538257 |
| 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: ‑ |