| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 298,960 | 0 | A | C | 44.4% | 79.1 / 17.8 | 54 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base A (25/5); new base C (0/24); total (25/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.09e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGC > NC_000913/298814‑299110 | ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcccc > 2:6503/1‑151 (MQ=255) acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCg < 1:869806/150‑1 (MQ=255) ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCa < 1:852611/150‑1 (MQ=255) gacaGCCGCGCTCACGCCGCATAGCCCCAGCTCACCGACACTTTTGG‑CCTTGCTCGGGGAGGATATGGGGTCGGT‑ATCATACGGGATAATTACCTCCGGTTTTGGCATATCCGCATGATCCAGCATCTCATACCCCGCCGTATCg > 2:743934/4‑145 (MQ=255) ccGCGCTTGCGCCGCACAGACCCAGTCAACAAACACGTTTGC‑CCTTCATGGGGGGGGAAATGGGGGGAAA‑ACCCCCCAGCAAAATCCCCTCCTTTTTGTGAAAATCCGCATAATCCGACGCTTGTAACTCCGCCATATCGTGATTAACGaa < 2:1007896/151‑1 (MQ=255) ccGCGCTCACGCCGCACAGGCCCAGCTCCCCGCCCCCTTTGG‑CCTTCACCGGGGAGGATAGGGGGTCGTT‑ATCACCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAAc < 1:678523/148‑1 (MQ=255) ccggcTCGCGCCGCACGGGCCCAGCTCACCCACCCTTTTGG‑CCTTCATCGGGGGGGATATGGGGTCGTT‑ATCACCCAGGAAAACCCCCCCCTTTTTGGGGATACCCGCATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAg < 2:682961/148‑1 (MQ=255) cACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc > 2:726574/1‑151 (MQ=255) cgccgcACAGGCCCAGCTCACCGACCCCTTTGC‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa < 1:969979/150‑1 (MQ=255) ccgcACGGGCCCAGCTCCCCGAC‑CCTTTGGCCCTTCACCGGGGAGGATA‑GGGGTTGGTAATCACCAAGAAAAACCCCCCCCTTTTTTGGAATACCCACATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc < 2:414502/151‑1 (MQ=255) ccgcACAGGCCCAGCTCACCGACCCCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAACCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc < 1:241351/151‑1 (MQ=255) acaGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc > 1:1047176/1‑151 (MQ=255) cAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAcc > 1:722293/1‑151 (MQ=255) ccccccGACCCCTTTGG‑CCTTCAGCGGGGGGGATAGGGGGTCGGT‑ATCACCCAGGAAAATCCCCTCCTTTTGTGGTATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc < 2:200994/146‑1 (MQ=255) tCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc > 1:172266/1‑150 (MQ=255) cACCGACCCCTTTGG‑CCTTCATGGGGGAGAATATGGGGTCGTT‑ATCACCCAGAAAAACCCCCTCCTGTTTGGGGATACCCGCATGACCCGCCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAg < 1:120226/151‑1 (MQ=255) ccccTTTGC‑CCTTCAGGGGGGAGGATATGGGGTCGTT‑ATCATCCAGAAAAACCCCCTCCTTTTTTGGGATCCCCGCATGACCCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcct < 2:256963/148‑1 (MQ=255) ccTTGGG‑CCTTCAGGGGGGAGGATAGGGGGTCGGT‑ATCACCCAGAAAAACCCCCTCCTTTTTGGGGATACCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcct < 1:628710/148‑1 (MQ=255) tGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc > 1:689683/1‑149 (MQ=255) tCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc > 1:482218/1‑151 (MQ=255) agcGGGGAGGATATGGGGTCGTT‑ATCACCCAGAAAAACCCCCCCCTTTTTTGGGATCCCCGCAGGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCgcgc < 1:947199/149‑1 (MQ=255) gaggaTAGGGGGTCGGT‑ATCACCCAGAAAAACCCCCCCCTTTTTTGGGATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgc < 1:1057367/151‑1 (MQ=255) tatGGGGTCGGT‑ATCACCCAGAAAAACCCCCTCCTTTTTTGGGGTACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat < 1:1058009/150‑1 (MQ=255) atGGGGTCGGT‑ATCATCTAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat > 2:647701/1‑149 (MQ=255) tGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt > 1:11021/1‑151 (MQ=255) ggTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAt > 1:936479/1‑151 (MQ=255) ggTCGGT‑ACCACCCAGAAAAACCCCCCCCTTTTTTGGAATACCCGCAGAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAt < 2:149736/151‑1 (MQ=255) gTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAtt > 2:963425/1‑151 (MQ=255) gT‑ACCACCCAGAAAAACCCCCTCCTGTTTGGGGATACCCGCGTGACCCGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg < 2:519300/148‑1 (MQ=255) accgagaaaaaCCCCCCCCTTTTTGGGTATACCCACAAGACCCGGCCCCTCTTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCa < 1:376076/141‑1 (MQ=255) acccAGAAAAACCCCCCCCTTTTTGGGAATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCa < 1:843292/146‑1 (MQ=255) ccAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc > 1:24425/1‑151 (MQ=255) ccAGAAAAACCCCCTCTTGTTTTGGGATACCCGAATGACCCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc < 2:4784/151‑1 (MQ=255) gAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTgg > 2:742275/1‑151 (MQ=255) aaaaTCCCCTCCTGTTTTGGGATATCCGCATAAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg < 1:347501/149‑1 (MQ=255) aaaaCCCCCTCCTTTTTTGGAATACCCGAATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg < 1:540934/149‑1 (MQ=255) aaaTCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt > 1:679693/1‑151 (MQ=255) aagcccctcttTTTTTGGGAATCCCGCAGGCCCCGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg < 2:631713/146‑1 (MQ=255) cACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc > 2:922838/1‑151 (MQ=255) ttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg > 2:373425/1‑151 (MQ=255) gggATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGAt > 1:820774/1‑151 (MQ=255) tatCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc > 2:886311/1‑149 (MQ=255) cGCAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg < 1:212432/151‑1 (MQ=255) gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc > 1:385053/1‑151 (MQ=255) tGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc > 1:984552/1‑151 (MQ=255) acccgcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc < 2:17199/141‑1 (MQ=255) acccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc < 1:314237/148‑1 (MQ=255) aaCCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt > 1:1023196/1‑151 (MQ=255) ccGGCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg < 1:632662/150‑1 (MQ=255) ccGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc > 2:135773/1‑151 (MQ=255) cACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgca > 1:723344/1‑148 (MQ=255) ccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcaca < 2:698264/151‑1 (MQ=255) ccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcaca < 1:51208/151‑1 (MQ=255) tCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac > 2:460016/1‑150 (MQ=255) aCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGc > 2:614582/1‑151 (MQ=255) | ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGT‑ATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGC > NC_000913/298814‑299110 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |