| 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.9% | 94.1 / 16.9 | 57 | G219G (GGT→GGG) | nanM | N‑acetylneuraminic acid mutarotase |
| Reads supporting (aligned to +/- strand): ref base A (28/4); new base C (0/24); total (28/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.40e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT > NC_000913/4537960‑4538258 | aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCAc > 1:7556/1‑151 (MQ=255) cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCGCCtt < 1:797630/151‑1 (MQ=255) ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCttt > 1:495543/1‑151 (MQ=255) ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCttt < 1:19483/151‑1 (MQ=255) cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTAtt > 1:390278/1‑151 (MQ=255) ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCTTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACaa < 2:587227/151‑1 (MQ=255) cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAAc > 1:572145/1‑151 (MQ=255) gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAAcc < 2:313965/151‑1 (MQ=255) ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAAccg > 1:707824/1‑151 (MQ=255) cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAAccgc > 2:3771/1‑151 (MQ=255) cTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACcgccgc > 1:219517/1‑151 (MQ=255) aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCAc > 2:395921/1‑151 (MQ=255) gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCAcc > 2:664561/1‑150 (MQ=255) gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCAcc > 2:33553/1‑150 (MQ=255) aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGcc > 1:109949/1‑151 (MQ=255) tCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTAAAATTATTCCCGGTGAAATCAAGTTAAAATCCGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATAAATAGGCCG‑GTTTTTATCCCCTTTATT‑CACAACCGCCGCACCa < 2:123436/147‑1 (MQ=255) ggTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTc > 1:115680/1‑151 (MQ=255) cGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAggg > 2:767012/1‑151 (MQ=255) gggAGAAAGCTTATCCCTTTTAAAATTATTCCCGGTAAAATAAATTTAAAATCCGG‑CA‑TCCGTTCGCATCCCTGT‑TTTGGTTTCCCCATTAATAAGCCA‑GGTTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCa < 1:387753/147‑1 (MQ=255) gattATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCCGGGCGATTCGcc > 2:666948/3‑151 (MQ=255) ttATTCCATTTAAATTTATTCCCGGTAAAATAAATTTAAAATCCGG‑CA‑CCCGTTCGAATCCCGGG‑TTTGGTTTCCCCATAAAAAAGCAA‑GTTTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc < 1:463941/149‑1 (MQ=255) aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg > 1:698455/1‑151 (MQ=255) ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt > 1:400862/1‑151 (MQ=255) tCCATTTAAAATT‑CTACCCGGAAAACGAACTTAAAAAACGG‑CCTTCCGGTCGCAAGCCGGT‑TTGGGCTTCCCCATTAA‑AAGGCC‑GGTTTTCACCCCTTTTTTCCAACCCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg < 2:175010/149‑1 (MQ=255) tttaAATTTATTCCGGGTAAAATCAATTTAAAATCCGG‑CC‑TCCTTGCACAACCCTGT‑TTTGGTTGCGCCATAAAAAAGCCG‑GGTTTTATCCCTTTTTTA‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc < 2:41807/149‑1 (MQ=255) tttaAAATTATTCCCGGTAAAATCAAGTTAAAATCCGC‑CA‑CCCGTTCGCAATCCTGG‑TTTGCTTCCCCCATAAATAACCCA‑GTTTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa < 1:12857/147‑1 (MQ=255) tttaAAATTATTACCGGGAAAATAAATTAAAAATACGC‑CA‑TCCGTTCCCAATCCT‑G‑TTTTGCTTCCCCCATAAAAAGGCC‑GGTTTTTACCCCCTTATTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa < 2:728829/147‑1 (MQ=255) aaaTTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt > 1:51968/1‑151 (MQ=255) aaTTATTCCGGGTAAAATAAATTTAAAATACGG‑CA‑CCCGTCCCAAACCTTGT‑TTTGGTTTCGCAATAAATAACCCA‑GTTTTTATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg < 1:763481/151‑1 (MQ=255) ttaCGGGTAAAAT‑AAATTCAAAATCGGCCA‑TCCTTCCGCAATCCTGG‑TTGGGCTCCCCCATAAAAAAGCCG‑GTTTTTCTCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 1:5191/149‑1 (MQ=255) taCCGGTGGAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt > 2:722182/1‑151 (MQ=255) aaaTCAA‑TTTAAATACCG‑CC‑CCCCTTCCGCATCCTGGTTTTGGCTCCGCCTTTAAAAACCCG‑GTTTTTACCCCTTTTATC‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa < 2:718695/150‑1 (MQ=255) aagtaaaaaTCCGC‑CA‑CCCGTTCGAATCCCTGT‑TTTGCTTTCCCCATTAAAAACCCA‑GGTTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaaa < 1:192810/145‑1 (MQ=255) aaGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa > 2:632230/1‑150 (MQ=255) aaaTCCGG‑CA‑CCCTTCCGCATTCCGGT‑TTTGCCTGCGCCA‑TAAAAAGCCAGGGTTTTACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa < 1:565100/147‑1 (MQ=255) aaaTACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac > 1:705023/1‑151 (MQ=255) aaTACGG‑CA‑CCCTTTCGCAATCCTGT‑TTTGGCTTCGCCATTAATAAGCCA‑GTTTTTATCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca < 1:258279/151‑1 (MQ=255) tCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt > 1:762152/1‑151 (MQ=255) cgttccgaatcccTGT‑TGTGGCTTCCCCA‑TAAAAACCCC‑GGTTTTTACCCCTTTTTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa < 2:87588/133‑1 (MQ=255) cgaattCCTGT‑TTTGCCTCCGCCATTAATAAGCCA‑GTTTTTATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga < 1:3771/146‑1 (MQ=255) aaTCCTGG‑TTTGGTTTCGCCATAAATAGCCCG‑GTTTTTACCCCTTTTATT‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACGGAAACTTATTgaag < 1:39110/147‑1 (MQ=255) accccggttttGCCTTCCCCATAAATAAGCAG‑GTTTTTACCCCCTTTCTC‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga < 2:777212/137‑1 (MQ=255) ctggg‑TTTGGCTCCCCAATTAAAAACCCA‑GTTTTTACCCCTTTTATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata < 1:481664/146‑1 (MQ=255) ccTGT‑TTTGGTTCCGCCATTAAAAAGCCA‑GTTTTTATCCCTTTTATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat < 1:624837/148‑1 (MQ=255) ccTGG‑TTTGGTTTCCCCATTAAAAACCCA‑GTTTTTACCCCCTTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa < 1:809051/150‑1 (MQ=255) ccTGG‑TTGGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACAGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa < 1:422448/150‑1 (MQ=255) tGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct > 2:319389/1‑151 (MQ=255) tttGGCTCCCCCATTAAAAAGCCA‑GTTTTTATCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc < 1:242149/150‑1 (MQ=255) tGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc > 2:626236/1‑151 (MQ=255) ggCTTCGTCATTAATAAGCCA‑GGTTTTATCAACTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCAAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt > 2:146947/1‑151 (MQ=255) cGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt > 1:773958/1‑151 (MQ=255) cataaaaacgccggtttttACCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt < 2:841497/135‑1 (MQ=255) ccA‑GTTTTTTCCCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg < 2:377203/149‑1 (MQ=255) ccA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGAGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt > 2:710685/1‑149 (MQ=255) ccggtttttATCCCCTTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa < 1:668801/146‑1 (MQ=255) ttttATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt > 2:842681/1‑151 (MQ=255) cACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCCGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt > 2:347900/1‑151 (MQ=255) | ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGG‑CA‑TCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCA‑GGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT > NC_000913/4537960‑4538258 |
| 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: ‑ |