| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,001,164 | 0 | A | C | 51.5% | 46.1 / 56.9 | 33 | T274P (ACC→CCC) | xdhA | xanthine dehydrogenase, molybdenum binding subunit |
| Reads supporting (aligned to +/- strand): ref base A (15/1); new base C (0/17); total (15/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.54e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-08 | |||||||||||
GCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCG > NC_000913/3001018‑3001314 | gCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCAGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑TCTCGCAACCCg > 1:715854/1‑151 (MQ=255) tACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTAc < 2:501735/151‑1 (MQ=255) ttttGTGGGTGGCGTTTTTGGAAAAAAACAGGATGTACGGAAAGAGCAAAGGGCGGCATTCCTGACCAGCAAGCTGGGCGG‑CATCCCGGTA‑AAATTTCCCCTTAGCCGT‑AAAGGGTGTTC‑CCCCGCACCCCGTCCCCGCCACGCTTTTAcc < 2:832221/150‑1 (MQ=255) ggtgtttttGATAAAAAACGGGAGGTAGGGAAAGAGCCAAGGGCGGCCTTCCTGCCCAGCAAGGTGGGCGG‑CATCCCGGTA‑AAATTTCCCCTTAGCGGT‑GAAGAGTTTTC‑CCTCCCACCCCGTCCCCGCCACGCTTTTACCATTGACggg < 1:492941/145‑1 (MQ=255) gTAATAACCAGGTTGTACGGGAAGGCCCAAGGGCGGCATTCCTGCCCAGAAAGCTTGGCGG‑CTTTCGGGTA‑AAAGTTCCCCTGAGCCGA‑GAAGGTTTTTT‑CCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCg < 2:315697/150‑1 (MQ=255) aaaCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAc > 1:233993/1‑150 (MQ=255) aaaCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAc > 2:148055/1‑150 (MQ=255) gTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGaa > 2:241539/1‑151 (MQ=255) gTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGaa > 2:200384/1‑151 (MQ=255) gTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGaa > 2:385516/1‑151 (MQ=255) agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAgg > 1:313294/1‑151 (MQ=255) ggcggcATGCCAGCCCAGCAAGCTGGGCGG‑CTTCCCGGTA‑AAAGTTCCCCTTGCCGGA‑GAAGAGTGTTC‑CCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAg < 1:669349/150‑1 (MQ=255) ggccATTCCTGACCAGAAGGTTTGGCGC‑CTTTCCGGGA‑AAAGTTCCCCTTAGCGGT‑GAAGAGTGTTC‑CCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTc < 1:399234/147‑1 (MQ=255) ttCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg > 2:341255/1‑151 (MQ=255) ttCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg > 2:249090/1‑151 (MQ=255) cTGACCAGCAGTCTTGG‑GG‑CCTTCCCGTT‑AAAATTTCCCCTTACCGTGGAAGAGTTTTT‑CCCGACACCCCTCACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc < 2:344965/151‑1 (MQ=255) cTGACCA‑CAAGGTTGGGGGCCTTTCGGGTA‑AAATTTCCCCTTGCCCGT‑AAAGAGTTTTC‑CCCCGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc < 2:557328/151‑1 (MQ=255) aGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa > 2:771092/1‑151 (MQ=255) gCTGGGCGG‑CCTCCGGGTA‑AAATTTCCCCTTAGCGGT‑GAAGGTTGTTC‑CCTCGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAacac < 1:608715/150‑1 (MQ=255) ctgggcgg‑cATCCGGGTA‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCCCACCCCGTACCCGCCACGCTTTTACCATTGACGAGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAacac < 2:289802/146‑1 (MQ=255) ggtggggggcttCCGGGGAAAAGTCCCCCTTGGCCGG‑GAAGATTGTTC‑CCCCGCACCCCGGCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAacac < 1:523643/140‑1 (MQ=255) cgg‑cATTCGGGTA‑AAAGTTCCCCTTAGCCGA‑GAAGAGTTTCC‑CCCCCCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgcgc < 2:577274/150‑1 (MQ=255) cgggTA‑AAATTTCCCCTTACCGGA‑GAAGAGTTTTC‑CCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGGACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAt < 2:293226/149‑1 (MQ=255) gTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACggg > 1:80770/1‑150 (MQ=255) gTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACggg > 2:651465/1‑150 (MQ=255) cccTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACt > 1:866348/1‑150 (MQ=255) cccTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACt > 1:206461/1‑150 (MQ=255) ccTTGGCCGT‑GAAGGGTTTTC‑CCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTcc < 1:63180/151‑1 (MQ=255) tgagccggaaggagTTTTC‑CCCCCCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCAt < 2:668615/140‑1 (MQ=255) attgtTC‑CCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg < 1:316859/148‑1 (MQ=255) gtgtTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg < 1:423393/150‑1 (MQ=255) ggttttcccTCGCACCCCCTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg < 1:808892/143‑1 (MQ=255) cTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATa > 1:753232/1‑151 (MQ=255) aCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCg > 1:406775/1‑151 (MQ=255) aCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCg > 1:39649/1‑151 (MQ=255) | GCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGG‑CATTCCGGTG‑AAAGTTTCCCTTAGCCGT‑GAAGAGTGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCG > NC_000913/3001018‑3001314 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |