| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,346,675 | 0 | T | G | 54.4% | 42.4 / 68.8 | 48 | T218P (ACC→CCC) | fumB | anaerobic class I fumarate hydratase (fumarase B) |
| Reads supporting (aligned to +/- strand): ref base T (5/16); new base G (25/0); total (32/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.05e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-06 | |||||||||||
ACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGG > NC_000913/4346528‑4346820 | aCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACg < 1:933866/151‑1 (MQ=255) tGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGCGCCGTTCCCAGGGGACGCCttt > 2:607569/1‑151 (MQ=255) cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtctc < 2:639162/151‑1 (MQ=255) cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGAGCCAttttt > 1:317610/1‑148 (MQ=255) aTGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcg < 1:446022/151‑1 (MQ=255) ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgac < 2:341851/150‑1 (MQ=255) tcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgacg < 1:625742/150‑1 (MQ=255) cTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGCAGTACCGGGGGGACGCATTTTTTCGACGA‑GGAAggt > 1:983811/1‑149 (MQ=255) ttCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCCGTAACGCGGGTACGCCTTTTCTCGACGA‑GGAAGgtt > 1:1149962/1‑151 (MQ=255) ttCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGGACCGCGGGGAAGCATTTTTTCGACGA‑GGAAggt > 2:828213/1‑148 (MQ=255) ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGttttt < 1:313909/151‑1 (MQ=255) gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGTGCAGGACCGAGGGGACGCCATTTCTCGGCGA‑GGGAGgtttt > 2:204161/1‑149 (MQ=255) ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGtt < 2:355722/151‑1 (MQ=255) gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGttt < 2:414560/151‑1 (MQ=255) cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGGCGCATTTTTTTCACGG‑GGAAGGTTTTTAGttt > 2:635030/1‑150 (MQ=255) aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑c < 2:720389/151‑1 (MQ=255) aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAAGGCGCAGTACCGCGGGTACGCCTTTTTTCGCAGG‑GGAAGATTTTCAGTTTgt > 2:249638/1‑150 (MQ=255) cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCggg < 2:535348/151‑1 (MQ=255) cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGGACCGAGGGGACGGATTTTCTTGCCGAGGGAGGTTTTTCAGTTGG‑Cggggg > 2:817056/1‑149 (MQ=255) aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTAACGAGGGGACGCATTTTTTCGACGG‑GGGAGGTTTTTCGTTTT‑GCggg > 2:136925/1‑145 (MQ=255) aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTcagc‑a‑g < 2:1121330/151‑1 (MQ=255) tGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGAAGTTTTTTAGTTTG‑CCGGGAGGCAGC‑A‑ggg > 1:886786/1‑151 (MQ=255) aGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTACCGAGGGGACGCATTTTTTCGGCGA‑GGAAGGTTTTTAGTTTT‑TCGGGGGGCAGC‑C‑ggg > 1:200648/1‑149 (MQ=255) gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGA‑GGAAGTTTTTTAGTTTG‑CCGGGGGGCAGC‑A‑GGGCtt > 1:1019808/1‑151 (MQ=255) gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCCGGAACCAGGGGACGCATTTTTTCGACGG‑GGGAGGTTTTTAGTTTT‑GCGGGGGGCAGC‑C‑GGGCttt > 2:644868/1‑150 (MQ=255) tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGt < 2:387654/150‑1 (MQ=255) cTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGGGGCGGACCGGGGGGACGCGATTTTTTTACGA‑GGGAGGTTTTTAGGTTT‑TCGGGGGGCAGC‑AGGGGGTTTGGtt > 2:707696/1‑150 (MQ=255) ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCAGTACCGGGGGGACGCATTTTTTCGACGA‑GGAAGGTTTTCAGATTG‑GCGGGAGGCAGC‑A‑GGGCTTTGGttt > 1:866074/1‑145 (MQ=255) ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGGAGGTTTTCAGTTTG‑CCGGGGGTCAGC‑A‑GGGCGTTGGTTTCCTgg > 1:681209/1‑150 (MQ=255) gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCAGTACCGGGGGGAAGGATTTTTTTGACGA‑GGGAGGTTTTTAGGTTT‑CCGGGGGTCAGC‑A‑GGGGTTTGGgtttttgg > 2:370282/1‑143 (MQ=255) aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTaga < 1:930216/151‑1 (MQ=255) gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGTTAGGACCGCGGGGACGCCCTTTTTCCACGA‑GGGAGGTTTTTAGTTTT‑TCGGGAGGCAGC‑A‑GGGCTTTTGGTTTCTGGTTGGGGTAAGTTttgtt > 1:357546/1‑149 (MQ=255) cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGGACCCAGGGGACGCCTTTTTTCGACGG‑GGGAGGTTTTTTGTTTG‑CCGGGGGGGAGC‑A‑GGGCTTTGGTTTTCTGGGAGAGATACGTTTTTTTTGCa > 2:540246/1‑148 (MQ=255) aGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGGGCAGACTCCTGGGGAAGCATTTTCTCGTCGG‑GGGAGTTTTTGCGTTTT‑CCGGGAGTCAGA‑C‑GGGGGTTGGGTTTCTCGTATAGATAAGTTTTTTTGGcggcagc > 1:1047885/1‑144 (MQ=255) gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGGAGTTTTTTAGTTTG‑CCGGGGGGCAGC‑A‑GGGGTTTGGGTTCCTGGGAGAGATACGGTTTGTTGGGAGAGCCGCCGCCttt > 1:719149/1‑150 (MQ=255) ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCAAACCCAGGGAAAACATTTTTTCGTCCA‑AGAAGTTTTTTAGGTTT‑CCGGGGGGCCAG‑A‑GGGGGTTGGGTTCCGTGTAGAGAGAAATTTTTTTGGCggagggcccgccct > 1:1091091/1‑138 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCCAAACCCCGGGAAAGCATTTTTTCGACGA‑GGGAGTTTTTCAGCTTG‑CCGGGGGTCAGC‑A‑GGGGTTTTGTTTCCTTGTAGGGATACGTTTTTTTGGCCGAGGCGCCGCCctt > 1:52740/1‑148 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgc < 1:986539/151‑1 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGAAGGTTTTTAGTTTG‑CCGGGGGGCAGC‑A‑GGGCTTTTGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgc > 1:942601/1‑151 (MQ=255) acaAACGCGATATGGTACGGCGGGGCGGGCGCAGGAACCAGGGGACGCCTTTTTTCGTCGA‑GGGAGTTTTTTTGATTGCCCGGGGGTAAG‑‑A‑GGGCGTTTGTTTCCTGGGAGGGATAAGTTTTGTTGGGAGAGCCGCCGCCCTTtgcgac > 2:839544/1‑149 (MQ=255) acaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTTCACGA‑GGGAGGTTTTTAGTTTG‑CCGGGGGTCAGC‑A‑GGGCTTTTGGTTTCTGGTAGAGATACGTTTTGTTGGGAGGGGCGCCGCCctttcccacgc > 2:440683/1‑143 (MQ=255) caAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgca < 1:154144/151‑1 (MQ=255) aaCGCGATATGGTACGGCGGGCAGGCTCCAGTACCGAGGGGACGCAATTTCTTGACGC‑GGCAGGTTTTTAGGTTG‑CCGGGGGTCAGCAA‑GGGCTTTGGTTT‑CTGGCAGAGATACGGTTTTTTTGCAGGGGCGCCGCCCTTcgcaac > 2:261684/1‑146 (MQ=255) cgcgATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGGAGTTTTTCAGTTTG‑CCGGGGGGCAGC‑A‑GGGGTTTTGTTTTCTGGGAGAGGTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAgg > 2:372485/1‑151 (MQ=255) atGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAttt < 1:1021836/150‑1 (MQ=255) gTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGGTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTAc > 1:675787/1‑151 (MQ=255) cggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTc > 2:948269/1‑150 (MQ=255) gAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg > 2:850337/1‑151 (MQ=255) | ACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTG‑CCGGGAGTCAGC‑A‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGG > NC_000913/4346528‑4346820 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |