| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,538,109 | 0 | A | C | 46.5% | 53.0 / 16.4 | 45 | G219G (GGT→GGG) | nanM | N‑acetylneuraminic acid mutarotase |
| Reads supporting (aligned to +/- strand): ref base A (19/4); new base C (0/20); total (19/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.65e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG > NC_000913/4537966‑4538253 | gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑GTTTATCACCTTTAt < 2:281111/151‑1 (MQ=255) ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑caca > 2:561308/1‑150 (MQ=255) tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACcgccgc > 2:19988/1‑150 (MQ=255) aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCAc > 2:476105/1‑151 (MQ=255) tgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTAc > 2:259963/1‑151 (MQ=255) gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTAcc > 1:563502/1‑151 (MQ=255) gggAGCAAGTTTTTCCCTTTTAAAATTATCCCCGGTAAAAC‑AAATTAAAAATACGCCA‑CCCTTCCCAACTCTTGTTTTGC‑CTCCCCCATAAAAAGCCCAGT‑TTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAggg < 2:55223/150‑1 (MQ=255) gggAGCAAGTTTATTCCTTTTAAAATTATTCCCGGTAAAAT‑AAAGTTAAAATCCGGCATCCGTTCGCAA‑CCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGT‑TTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAg < 1:119450/148‑1 (MQ=255) gggAGCAAGTTTATTCCTTTTAAAATTATTACCGGTGAAAT‑AAAGTTAAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATAAAAAACCCAGT‑TTTTATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAgg < 1:586003/149‑1 (MQ=255) gggAGCAACTTTATTCCTTTTAAATTTTTTCCCGTTAAAAT‑CAAGTTAAAATCCGGCATCCTTCCGC‑A‑TCCCTGTTTGG‑GCTCCCCAATAAATAGCCAGGTTTTTATC‑CCTTTATTCCACACCCGCCGCACCAGCCGTTCCGTACCAggg < 2:624051/150‑1 (MQ=255) ggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg > 2:315543/1‑151 (MQ=255) ggAGCAA‑CTTTTTTCATTTTAAATTTTTACCGGGAAAATCCAAGTTAAAAGCGGGCACCCGTTCCCAA‑‑CCTTGGTTGG‑GTTCCCCCTTAAAAAACCCGGTTTTTACCCCTTTTTTT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg < 2:7765/151‑1 (MQ=255) gAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGc > 2:488147/1‑149 (MQ=255) aGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt > 1:433282/1‑151 (MQ=255) tCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt > 2:15136/1‑150 (MQ=255) tttAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc > 2:651618/1‑151 (MQ=255) ccGGTAAAAT‑CAAGTTAAAATACGGCATCCGTTCGCAA‑TCCGGTTTTGG‑CTTCCCCATTAAAAAGCCAGG‑TTTTATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt < 1:109589/150‑1 (MQ=255) tGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg > 1:426669/1‑151 (MQ=255) gttaaaaaaaattCAAAACCCGCCTCCCTTCCAAA‑ACCTGGTTTGGCCTGCGCCATAAAAAAGCCGGG‑TGTTATCCCTTTAATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt < 2:339532/137‑1 (MQ=255) aaaT‑CAAGTTAAAATCCGGCACCCTTTCGCAA‑TCCTGTTTTGG‑TTTCGCCATAAAAAACCCAGT‑TTTTACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt < 1:539650/151‑1 (MQ=255) gtacAATTACGGCATCCGTCCGAAT‑TCGTGGTTGGG‑CTCCGCCATTAAAAAGCCGGT‑TTTTCCCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg < 1:212667/147‑1 (MQ=255) ccgcatccctttCGCAT‑CCTTGTTTTGC‑CTCCCCCATAAATAACCCGGG‑TTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa < 2:528100/141‑1 (MQ=255) actgtaTCCCTTCGAAA‑TCTTGTTTTGG‑TTGCCCCATAAAAACGCCAGT‑TTTTACCCCTTTTATC‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg < 1:573871/144‑1 (MQ=255) ccGTCCGCAA‑TCCTGTTGTGG‑CTCCGCCATAAAAAAGCCGGG‑TTTTCCCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt < 2:543384/150‑1 (MQ=255) gctCGCAC‑TCC‑GGTTGGGCCTTCCCCATTAAAAAGCCAGT‑TTTAATCCCTTTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga < 2:100224/149‑1 (MQ=255) cAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa > 1:348421/1‑150 (MQ=255) aa‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa > 2:248544/1‑150 (MQ=255) aacccggttttGG‑CTCCCCCATAAAAAAGCCGGT‑TTTTATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga < 1:709403/140‑1 (MQ=255) cTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc > 1:140674/1‑151 (MQ=255) cTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc > 2:530934/1‑151 (MQ=255) ttttGC‑CTTCGCCATAAAAAAGCCGGT‑TTTTACCACCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt < 1:343076/149‑1 (MQ=255) gg‑CTTCCCCATTAATAGGCCGGT‑TTTTATCCCTTTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg < 1:316061/149‑1 (MQ=255) cGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt > 1:491488/1‑151 (MQ=255) cGCCATAAAAACGCCAGT‑TTTAACCCCCTTTTTC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt < 2:428070/151‑1 (MQ=255) aaaaaGCCGGT‑TTTTATCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa < 1:299644/146‑1 (MQ=255) aaaaggccggtttttAGCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa < 2:210635/139‑1 (MQ=255) aaaaGCCAGT‑TTTTATCACTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt < 2:499143/149‑1 (MQ=255) aaccccAGT‑TTTTCCCCCCTTTATC‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt < 2:70576/146‑1 (MQ=255) aaGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt > 1:574799/1‑149 (MQ=255) ccAGT‑TTTTATCCCCTTAATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa < 2:246137/147‑1 (MQ=255) aGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc > 2:367066/1‑151 (MQ=255) aGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc > 1:397627/1‑151 (MQ=255) gtttttATCTCTTTTATT‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc < 1:669663/148‑1 (MQ=255) tttATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt < 1:215474/150‑1 (MQ=255) ttctCACTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg < 2:55831/147‑1 (MQ=255) | GCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAAT‑CAAGTTCAAATACGGCATCCGTTCGCAA‑TCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG > NC_000913/4537966‑4538253 |
| 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: ‑ |